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MATHEMATICAL STUDY OF THE BETAPLANE MODEL: 
EQUATORIAL WAVES AND CONVERGENCE RESULTS 



Isabelle Gallagher, Laure Saint-Raymond 



Abstract. — We are interested in a model of rotating fluids, describing the motion of the ocean 
in the equatorial zone. This model is known as the Saint- Venant, or shallow-water type system, to 
which a rotation term is added whose amplitude is linear with respect to the latitude; in particular 
it vanishes at the equator. After a physical introduction to the model, we describe the various waves 
involved and study in detail the resonances associated with those waves. We then exhibit the formal 
limit system (as the rotation becomes large), obtained as usual by filtering out the waves, and prove 
its wellposedness. Finally we prove three types of convergence results: a weak convergence result 
towards a linear, geostrophic equation, a strong convergence result of the filtered solutions towards 
the unique strong solution to the limit system, and finally a "hybrid" strong convergence result of the 
filtered solutions towards a weak solution to the limit system. In particular we obtain that there are 
no confined equatorial waves in the mean motion as the rotation becomes large. 

Resume. — On s'interesse a un modcle dc fluides en rotation rapide, decrivant le mouvement 
de I'ocean dans la zone equatorialc. Ce modcle est connu sous Ic nom de Saint- Venant, ou systeme 
"shallow water" , auquel on ajoute un tcrme de rotation dont I'amplitudc est lineaire en la latitude ; en 
particulier il s'annule a I'equatcur. Apres une introduction physique au modele, on decrit les differentes 
ondes en jeu et Ton etudie en detail les resonances associees a ces ondes. On exhibe ensuite un systeme 
limite formel (dans la limite d'unc forte rotation), obtenu comme d'habitude en filtrant les ondes, 
et Ton dcmontrc qu'il est bien pose. Enfin on dcmontrc trois types de resultats de convergence : un 
theoreme de convergence faible vers un systeme geostrophique lineaire, un theoreme de convergence 
forte des solutions filtrees vers la solution unique du systeme limite, et enfin un resultat "hybride" 
de convergence forte des solutions filtrees vers une solution faible du systeme limite. En particulier 
on demontre I'absence d'ondes cquatoriales confinees dans le mouvement nioyen, quand la rotation 
augmente. 
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CHAPTER 1 



INTRODUCTION 



The aim of this paper is to obtain a description of geophysical flows, especially oceanic flows, in 
the equatorial zone. For the scales considered, i.e., on domains extending over many thousands of 
kilometers, the forces with dominating influence are the gravity and the Coriolis force. The question 
is therefore to understand how they counterbalance eachother to impose the so-called geostrophic 
constraint on the mean motion, and to describe the oscillations which are generated around this 
geostrophic equilibrium. 

At mid-latitudes, on "small" geographical zones, the variations of the Coriolis force due to the curvature 
of the Earth are usually neglected, which leads to a singular perturbation problem with constant 
coefficients. The corresponding asymptotics, called asymptotics of rotating fluids, have been studied 
by a number of authors. We refer for instance to the review by R. Temam and M. Ziane [34], or to 
the work by J.-Y. Chemin, B. Desjardins, I. Gallagher and E. Grenier [4]. 

In order to a get a more realistic description, which allows for instance to exhibit the speciflcity of the 
equatorial zone, one has to study more intricate models, taking into account especially the interaction 
between the fluid and the atmosphere (free surface), and the geometry of the Earth (variations of 
the local vertical component of the Earth rotation). The mathematical modelling of these various 
phenomena, as well as their respective importance according to the scales considered, have been studied 
in a rather systematic way by A. Majda [24], and R. Klein and A. Majda [19]. 



Here we will focus on quasigeostrophic, oceanic flows, meaning that we will consider horizontal length 
scales of order 1000km and vertical length scales of order 5km, so that the aspect ratio is very small and 
the shallow- water approximation is relevant (see for instance the works by D. Bresch, B. Desjardins 
and C.K. Lin [3] or by J.-F. Gerbeau and B. Perthame [13]). In this framework, the asymptotics of 
homogeneous rotating fluids have been studied by D. Bresch and B. Desjardins [2]. 

For the description of equatorial flows, one has to take further into account the variations of the Coriolis 
force, and especially the fact that it cancels at equator. The inhomogeneity of the Coriolis force has 
already been studied by B. Desjardins and E. Grenier [6] and by the authors [10] for an incompressible 
fluid with rigid lid upper boundary (see also [7] for a study of the wellposedness and weak asymptotics 
of a non viscous model) . The question here is then to understand the combination of the effects due 
to the free surface, and of the effects due to the variations of the Coriolis force. 



CHAPTER 1. INTRODUCTION 



Note that, for the sake of simplicity, we will not discuss the effects of the interaction with the boundaries, 
describing neither the vertical boundary layers, known as Eknian layers (see for instance the paper by 
D. Gerard- Varet [12]), nor the lateral boundary layers, known as Munk and Stommel layers (see for 
instance [6]). We will indeed consider a purely horizontal model, assuming periodicity with respect 
to the longitude (omitting the stopping conditions on the continents) and and infinite domain for the 
latitude (using the exponential decay of the equatorial waves to neglect the boundary) . 



1.1. Physical phenomena observed in the equatorial zone of the earth 

The rotation of the earth has a dominating influence on the way the atmosphere and the ocean respond 
to imposed changes. The dynamic effect is caused (see [14], [28]) by the Coriolis acceleration, which 
is equal to the product of the Coriolis parameter / and the horizontal velocity. 

An important feature of the response of a rotating fluid to gravity is that it does not adjust to a state of 
rest, but rather to an equilibrium which contains more potential energy than does the rest state. The 
steady equilibrium solution is a geostrophic balance, i.e. a balance between the Coriolis acceleration 
and the pressure gradient divided by density. The equation determining this steady solution contains a 
length scale a, called the Rossby radius of deformation, which is equal to c/\f\ where c is the wave speed 
in the absence of rotation effects. If / tends to zero, then a tends to inflnity, indicating that for length 
scales small compared with a, rotation effects are small, whereas for scales comparable to or larger 
than a, rotation effects are important. Added to that mean, geostrophic motion, are time oscillations 
which correspond to the so-called ageostrophic motion. The use of a constant-/ approximation to 
describe motion on the earth is adequate to handle the adjustment process at mid-latitudes : Kelvin 
[35] stated that his wave solutions (also known as Poincare waves) are applicable "in any narrow lake 
or portion of the sea covering not more than a few degrees of the earth's surface, if for ^f we take the 
component of the earth 's angular velocity round a vertical through the locality, that is to say 

-f = nsiiKp, 

where il denotes the earth 's angular velocity and (j> the latitude. " 

The adjustment processes are somewhat special when the Coriolis acceleration vanishes : the equatorial 
zone is actually found to be a waveguide: as explained in [14], there is an equatorial Kelvin wave, and 
there are equatorially trapped waves, which are the equivalent of the Poincare waves in a uniformly 
rotating system. There is also an important new class of waves with much slower frequencies, called 
planetary or quasi-geostrophic waves. These owe their existence to the variations in the undisturbed 
potential vorticity and thus exist at all latitudes. However, the ray paths along which they propagate 
bend, as do the paths of gravity waves, because of the variation of Coriolis parameter with latitude, 
and it is this bending that tends to confine the waves to the equatorial waveguide. 



1.2. A mathematical model for the ocean in the equatorial zone 

In order to explore the qualitative features of the equatorial fiow, we restrict our attention here to a very 
simplified model of oceanography. More precisely, we consider the ocean as an incompressible viscous 
fluid with free surface submitted to gravitation, and further make the following classical assumptions : 

(HI) the density of the fluid is homogeneous. 
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{H2) the pressure law is given by the hydrostatic approximation, 

{H3) the motion is essentiahy horizontal and does not depend on the vertical coordinate, 

leading to the so-called shallow water approximation. 

We therefore consider a so-called viscous Saint- Venant model, which describes vertically averaged flows 
in three dimensional shallow domains in terms of the horizontal mean velocity field u and the depth 
variation h due to the free surface. Taking into account the Coriolis force, a particular model reads as 

dth + V-{hu) ==0 

(1.2.1) , 1 

dt{hu) + W-{hu(g)u) + f{hu)^ + ^hVh-hVK{h)-Aih,u)^0 

Fr 

where / denotes the vertical component of the earth rotation, Fr the Froude number, and K and A 

are the capillarity and viscosity operators. We have written w^ for the vector (u2, —ui). 

Note that, from a theoretical point of view, it is not clear that the use of the shallow water approxima- 
tion is relevant in this context since the Coriolis force is known to generate vertical oscillations which 
are completely neglected in such an approach. Nevertheless, this very simplified model is commonly 
used by physicists [14, 29] and we will see that its study already gives a description of the horizontal 
motion corresponding to experimental observations. 

Of course, in order that the curvature of the earth can be neglected, and that latitude and longitude can 
be considered as cartesian coordinates, we should consider only a thin strip around the equator. This 
means that we should study (1.2.1) on a bounded domain, and supplement it with boundary conditions. 
Nevertheless, as we expect the Coriolis force to confine equatorial waves, we will perform our study 
on R X T where T is the one-dimensional torus 'R./2ttZ, and check a posteriori that oscillating modes 
vanish far from the equator, so that it is reasonable to conjecture that they should not be disturbed 
by boundary conditions. 



1.3. Some orders of magnitude in the equatorial zone 

For motions near the equator, the approximations 

siiKJ) ^ (j), coscf) ~ 1 

may be used, giving what is called the equatorial betaplane approximation. Half of the earth's surface 
lies at latitudes of less than 30° and the maximum percentage error in the above approximation in 
that range of latitudes is only 14%. In this approximation, / is given by 

f = (3x1, 

where xi is distance northward from the equator, taking values in the range 

xi e [-3000 km, 3000 km], 

and /3 is a constant given by 

/3^^ = 2.3xl0-"m-is-i. 
r 

A formal analysis of the linearized versions of the equations shows then that rotation effects do not 
allow the motion in each plane xi — const to be independent because a geostrophic balance between 
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the eastward velocity and the north-south pressure gradient is required. Equatorial waves actually 
decay in a distance of order a^, the so-called equatorial radius of deformation, 

/ \ 1/2 

where c is the square root of gH , H being interpreted as the equivalent depth. For baroclinic ocean 
waves, appropriate values of c are typically in the range 0.5nis~^ to 3nis~^ , so the order of the equatorial 
Rossby radius is 

Qe ^ 100 km, 

which is effectively very small compared with the range of validity of the betaplane approximation. 



1.4. The Cauchy problem for the betaplane model 

Before describing the equatorial waves and the asymptotic behaviour of the ocean in the fast rotation 
limit, we need to give the mathematical framework for our study, and therefore to specify the dissipative 
operators A and K oceuring in (1.2.1). 

From a physical point of view, it would be relevant to model the viscous effects by the following 
operator 

A{h,u) ^ i/V- (hVu), 

meaning in particular that the viscosity cancels when h vanishes. Then, in order for the Cauchy 
problem to be globally well-posed, it is necessary to get some control on the cavitation. Results 
by Bresch and Desjardins [2] show that capillary or friction effects can prevent the formation of 
singularities in the Saint- Venant system (without Coriolis force). On the other hand, in the absence 
of such dissipative effects, Mellet and Vasseur [26] have proved the weak stability of this same system 
under a suitable integrability assumption on the initial velocity field. All these results are based on a 
new entropy inequality [2] which controls in particular the first derivative of v/i. In particular, they 
cannot be easily extended to (1.2.1) since the betaplane approximation of the Coriolis force prevents 
from deriving such an entropy inequality. 

For the sake of simplicity, as we are interested in some asymptotic regime where the depth h is just a 
fluctuation around a mean value H, we will thus consider the following viscosity operator 

A(h, u) — vAu, 

and we will neglect the capillarity 

K{h) = 0, 
so that the usual theory of the isentropic Navier-Stokes equations can be applied (see for instance [23]). 

Theorem 1 (Existence of weak solutions). — Let [h^^vP) he some measurable nonnegative func- 
tion and vector-field onH x T such that 

(1-4.1) .»S//(fil^ + ^|..|A ..<+=.. 



2Ft^ 2 

Then there exists a global weak solution to (1.2.1) satisfying the initial condition 

h\t=o = h°, "|t=o = w", 
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and which furthermore satisfies for almost every t > the energy estimate 



(1.4.2) I [^ -^ + -\u\^]{t,x)dx + i^ I I \Wu\^{t',x)dxdt' <£° 



In this paper we are interested in describing the behaviour of the ocean in the equatorial zone. We 
thus expect the Froude number Fr, which is the ratio of the fluid speed to a measure of the internal 
wave speed, to be small. More precisely we will consider depth variations 

h^H{l+er]) 

where e stands for the order of magnitude of the Froude number. 

As seen in the introduction, in order for gravity waves to be notably modified by rotation effects, the 
Rossby radius of deformation has to be comparable to the typical length scales. In order to derive the 
quasi-geostrophic equations with free-surface term used in oceanography, we will assume that e is also 
the order of magnitude of the Rossby number. 

In non-dimensional variables, the viscous Saint- Venant system (1.2.1) can therefore be rewritten (nor- 
malizing H to H = 1 for simplicity) 



dtV + -V- ((1 + £77)1*) == 0, 

dt({l+er])u] +V ■ ({l + er])u'g)u] + —^{l+er])u^ + -{1 + e'q)Vr] - i^Au ^ 0, 

(1.4.3) V\t=o^v'', "|t=o = "°- 

In such a framework, the energy inequality (1.4.2) provides uniform bounds on any family (?7e,Ug)e>o 
of weak solutions of (1.4.3). 

In all the sequel we will denote respectively -ff " and -ff * the homogeneous and inhomogcneous Sobolev 
spaces of order s, deflned by 

H^R X T) = {/ e 5'(R XT)/ Tf e lI,{R x T) 

and ii/iiif.(RxT) ^J2jj^' + fc'n-^/(^' ^)i' d^ < ^}' 



fcez"^^ 



and 



H^{K X T) = I / e 5'(R XT)/ ||/|||,.(r,t) = E ^ I^ + ?' + ^'11-^/(^,^)1' d( < ^\ 



where J-' denotes the Fourier transform 



Vfc e Z, Ve e R J'fitk) = f e-'''''^e-'^'''f{xi,X2) dx2dxi. 

We will also denote, for all subsets 17 of R x T and for all s > 0, by HQ{il), the closure of 15(17) for 
the H'^ norm, and by H~^{i^) its dual space. 

The following result is a consequence of Theorem 1 . 
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Corollary I.4.I. — Let {r)",u") G ^^(R x T) and (r/0,u^) such that 
(1.4.4) 2 J 



Then, for all e > 0, System (1.4--3) has at least one weak solution {r]^,u^) with initial data (77", u°) 
satisfying the uniform bound 

(1.4.5) ^ firj^^ + {l + eri,)\u,f)it,x)dx+ f j v\Vu,\''{s,x)dxds < E'\ 

Furthermore Ue is uniformly bounded in _L^q^(R+; iy^(R x T)). 

In particular, there exist rj G L°° (11+ ; L'^ {R x T)) and u e L^ {R+ : L'^ (K x T)) D L"^ (11+ ; H^ (R x T)) 
such that, up to extraction of a subsequence, 

{r]e,Ue) -^ (?7, u) in w-L^i^^(R+ x R x T). 

Proof. — Replacing /i by 1 + £77 in the energy inequality (1.4.2), we get (1.4.5) from which we deduce 
that there exist 77 e L°^ {R+ ; L"^ (R x T)) and u e L'^(R+;H^(R x T)) such that, up to extraction of 
a subsequence, 

iVe.u,) - (77, u) in w-L2„,(R+ X R X T). 
Furthermore we have the foUowing inequality: 

\ug\^{t,x)dx < / (1 +e?7g)|wgp(t, a;)dx + e ( / ril(t,x)dx\ ( / \us.\'^(t,x)dx 



In Jn 

from which we deduce that 



u,eLUR+;L^(RxT)), 



where we have used the interpolation inequality 

\u,\*{t,x)dx <C / \u^\'^(t,x)dx / \Vu^\'^(t,x)dx. 
iQ. Jn Jn 

By Fatou's lemma we get that u belongs to i°°(R+; L^(R x T)). That concludes the proof. D 
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EQUATORIAL WAVES 



The aim of this chapter is to describe precisely the various waves induced by the singular perturbation 

(2.0.1) L : (?7, u) e L2(R X T) H^ (V- m, I3xiu^ + Vry). 

In the first paragraph we study the kernel of the operator, which describes the mean flow as we will 
see in Chapter 4. In the second paragraph we describe all the other waves, using the Hermite functions 
in xi and the Fourier transform in X2', this enables us to recover results which are well-known from 
physicists (see for instance [14], [28], [30]- [32], as well as [7] for a mathematical study). Finally in the 
last paragraph we study the possible resonances between all those waves; that result will be useful in 
Chapter 3 to prove regularity estimates for the limit system introduced in Paragraph 2.3 below. 



2.1. The geostrophic constraint 

In this section we are going to study the kernel of the singular perturbation L defined in (2.0.1). 

Proposition 2.1.1. Define the linear operator L by (2.0.1). Then [rj^u) E i^(R x T) belongs 
to KerL if and only if {r],u) belongs to L'^(Rxi) cmd 

(2.1.1) ui = 0, (3xiU2 + dif]^0. 

Proof. — If (ry, u) belongs to L^(R x T) n KerL, then we have 

V- u = 0, f3xiu^ + Vr] = 0, 

in the sense of distributions. Computing the vorticity in the second identity leads to 

V-^ • {f3xiu^ + Vt]) = l3{xiV • u + ui) = 0, 

from which we deduce, since V-u = 0, that ui = 0. Plugging this identity respectively in the divergence- 
free condition and in the second component of the vectorial condition gives 

d2U2 = 0, d2T] = 0, 

meaning that rj and u depend only on the xi variable. The last condition can then be rewritten 

(3xiU2 + dirj = 0. 
Conversely, it is easy to check that any (?7, u) G -^^(R) satisfying (2.1.1) belongs to KerL. D 
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In the following wc will denote by Hq the orthogonal projection of L^(R x T) onto KerL. It is given 
by the following formula. 

Proposition 2.1.2. — Define the linear operator L by (2.0.1). Denote by flp the orthogonal projec- 
tion of L2(R X T) onto KerL. Then, for all (77, u) G ^^(R x T) 

(2.1.2) no(?7,u)= ( f {DD'^ + Id)-\r] + Du2)dx2,0, f D^iDD'^ + Id)-\r] + Du2)dx2] , 



where D is the differential operator defined by D- = di — — 

\pxi 

Proof — By Proposition 2.1.1, for all {t],u) e L^(R x T), {ri*,u*) = no(?7,M) belongs to L^(R) and 
satisfies 

< = 0, I3xiu*2+dif]* =Q. 

Averaging with respect to the X2-variable, one is reduced to the case when (77, u) G i^(R). 

By definition {rj — rj* ,u — u*) is orthogonal in L^ to any element (p, v) of KerL : that implies that 

/ ((j? - V*)p + ("2 - u*)v2) dxi^ Ur]- r]*)p - {u2 + l^^^'^l^^^p) ^-^i = 0. 
An integration by parts leads then to 

1 ^ , r ,\ * . ^ U2 



^i75F^'^i +Id]r]* ^ri + di- — 
p'^x^ ) \ix\ 

Plugging this identity in the second constraint equation gives the expected formula for u\ . 

That proves Proposition 2.1.2. D 

2.2. Description of the waves 

In this section we are going to describe precisely the various waves created by L. In the first paragraph 
of this section (Paragraph 2.2.1) wc compute the eigenvalues of L and present its eigenvectors, which 
constitute a Hilbertian basis of L^(R x T) (that is proved in Paragraph 2.2.2). That basis enables us 
in the last paragraph to introduce the filtering operator and formally derive the limit filtered system, 
in the spirit of S. Schochet [33] (see also [17]). 

2.2.1. Precise description of the oscillations. — In this paragraph, we are going to explain how 
to obtain the various eigenmodes of L. The crucial point is that the description of these eigenmodes 
can be achieved using the Fourier transform with respect to X2 and the decomposition on the Hermite 
functions ('!/'n)neN with respect to x\. Here the Hermite functions are conveniently rescaled so that 

_M r- 

■^n{xx) = e 2 P„(^/3a;i), 

where P„ is a polynomial of degree n, and (V'n)neN satisfy 

-C + /3'a;?Vn = /3(2n+l)V'„. 

We recall that (V'Ti)neN constitutes a Hermitian basis of L^(R). 
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Moreover we have the identities 
(2.2.1) 



i^'aixi) - f3xiip.a{xi) = -A/2^(rr+T)V'n+i(a::i)- 
We have used the convention that ■(/;„ = if n < 0. 

In the foUowing we will then denote by f{n, fe), for n G N and fc € Z, the components of any function / 
in the Hermite-Fouricr basis {2T:)^^^'^ijJn{xi)e^'"^'^ . In other words we have 

V(n, fc)eNxZ, f{n,k) = —^l i;.a{xi)e'''''''^ f{xi,X2) dxidx2, 



along with the inversion formula 

y{xi,X2)eIixT, f(x,,X2)^^ Y. Mxi)e"''''f{n,k). 



27r ^, 
feez 

In order to investigate the spectrum of L (which is an unbounded skew-symmetric operator), we are 
interested in the non trivial solutions to 

(2.2.2) L{7j,u)^it{7^,u). 

If one looks for the L^ solutions of (2.2.2) with ui non identically zero, one gets as a necessary condition 
that the Fourier transform of ui with respect to X2 (denoted by ^^2^1) satisfies 

(^2«i)" + ("r^ - fc2 + ^ „ p^A T2U, = 0, 

from which we deduce that JF2W1 is proportional to some ipn and that 

(2.2.3) T^ - {k^ + I3{2n + 1))t + f3k ^ 0, 
for some n G N. 

The following lemma is proved by elementary algebraic computations. 

Lemma 2.2.1. — For any /3 > and any {n,k) G N* x Z, the polynomial 

P{t) = T^ - (fc2 + f3{2n + 1))t + (3k 
has three distinct roots in R, denoted in the following way: 

(2.2.4) T(n, fc, -1) < T(n, k, 0) < T(n, k, 1). 

Moreover if T{n,k,j) — T{n' ,k,j') ^ for some {n,n') G N^ with n ^ and {j,j') G { — 1,0,1}^, then 
necessarily n ~ n' and j = j' . 

Finally the following asymptotics hold if n or \k\ goes to infinity: 

Tin,k,±l)^±^k^+P{2n+1), and r(n,fc,0)- k^ + mn+1) ' 

Proof. — To prove that the polynomial has three distinct roots we simply analyze the function P{t). 
Its derivative P'{t) vanishes in ±a, where 



k^ +f3{2n+l) 



10 
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It is then enough to prove that P{a) < and P{—a) > 0. Let us write the argument for P{a). We 
have 

P{a) = -2a^ + 13k. 
But for n 7^ 0, one checks easily that 2a'^ > /3|fc|. Indeed one has 

4a6 = ^(fc' + I3{2n + l)f > P/32 

as soon as n > 1. So the first rcsuh of the Icmnia is proved. 

To prove the second resuh, we notice that if T{n,k,j) — T{n',k,j') = t 7^ 0, then 2(n — n')T ~ 
from which we deduce that n = n' , and therefore that j ~ j' since the polynomial (2.2.3) admits three 
separate roots for n ^ 0. 



Finally the asymptotics of the eigenvalues is an easy computation. The lemma is proved. 



D 



Remark 2.2.2. — In the case when n~Q, the three roots of P are 

1 



k 1 



(2.2.5) T(0,fc,-l) = ----v/fc2+4/3, r(0,fc,0)==fc, and t(0, k,l) ^ -- + -^Jk"^ + 4/3. 

It follows that in the case when j3 — 2k^ , the roots become k (double) and —2k. 

Now let us study more precisely the waves generated by L. 



• If fc 7^ and n y^ 0, (2.2.3) admits three solutions according to Lemma 2.2.1, and one can check (see 
Paragraph 2.2.2 below) that these solutions are eigenvalues of L associated to the following unitary 
eigenvectors 



/ iT{n,k,j) , 



*n,fcj(a;i,a;2) = Cn,k,je 



ikX2 



ik 



k^ — T{n, k, jy 

ik 
\ k'^ — T{n, k,j)'^ 



k^ — T{n, k, j)^ 



Pxi1pn{xi) 



\ 



Vn a^l - 71 7 p-T^PXlV'n Xl) 

kr — T[n, fc, j)^ 



which can be rewritten 
(2.2.6) 



'^n,k,jixi,X2) = Cn^k^e 



ikx2 



T{n,k,j) + k\l 2 



— -(/'„-i(a;i) + 



— -0„_i(xi) + 



T{n,k,i) - k 

i Il3{n + 1 



/3(» + l) , f , \ 

Ipn+iixi) 



■rpn+i{xi) 



V Tin,k,j) + k\ 2 "" ^^ '' r(n,fc,j)-fcV 2 

because of the identities (2.2.1). The factor Cn.k.j ensures that '^n,k,j is of norm 1 in L^(R x T), its 
precise value is given in (2.2.12) below. 

The modes corresponding to T(n, fc, —1) and T(n, k, 1) are called Poincare modes, and satisfy 

T{n, k, ±1) -- ±^fc2 + P{2n+ 1) as |fc| or n ^ 00, 
which are the frequencies of the gravity waves. 

The modes corresponding to T(n, k, 0) are called Rossby modes, and satisfy 

I3k 



T{n, fc,0) 



fc2+/3(2n+l) 



as |fc| orn ^- 00, 



2.2. DESCRIPTION OF THE WAVES 



11 



meaning that the oscillation frequency is very small : the planetary waves ^n,fc,o satisfy indeed the 
quasi-geostrophic approximation. 



If fc = and n 7^ 0, the three distinct solutions to (2.2.3) are the two Poincarc modes 



r(n,0,±l) = ±v//3(2n + l) 

and the non-oscillating, or geostrophic, mode T(n, 0,0) ~ 0. The corresponding eigenvectors of L are 
given by (2.2.6) if j 7^ and by 



(2.2.7) 



*n, 0.0(2:1) 



^27r(2n+l) 







— ^■0„_l(xi) - J-lpn+lixi) 



• If n = 0, the three solutions to (2.2.3) are the two Poincarc and mixed Poincare-Rossby modes 
(2.2.8) 



-(0,fc,±l) = -^±iv/FT4^ 



with asymptotic behaviours given by 

t(0, k, — sgn(fc)) ^ — fc as |fc| -^ 00, 

t(0, /c, sgn(fc)) ^ — as \k\ -^ 00, 
k 

and the Kelvin mode r(0, fc, O) — k. The corresponding eigenvectors of L are given by (2.2.6) if j ^ 
and by 



(2.2.9) 



*o.fc.o(a;i,a;2) 



1 



Akx2 



Att 



Mxi) 


V'o(a;i) 



Note that in the case when the fluid studied is the atmosphere rather than the ocean, the mixed 
Poincare-Rossby waves are known as Yanai waves. 

We recall that the functions ■0„ are defined by 

'ipnixi) = e ~Pn{\/Pxi), 

where P„ is a polynomial of degree n. We therefore have an exponential decay far from the equator. 

As mentioned in the introduction, the adjustment processes are somewhat special in the vicinity of 
the equator (when the Coriolis acceleration vanishes). A very important property of the equatorial 
zone is that it acts as a waveguide, i.e., disturbances are trapped in the vicinity of the equator. The 
waveguide effect is due entirely to the variation of Coriolis parameter with latitude. 

Note that another important effect of the waveguide is the separation into a discrete set of modes 
n = 0, 1, 2, ... as occurs in a channel. 
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2.2.2. Diagonalization of the singular perturbation. — In this paragraph we are going to show 
that the previous study does provide a Hilbertian basis of eigenvectors. 

Proposition 2.2.3. — For all (n,k,j) e N x Z x {—1,0,1}, denote by T{n,k,j) the three roots 
of (2.2.3) and by '^n,k,j the unitary vector defined in Paragraph 2.2.1. 

Then (^n.fe,j)(n,fc,j)eNxZx{-i,o.i} ** O' Hilbertian basis o/_L^(Rx T) constituted of eigenvectors of L : 
(2.2.10) V(n,fc,j)eNxZx {-1,0,1}, L*„,fc,, == ir(n, fc, j)*„,fc,j. 



Furthermore we have the following estimates : for all s > 0, there exists a nonnegative constant Cs 
such that, for all (n, fc, j) e N x Z x {—1, 0, 1}, 

||^«,fej|U-(RxT) <Co, ||*„,fcj||M/».oo(RxT) <C,(l + |fc|2 + n)''/^ 

C-\l + |fc|2 + n)*/2 < ||*„,fe,,||Ha(RxT) < C,(l + |fc|2 + nY'\ 

where W^'°° denotes the usual Sobolev space. Moreover the eigenspace associated with any non zero 
eigenvalue is of finite dimension. 



Proof. — In order to estabhsh the diagonahzation resuh, the three points to be checked are the 
identity (2.2.10), the orthonormaUty of the family (^„,fcj), and the fact that it generates the whole 
space i2(Rx T). 

• ^n.k.j is an eigenvector of L 

We start by establishing the identity (2.2.10), where T(n,k,j) is defined by (2.2.4) and (2.2.5) and 
^n.fe.j is defined either by (2.2.6) (for the Poincarc and Rossby modes) or by (2.2.7) (for the non- 
oscillating modes), or by (2.2.9) (for the Kelvin modes). 

For the Poincarc, Rossby and mixed Poincarc-Rossby modes, we start from formula (2.2.6) 

f -i f]3^ 



^n,k,j — ^n,k,j^ 



ikx2 



1pn+l(Xl) 



—I pn I 

T{n,k,j) + k\ 2 " T{n,k,j) — k 

tpn{xi) 

i I /3n , , , i I P{n + 1) , , , 

V r(n,fc,,) + fcVT^"-^("^)+r(n,fc,,)-fcV^^— '^"+^^"^) 



We have L'i'n.kj — Cn,k,je^''^^Vn^kj where Vn,k,j denotes 



i^'ni^l) + ^kJ 7- 



ix/n zVnTT 

Wn-l{Xl) + —. ; rV'n+l(a;i) 



2 \T{n,k,j) + k 



T{n, k,j) — k 

-{Pxi1p.a^l{xi) - -Ip'n^iixi)) + — 7— -{(3xi1pn+l{xi) + 1p'^^^{xi)) 



V 



r(n, fc, j) + fc 

-/3xiV'n(a;i) +ik\j ^ 



T{n,k,j) - fc 



2 \T{n,k,j) + k 
which can be rewritten using the identities (2.2.1) 

/ T{n,k,j) 



, , , t^/n + 1 

i)n-l[Xl) + — J— -1p„+l{Xi) 

r{n,k,j) - k 



l]3n T{n,k,j) f3{n + 1) \ 

r(n,fc,,)+fcVT^"-^("^)' r(n,fc,,)-fc V^^^"+^^"^^ 



r{n,k,j) + k 
'r{n,k,j) fpn 
\ T{n,k,j) + k 



PnTpn{xi) + 



T{n,k,j) 



-/3(n + l)V'„(xi) 



— -0„-i(a;i) 



T{n,k,j) Jl3{n + 1) 
T{n, k,j) — fc 



-'(/;„+i(xi) 
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As T{n,k,j) satisfies (2.2.3), we have 
i „ i 



-(3n- 



T{n,k,j)+k T{n,k,j)-k 

from which we deduce that 



I3{n + 1 = I — j—-^ r^ = iT{n, fc, J 



L^nM.] =«T(n, fc,j)*„^fej for all in,k,j) G N x Z x {-1,1} UN* x Z* x {0}. 



For the Kelvin modes we start from formula (2.2.9) 

1 



* 



O.fc.O 



47r 



■00 

ikx2 I Q 
■00 



We have 



or equivalently 



i*o,fc,o 



Akx2 



ikxpQ 



ik 



Att 



^tkX2 I Q 



ikipo 
L^o,k,Q = 'ik^o,k,Q for all fc G Z 



■00 





For the non-oscillating modes we start from formula (2.2.7) 



*n,0,0 



A/27r(2n + l) 
An easy computation shows that 







/ 



n+1 n 

-^—Vn-l - \Jl^Vn+l 







L^nflfi = 



^27r(2n + l) 
which is zero by (2.2.1). Thus, 



^ (XiV-n-l -1pn-l)- \^ ^{Xilpn+l + -ipn+l) 





L'^nfifi ^ for all n e N* 



• (^n.fe.j) is an orthonormal family 

By identity (2.2.10) and the fact that {{2T:)^^^'^e^'"^^)kez and (V'n(a;i))neN are respectively Hilbertian 
basis of L^(T) and i^(R) we are going to deduce that ('^n,k,j) is an orthonormal family. 

In formula (2.2.6), we choose 



(2.2.12) 
so that 



CnMj = (2^)-l/2 f 



f3n 



/3(n + l) 



(T(n, fc,j) + fc)2 {T{n,kJ)-k)^ 



-1/2 



ll*n,fe,illL2(RxT) — Ij 

for aU in,k,j) G Nx Z x {-1,1} UN* x Z* x {0}. 
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In the same way, it is immediate to check that 

ll*n,0,ollL2(RxT) = 1' 

for all n E N*, and that 

ll*0,fe,o||L2(RxT) = 1' 

for aU fc e Z. 

In order to establish the orthogonality property we proceed in two steps. 

If T(ri, k,j) 7^ T(ri', k',j'), as L is a skew-symmetric operator, we have 

iT(n,fc, j)(*„,fej|*„-,fe-y) = -(i*„,fcj|*„',fe'j') 

= (*„,fcj|L*„-,fc/j-) 
= iT(n',fc',j')(*«,fcj|*«',fe'j') 
from which we deduce that 

(*n,fcj|*«%fc'j') =0. 

If r(n, k,j) — T{n' , k',j'), we first note that 

iork^k', (*„,fej|*„',fc'j') =0 

using the orthogonality of e''^^^ ande**^^^. So we are left with the case when fc = fc'. First, if T(n, fc,j) = 
T{n',k,j') = T 7^ and n ^ 0, then Lemma 2.2.1 implies that n = n' and j = j' . Then in the case 
when T{0,k,j) = t(0, fc,j') = t 7^ 0, with j ^ j' , we just have to consider the explicit definition 
of '^o,k,j and ^o,fcj' given in Paragraph 2.2.1 to find that 

Finally, if T{n,k,j) = T{n',k',j') ~ 0, we have k = k' = and j ~ j' = and we deduce from 
formula (2.2.7) that 

for n ^ n', (*„,o,o|*«%o,o) - 0. 



We thus conclude that 



as soon as {n,k,j) ^ {n',k',j'). 



(*«,fcj|*«%fe'j') =0, 



• (*«,fcj) spans L2(Rx T) 

It remains therefore to see that any vector $ of L^(R x T) can be decomposed on the family {'i'n.k.j) 

We first decompose each component on the Hermite-Fourier basis 

/ $o(fc,«)^n(a;i) 

^ ^ n,k \ ^2{k,n)ijn{xi) 
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which can be rewritten 

($o(fc, n + 1) + *2(fc, n + l))V'„+i(xi) + ($o(fc, n-l)- $2(fc, n - l))V'«-i(a;i) 

V e'^''^ I 
2\/27r 



($o(fc, n + 1) + $2(fc, J^ + l))^n+i(a;i) - ($o(fc, J^ - 1) - $2(fc, n - l))V'„-i(a;i) 
/ ($o(fc, 0) + $2(fc, O))^o(xi) + ($o(fc, 1) + $2(fc, l))^i(xi) 



n>0 
k 



T^T.' 



2V2tt 



2$i(fc,0)7Ao(xi) 
y ($o(fc, 0) + $2(fc, O))^o(xi) + ($o(fc, 1) + $2(fc, l))^i(xi) 



We then introduce for all (n, fc) e N x Z the matrix Mn.k G ^3(R-) defined by 

f -iCn^k.-i\/l3n/2 -iCn,kfi^/(in/2 -iCn,k,i^/ (in/2 \ 

T(n, fc,— l) + fc T(n, fc,0)+fc T(n, fc,l) + fc 

iCn,k,-iyJP{n + l)/2 iCn,kfi^P{n + l)/2 iCn,k,i^/ P{n + l)/2 

T{n, k,l) — k / 



(2.2.13) M„,fc 
if n 7^ and k ^ 0,by 

(2.2.14) M„,o 

if n 7^ and by 

(2.2.15) 



\ T{n,k,—1) — k T{n,k,0) — k 



r(n,0,-l) 

C'n,0,-1 

zC„,o,-iV/3(n+l)/2 
V r(n,0,-l) 



T(n, 0,1) 

C'n. 0,1 



-C„,o,oV^^ 



T(n, 0,1) / 



Mo,fc = 



/ Vl/47r \ 

Ccfc,-! Co,fc,l 

iCo.k,^i\/P/2 iCQ.k,i\/ P/2 

V r(0,fc,-l)-fc T(0,/c,l)-fc / 



As the eigenvectors ^n,fc,-i, '^n,kfi and ^„,fc,i are orthogonal in L^(R x T), these matrices are neces- 
sarily invertible. 



We conclude by checking that one can write 

n,k,j 

where fn,k,j is defined by 



¥'ri,fe,-l 
Vn,kfi 
'Pn,k,l 



1 



/2^ 



■^^-^n,k 



( i($o(fc,n-l)-$2(fc,n-l)) \ 

$i(fc,n) 
\^ -($o(fc,n+l) + d2(fc,n + l)) y 



for n ^ 0, and by 



</50,fe,0 
</50,fc,l 



/2^ 



^'-'o.fe 



( ^(<l>o(fc,0) + <l>2(fc,0)) \ 



1 



$i(fc,0) 



y -($o(fc,i) + $2(fc,i)) y 
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• The regularity estimates are obtained using the exphcit formulas (2.2.6), (2.2.7) and (2.2.9), as well 
as the following bounds on the elementary Fourier and Hermite functions : 



^^'"^'li^=(T) = Ifcnie^'^llL^(T), l|e^'^1|^„»(T) = |A:n|e^'=-^|U»(T) 



e- 



and 

ll'0n||//=(R) ^ (1+")''/^SUp||-0„||l2(r,), ||-0n||H'=.-(R) < CsC^ + ")''/^ SUp || V'™ || L- (R) 

n n 

coming from identities (2.2.1) by a simple recurrence. The crucial point is therefore to have a uniform 
i°°-bound on the Hermite functions, which is stated for instance in [20]: 

(2.2.16) VneN, ||^/>„||ioo(R,) < Coo with Coo -- 1.086435. 

Finally let us prove that the eigenspace associated with a nonzero eigenvalue is of finite dimension. 
Suppose by contradiction that there is A 7^ and a sequence (tt-m, ^m, Jm)mgn in N x Z x {—1, 0, 1} 
such that 

T{nM,kM,JM) ~ A and um + I^mI —> 00, as M —> 00. 

By Lemma 2.2.1, as n or |fc| goes to infinity, the eigenvalue T{n,k,j) goes to zero or to ±00, which 
contradicts the assumption that T{nM,kM,JM) = A. 

This concludes the proof of Proposition 2.2.3. D 



As the behaviour of the eigenmodes are expected to depend strongly of their type, i.e. of the class of 
the corresponding eigenvalue, we split L^(Rx T) into five supplementary subsets, namely the Poincare 
modes, the Rossby modes, the mixed Poincarc-Rossby modes, the Kelvin modes and the non-oscillating 
modes. 

Definition 2.2.4- — With the above notation, let us define 

P = yect{*„,fc,, / (n, k,j) e N* X Z X {-1, 1} U {0} x {(fc, -sign{k))ke7,'}^ {0} x {0} x {-1, 1}}, 

R = Vect{-^nMfi I (n, fc) e N* X Z*}, 
M = Fect{*o,fcj / fc e Z*, j = sign{k)}, 

K = Vect{^o^kfi/keZ*}, 

so that L^{IixT) = P® R<S) M ®K QKerL. Then we denote by Up (resp. IIii,IIm,TIk and Uq) the 
L^ orthogonal projection on P (resp. on R, M , K and KerLj. 

Moreover we define & the set of all eigenvalues of L, as well as the following subsets of IS! x X x 
{-1,0,1}.- 

6p = {rin, fc, j) / {n, fc, j) G N* x Z x {-1, 1}} U {t(0, fc, -sign{k)) / fc e Z*} U {±v^}, 

6i? = {T(n,fc,0) /(n,fc,j)eN* X Z*}, and 6/^ == Z*. 

Finally it can be useful for the rest of the study to sum up the previous notation in the following 
picture. 



2.2. DESCRIPTION OF THE WAVES 



ir 



wave 


n 


k 


j 


definition of *n,fcj 


definition of T(n, k, j) 


Poincare 


N* 
{0} 
{0} 


Z 

z* 

{0} 


{-1,1} 

-sign(A:) 
{-1,1} 


(2.2.6) page 10 
(2.2.6) page 11 
(2.2.8) page 11 




r(n, fc, ±1) - ±v/fc2 + (3{2n + 1) 

r(0,fc,-^sign(fc)) k 

t(0,0,±1)=±v^ 


Mixed 


{0} 


z* 


sign(fc) 


(2.2.6) page 11 


T{0,k,s[gn{k))^'j 


Kelvin 


{0} 


z* 


{0} 


(2.2.9) page 11 


r(0,fc,0) = fc 


Rossby 


N* 


z* 


{0} 


(2.2.6) page 10 


-(■n /" n^i - ^ 


r{n,L,0)^- fc2 + ^(2n+l) 


non oscillating 


N 


{0} 


{0} 


(2.2.7) page 11 


T(n,0,0) = 



Table 1. Description of the waves 



2.2.3. Orthogonality properties of the eigenvectors. — In this section we are going to give 
some additional properties on the ^n,fej defined above, which will be useful in the next chapters. We 
will write Tln.k.j for the projection on the eigenmode "^n.k.j of L, and IIa for the projection on the 
eigenspace associated with the eigenvalue iX of L. The main result is the following, which states an 
orthogonality property for the ageostrophic modes (meaning the eigenvectors in [KerL)^). Note that 
there is no analogue of that result for gcostrophic modes. 

Proposition 2.2.5. — Let s > be a given real number. There is a constant C > such that for 
any non zero eigenvalue i\ of L and for any three component vector field $ in (Keri)^, we have 

2 ^ IITT. <Tvl|2 



(2.2.17) C7I 



r{n,k,j) — X 



<Cs Y^ I 



nn,fc,j$||^,(R^T) 



Proof. — Let $ in (KerL)-'- be given and let s be any integer (the result for all s > will follow by 
interpolation). We have 



r(n.k,j) — X 



E 



(a*(n„,fe,,ci>)|a^(n„.,fe-,,.ci>))^.(j,^T) 



(n,k,j)^(n» ,k',j') 

Of course, 

(2.2.18) (9^*„,fcj|9^*„*,fe.j-)L2(KxT) = if fc ^ fc*. 

Moreover we know by Proposition 2.2.3, page 12 that if T(n,k,j) = T{n* ,k,j*) = A 7^ and n 7^ 0, 

then necessarily n — n* and j = j* . Therefore one has in fact 



E 



T(n,k,j) = ^(n'-,k*,j*) = \ 
(n,fc,j)^(n*,fc*,3*) 



E (5*(no.fe,,<i>)|a^(no,fe,,.$))^.(HxT)' 



T(0,k,j) = T(0,k,j*) = \, 



But according to Remark 2.2.2 page 10, such a situation occurs only if 2k'^ = f3, in which case t(0, k, j) 
is equal to k. So there is at most one possible value for k (k — X) which occurs only in the case 



18 CHAPTER 2. EQUATORIAL WAVES 

when A = -ii^J j3/2. In this last case, we have obviously 

The result follows. D 

Remark 2.2.6. — Note that the same argument allows actually to prove similar estimates for the 
components separately : 

T(n,k,j) — X T{n.k.j)—X 

T{n,k,j) — X T{n,k,j)—\ 

denoting by $0 the first coordinate and by $' the two other coordinates of^. 



2.3. The filtering operator and the formal limit system 

In the previous paragraph wc have presented a Hilbertian basis of L^ (R x T) consisting in eigenvectors 
of the singular penalization L. We are then able to define, in the spirit of S. Schochet [33], the "filtering 
operator" associated with the system. 

Let C be the semi-group generated by L : we write C{t) — exp (—tL). Then, for any three component 
vector field $ G i^(R x T), we have 

(2.3.1) £(i)$ = E e"**^nA$, 

where Il\ denotes the L'^ orthogonal projection on the eigenspace of L corresponding to the eigen- 
value iX, and where & denotes the set of all the eigenvalues of L. 

Now let us consider (j^^^Ue) a weak solution to (1.4.3), which is formally equivalent to 



dtVe 


e 


((1 + e% 


k) = 


-0, 


4 + 


-Viye 


V 


■AWe : 


= 0, 




Vrit^n 


= ^7° u.^ 


Hf=n = 


--ul 



and let us define 

(2.3.2) $, =/:('-^V77„«e). 

Conjugating formally equation (2.3.2) by the semi-group leads to 

(2^3^3) a*. +£ (4) Q [C (i) *.,£ (i) *,) - .£ (-£) A'C (i | *. . R, 

where A' and Q are the linear and symmetric bilinear operator defined by 
(2.3.4) A'$ = (0, A$') and Q($, $) = (V • ($o$')> (*' ' V)$') 



and 
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e) l + e?7e 



We therefore expect to get a bound on the time derivative of <I>e in some space of distributions. A 
formal passage to the hmit in (2.3.3) as e goes to zero (based on formula (2.3.1) and on a nonstationary 
phase argument) leads then to 

(2.3.5) at$ + QL($,$)-j/A'^$ = 0, 

where A'^ and Qj^ denote the linear and symmetric bilinear operator defined by 

(2.3.6) A^$= ^ nAA'nA$andOL($,$)= ^ nAO(n^$,n^$). 

X = fj. + ^i 

The study of (2.3.5) is the object of Chapter 3. The proof that (2.3.5) is indeed the limit system 
to (2.3.3) is the object of Chapter 4. 

In the next section we study the resonances associated with the operator L: more precisely we describe 
in what cases the equality 

T{n,k,j) +T{n*,k*,j*) = T{m,k + k*,£) 

can hold. That will be very important in the rest of the study, to understand the structure of the 
nonlinear terms in (2.3.5). 



2.4. Interactions between equatorial waves 

In this section we will study the nonlinear term in (2.3.5). We will first study the resonances of L, and 
then prove that the projection of (2.3.5) onto the kernel of i is a linear equation. 

2.4.1. Study of the resonances. — Let us prove the following result. 

Proposition 2.4-1- — Except for a countable number of (3 and with the notation of Section 2.2.1, the 
following condition of non resonance holds for all n,n* ,m C^ N, all fc, /c* G Z and all j, j* , ^ G { — 1, 0, 1}; 

T{n,k,j) +T{n*,k*,j*) =T{m,k + k*,i) 

implies 

either T{n, k,j) ^ or T{n* , k* ,j*) ~ or T{m, k + k*,£) ~ 0, 

or T{n,k, j),T{n* , k* , j*),T{m, k + k* ,£) G &k, 

meaning that, among the ageostrophic modes, only three Kelvin waves may interact. 

Proof. — Let us start by noticing that by definition of Kelvin waves, Kelvin resonances necessarily 
take place simply because they correspond to convolution in Fourier space. 

Before starting with technical results, let us describe the main ideas of the proof. The crucial argument 
is that the eigenvalues of the penalization operator L are defined as the roots of a countable number of 
polynomials whose coefficients depend (linearly) on the ratio (3. In particular, for fixed n,n* ,m G N 
and k,k* G Z, the occurence of a resonant triad 

T{n,k,j) +T{n*,k*,j*) ^T{m,k + k*,i) 
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is controlled by the cancellation of some polynomial Pn,n*,m.k.k*{l3)- Therefore, either this polynomial 
has a finite number of zeros, or it is identically zero. The difficulty here is that we are not able to 
eliminate the second possibility using only the asymptotics /3 ^ oo. We therefore also study the 
asymptotics /3 — > 0, and in the case when n = or n* = 0, we have to refine the previous argument 
introducing an auxiliary polynomial. 



• Definition of the polynomial Pn,n*.m,k,k-'{P) 

For fixed n, n*, to e N and k, k* G Z, it is natural to consider the following quantity 

Pu,n',m.k,k'{P)= n {T{n,k,j)+T{n*,k*,j*)~'T{m,k + k*,i)). 

j,j*./e{-i..o,i} 

Considerations of symmetry show that this quantity can be rewritten as a polynomial of the sym- 
metric functions of (T(n, fc, j))jg{_i^o.i}7 the symmetric functions of (r(n*, fc*, j*))j.g{_i^o,i} and the 
symmetric functions of (r(TO,, /c + k* ,£))£^^_iqij. 

Therefore, as the eigenvalues (T(ri, fc, j))je{-i.o,i} of the linear penalization L are defined as the three 
roots of a polynomial (2.2.3) whose coefficients depend (linearly) on (3 

T^ - {k^ + /3(2n + l))r + (3k = 0, 

the symmetric functions of (T(n, fc, j))j£|_i 0,1} satisfy 

JG{-1,04} 

(2.4.7) E X{r{n,k,f)^-{e + {2n+l)P), 

je{-ifl.i}3'i^j 

Jl T{n,k,j) ^-Pk, 
je{-i,o,i} 

from which we deduce that Pn,n*,m.k.k*{P) is a polynomial (of degree at most 27) with respect to /3. 

In particular, for fixed n,n* ,m € N and k, k* € Z, either -Pn,n*,m,fc,fc* (/3) is identically zero or it has a 
finite number of roots. In other words, that means that 

(a) either, for all f3 E R*, there is a resonance of the type 

T{n,k,j) +T{n*,k*,j*) =T{m,k + k*,i) 
for some j,j*,£ E {-1,0,1}, 

(b) or, except for a finite number of /3, such resonances do not occur. 



• Asymptotic behaviour of Pn^n*^m.k.k*{P) as (3 ^ 00 

In order to discard one of these alternatives, we are interested in the asymptotic behaviour of the 
polynomial Pn,n'-,mMM*{l3) as /3 ^ cx). 

We start by describing the asymptotic behaviour of each root {r^n, fc, j))je{-i.o,i} as /3 ^ cx3. 
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Lemma 2.4.2. — With the notation of Paragraph 2.2.1, for all k ^ Z and all n € N, the following 
expansions hold as /3 ^ oo ; 



-K k, 1) = v/(2n+l)/3 - „./^^. + o(l), 

2,{2,n + Ij 



(2.4.8) T(n, fc, -1) = -v/(2n + l)/3 - ^^^^ ^ ^^ + o(l), 



fc 4n(n + l)fc-'' /I 

2nTT " (2n+ 1)4/3 ^ " 1^^ 



^("^fe^O)-:^^^- /ol , .^4^ +0 



Proof. — We start with the most general case, namely the case when k ^ Q. We proceed by successive 
approximations. As the product of the roots — /3fc tends to infinity as /3 ^ oo, there is at least one 
root which tends to infinity. Therefore, we get at leading order 

T3-/3(2n+l)T = 0, 

which implies that the Poincarc and mixed Poincare-Rossby modes are approximately given by 

T{n, fc, ±1) ~ ±V(2n+l)/3. 

Plugging this Ansatz in the formula 

T2 = (2n + l)/3 + fc2-^ = (2n+l)/3('l ^ ' ^' 



T ' " \ (2n+l)r /3(2n+l) 

provides the next order approximation of the Poincarc modes, namely 

fc 



-(n,fc,j)^jV(2«+l)/3- 

e roots is 
Rossby mode, satisfies 



2(2n+l) 
Then, as the sum of the roots is zero (see (2.4.7)), we deduce that the third mode, i.e. the Kelvin or 



r(fc,n,0) = — ^+o(l). 



2n+ 1 
Plugging this Ansatz in the formula 

/3fc + T 



(2n+l)/? + fc2 

leads then to 

fc / 1 fc3 fc2 \ /I 



"("'"'') = 2;^ Iv' + ^ (^^^TIF " (^^^TT)^ J + %^ 



The other case (when fc = 0) is dealt with in a very simple way. The Poincare modes are exactly 

r(n,0,±l) = ±V(2n+l)/3, 

whereas the third mode is zero 

T(n,0,0) = 0, 
and thus they satisfy the general identities (2.4.8). 

The result is proved. D 
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Equipped with this technical lemma, we are now able to characterize the asymptotic behaviour of most 
of the factors 

T{n, k,j) + T{n*,k\f) - rim, k + fc*, £) 

in Pn,n' .mMM* (/?) as /3 ^ oo. 

Lemma 2.4-3. — With the notations of Paragraph 2.2.1, any triad of non zero modes 

{T{n,k,j),T{n*,k*,j*),T{m,k + k*,£)) 

with k,k* € Zi and n,n*,m G N, which is not constituted of three Kelvin or three Rossby modes, is 
asymptotically non resonant as P ^ oo. 

More precisely the following expansions hold as f3 —^ oo : 

(i) for three Poincare or mixed Poincare-Rossby modes (j 7^ and j* 7^ and £ j^ 0) 

T{n,kJ) + T{n*,k*,j*)-T{m,k + k*,£) ~ ^/S (jV2^rn" + j*\/2^V^TT- ^\/2^^^H~T) ; 

(a) for one Poincare or mixed Poincare-Rossby mode and two Rossby or Kelvin or zero modes 
T{n,k,j) + T{n*,k*,j*)-T{m,k + k*,i) - 7/3 (jV2^r+T + j*V2^V^TT - £\/2^^n^) ; 

(Hi) for two Poincare or mixed Poincare-Rossby modes and one Rossby or Kelvin mode 

3C= C{n,n\m,k,k*) > 0, \T{n,k,j)+T{n*,k*,j*)-T{m,k + k*,e)\ > C; 

(iv) for two Kelvin modes and one Rossby mode 

C 

3C= C{n,n*,m,k,k*) > 0, \T{n,kJ) -\- T{n* ,k* J*) - T{m,k + k* ,£)\ > — ; 

(v) for two Rossby modes and one Kelvin mode 

C 

3C = C{n, n* ,m,k,k*) > 0, \T{n, k,j) + T{n* , k* ,]*) — T{m, k + k* ,£)\ > —■ 

Proof — The proof of these results is based on Lemma 2.4.2. 

(i) In the case of three Poincare or mixed Poincare-Rossby modes, Lemma 2.4.2 provides 

T{n,k,j)+T{n*,k*,f)-T{m,k + k*,£) = V/? (.?\/2^n~T + jW2n* + 1 - £^2m + l) + o{^/]3), 
and it is easy to check, using considerations of parity, that the constant 

{jV2n + 1 + j* V2n* + 1 - £^/2m + l) 
cannot be zero. 

(ii) In the case of one Poincare or mixed Poincare-Rossby mode, we have one term which is exactly 
of order y^ whereas the others are negligible compared with \/]3, thus the sum is equivalent to the 
Poincare mode, and the same formula holds 

T{n,k,j)+T{n*,k*,j*)-T{m,k + k*,£) = ^ {jV2^r+T + jW2n* + 1 - £V2m + l)+o{^). 

(iii) The third case is a bit more difficult to deal with, since the leading order terms can cancel each 
other out. Without loss of generality, we can assume that £ = and j,j* ^ (the other cases being 
obtained by exchanging j,j* and ~£). 
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If j = j* , OT ii j + j* — and n ^ n* , the same arguments as previously show that the same formula 
holds 

T{n,k,j)+T{n*,k*,f)-T{m,k + k*,£) - 7/3 {jV2n + 1 + jW2n* + 1 - £V2m + 1) , 

since the factor of ^/]3 is not zero. 

If .1 + j* =0 and n = n*, the factor of V/? cancels and we have to determine the next term in the 
asymptotic expansion : 

rin,k,j)+Tin*,k*,f)-rim,k + k*,e) = ^ f-^^* -^±^ + o{l). 

2{2n + 1) 2m + 1 

Considerations of parity show therefore that the limit cannot be zero if fc + fc* 7^ 0, or equivalently if 
T(m,fc + fc*,0) ^0. 

(iv) In the case of one Rossby and two Kelvin modes, we are not able in general to prove that the 
leading order term, i.e. the limit as /3 ^ 00 of T{n, k,j) + T(n* , k* ,j*) — T{m, k + k* , (.) is not zero. 
But we can look directly at the second term of the expansion, i.e. the factor of /3~^ : 

_ 4fc'^n(n + l) ■i{k*fn*{n* + l) A(k + k*fm{'m + I) 
^^ ~ (2n+l)4 (2n* + 1)4 ^ (2m +1)4 

Considering one Rossby and two Kelvin modes means that k, k* and k + k* are not zero, and that 
exactly two indices among n, n* and m are zero. Thus lui y^ and 



\T{n,k,j)+T{n\k*,r)-T{m,k + k*,£)\>^ 



for P large enough. 



(v) The last situation is the most difficult to deal with, since the only thing wc will be able to prove is 
that the two first terms of the asymptotic expansion of T(n, k, j) + T{n* , k* ,j*) — T{m, k + k* ,£) with 
respect to (3 cannot cancel together. By Lemma 2.4.2, we deduce that for one Kelvin and two Rossby 
modes 



with 



and 



T{n, fc, j) + T{n*,k*,j*) - T{m, k + k* ,£) = ujo + "^ + o (- 

k k* k + k* 

Wo = 



wi 



2n + 1 2n* + 1 2m + 1 
4k^n{n + l) 4:{k*fn*{n* + l) 4:(k + k*fm{m + 1) 



(2n+l)4 (2n* + l)4 (2m +1)4 

Recall moreover that k, k* and k + k* are not zero, and that exactly one index among n, n* and m. 
is zero. Without loss of generality, we can assume that m — and n, n* 7^ (the other cases being 
obtained by exchanging n, n* and m). 

Then, ii luq = lui = 0, 

kn k*n* 

2n+ 1 2n* + 1 ^ ' 

k^n{n + l) {k*fn*{n* + l) 

(2n+l)4 ^ (2n* + 1)4 

from which we deduce that 

(n + 1) (n* + l) 



0, 



n2(2n + l) (n*)2(2n* + l) 



0. 
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Therefore, as the function 

x + 1 



a;2(2a; + l) 
decreases strictly on R+, we get n = n* and thus k ~ — fc*, which contradicts the fact that k + k* y^ 0. 

We conclude that cither cjg 7^ or t^i 7^ 0, so that 

\r{n,k,j)+r{n*,k*,r)-T{m,k + k*,e)\>^-^ 

for P large enough. Lemma 2.4.3 is proved. D 



Let us go back to the proof of Proposition 2.4.1, and first consider the case when k j^ 0, k* ^ 
and k + k* y^ 0. In view of Lemma 2.4.3, the asymptotic behaviour of Pn,n*,m.k.k*{P) as /3 — > 00 is 
completely determined by the behaviour of the factor 

<Jn,n',m,k,k'W) = T(n, fc, 0) +T{n*,k*,0) - T{m,k + fc*,0). 

Indeed, Pn,n-',m,k,k*{P) is defined as a product, eight factors of which involve triads of type (i), six 
of which involve triads of type (ii) , twelve of which involve triads of type (iii) and the last factor of 
which is crk,k* ,n,n*,m{P)- By Lemma 2.4.3 we then deduce that there exists a nonnegative constant C 
(depending on fc, k* , n, n* , m) such that 

\Pn,n* ,m,k,k*{P)\ ^t^P \0'n,n* ,m,k,k* {P)\- 

If one or two among n, n* and m are zero, properties (iv) and (v) in Lemma 2.4.3 allow to conclude 
that for /3 large enough 

and thus Pn.n*.m,fc,fc* has a finite number of roots. 

If n,n*,m are all equal to zero or n,n*,m are all non zero, we cannot conclude as no estimate on 
'^n,n'',m,k,k-'{P) at infinity is available. Therefore, cither (Jn^n* ,m.k.k*{P) is identically zero for j3 large 
enough, or Pn,n*,m,k,k*{P) has a finite number of roots. 

Thus at this stage, in order to prove Proposition 2.4.1, it remains 

(1) to consider the case when k,k*,k + fc* ^ and <Jn,n»,m,k,k»{[i) is identically zero for j3 large 
enough, with n,n*,m all zero or all non zero; 

(2) to study the case when fc or fc* or fc + k* is zero (in order to establish that only the triads 
involving a zero mode may be resonant). 

• Conclusion in the case of (1) 

In the case when n,n*,m are all zero, then the resonances corresponding to cro,o.o,fc,fc' (/?) — are 
precisely Kelvin resonances, which cannot be removed. 

In the case when n,n* ,m are all non zero, then CTn,n*,m,fc,fc* (/3) is an analytic function of (3 (the roots of 
(2.2.3) - defined explicitely with Cardan's formula - do not cross each other according to Lemma 2.2.1, 
and thus depend analytically on /?) : in particular, if (T„^„»^m,fc,fc*(/3) cancels for /3 large enough, then 
it is identically zero. Let us describe the asymptotics of the roots as (3 goes to zero. 
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Lemma 2.4.4- — With the notation of Paragraph 2.2.1, for all k (^ Z and all n € N*, the following 
expansions hold as P ^ : 

T(n,fc,0) = |+o(l). 

Proof — Since the product of the roots goes to zero as f3 goes to zero, we infer that at least one root 
goes to zero with /3. Let us consider that root. Since f3{2n + 1) is neghgible with respect to fc^ and r^ 
is negligible with respect to k^r, we find that 



fcV - fc/5 - 0, 



so that one root is equivalent to '^ as /? goes to zero. It is easy to sec that the two other roots are then 
equivalent to ±fc, so that we do have T{n, fc, 0) ^ £ (we recall that for n ^ 0, the roots are numbered 
in increasing order). The lemma is proved. D 

Now going back to the study of case (1), in view of Lemma 2.4.4 it is obvious that CTn,n',m,k,k'iP) 
cannot vanish indcntically. 

• Conclusion in the case of (2) 

In this situation, we need to refine the previous analysis by introducing an auxiliary polynomial. We 
thus define 

IkM* - {{.h.fJ) e {-1,0, l}^/e^Oiik + k* = 0, j^Oiik^O and j* 7^ if fc* = O} 
and 

Qn,n',m,k,k'{P)= JJ {rin, k, j) + T{n* , k* , j*) ~ T{m, k + k* , £')) . 

That corresponds to the remaining possible resonances, where we have omitted the trivial case when 
one wave is geostrophic (r = 0). 

As previously, considerations of symmetry show that this quantity can be rewritten in terms of the 
symmetric functions of (T(n, fc, j))jg{_i g ij (or {T(n, fc, j))jg|_i_i} if A; = 0), the symmetric functions of 
(T(n*,A:*, j*))j.g|_i,o,i} (or (T(n*, /c*, j*))j.g|_i,i} if fc* = 0) and the symmetric functions of (T(m,fc + 
k* ,£))ie{-i.i} (or {T{m,k + fc*, ^))£g{-i,o.i} if fc + fc* = 0). Noticing that the symmetric functions 
of (T(n, 0, j))jgr 1 ij, arc affinc in /3, we conclude that (5n.ri*.m,fc,fc*(/3) is a polynomial in /3. 

The asymptotic analysis of the various factors as /? ^ 00 shows that 

|Qn,n*,m,fc,fc*(/3)| > C (3 

for /3 large enough. Therefore, Qn,n*,m,k,k'{(3) has a finite number of roots, meaning that there exist a 
finite number of /? such that resonant triads with fc = or fc* = or fc + fc* = (other than the triads 
involving a non-oscillating mode) can occur. 

We have therefore proved that 

(1) in the case when fc, fc*, fc + fc* 7^ and cr„^„*^m.fe.fc*(/3) is identically zero, only the triads involving 
three Kelvin modes are resonant for an infinite number of /?; 

(2) when fc or fc* or fc + fc* is zero, only triads involving zero modes are resonant for an infinite 
number of /3. 
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Combining this result with the conclusion of the previous paragraph achieves the proof of Proposi- 
tion 2.4.1. n 

2.4.2. The special case of KerL. — In this short section we are going to write an algebraic 
computation which in particular allows to derive the following proposition. 

Proposition 2.4-5. — Let $ and <i>* be two smooth vector fields. Then for every n € N, we have 

(*nAo|OL(<i>,$*))^.(j,,T)=0- 

Remark 2.4-6. — That proposition implies that the projection of the limit system (2.3.5) onto KerL 
can he formally written 

(9tno$ - vA'Jlo^ = 0. 

Proof. — We are going to prove a more general result, computing the quantity 

where $Aj *&/j and $^ are three eigenmodes of L associated respectively with the eigenvalues iA, i/i 
and ijl where X = fi + fi. The proposition corresponds of course to the case when A = 0. 

We have 

hence denoting by $a the complex conjugate of $a, we get 

V ' ' ^ ^ / L^(RxT) 

^\j (^A,oV • ($^,0$^ + $A,o*;) + $A • ('&; • v$;i + <^'f, ■ v$;,)) dx 

= \j (-v#A,o • ($^.0$;^ + *A.o*;) + ^'x ■ (v($;, • ^'f,) + ^'^^^ ■ ^'f, + ^f^"- ■ *;,)) dx 

= -\j {{p^i^'x + v*A,o) • ($^,0*;^ + $A,o*;,) + (V • $^)$; • <i>y dx 

Using the identities 

V • $^ = i\<Px.o, 

/3a;i$l^ + V$A,o = «A$^, 
as well as their combination 

/3$A,i - iA(V^ • $'^ - /3a;i$A,o) 
and similar formulas for $^ and $^, we get 

($a|Ol($m,$a)) 

V ^ '^ /L2(RxT) 

= ^ / (*Ai.^ • ($^,0$;^ + $^,0$;.) + *A#A,o$;. • $^) ^x 

" 2^ / ^""-^^'^ + ''^^^^^ ■ "^^ ~ /32^i'^M,o)(V^ • $;! - Px,^f,,o)dx 
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from which we deduce 

2(RxT) 






/" i'A,2(V^ • $; - /32;i$^,o)(V^ • $;, - /3a;i$^,o)dx 



(2.4.9) _ ^ 

2/3 

- ^ J(y^ ■ ^'x - /32;i$A,o)($M,2(V^ • ^'f, - Px,^^,o) + $A.2(V^ . $;, - /3a;i<I>^^o))rfa;. 
In particular for A = this quantity is always zero, which proves Proposition 2.4.5. D 



CHAPTER 3 



THE ENVELOPE EQUATIONS 



The aim of this chapter is to study the system {SWq) obtained formally page 19 as the limit of the 
filtered system (2.3.3) as e ^ 0. Let us recall the system: 

[bWo) \ 

where A'j^ and Ql denote the linear and symmetric bilinear operator defined by (2.3.6) page 19. 

Two different types of wellposedness results will be proved on (S'Wq): first we will prove the existence 
of weak solutions in L^ and of a unique, strong solution if the data is smooth enough (on a short time 
interval, which becomes infinite for small data). Then we will show that except for a countable number 
of /?, the strong solutions exists globally in time as soon as the initial data is only in L^, of arbitrary 
norm. 

The statements of both theorems can be found in Paragraph 3.2, and their proofs are respectively 
the object of Paragraphs 3.4 and 3.5. In order to establish those results we will need to define, in 
Paragraph 3.1, suitable function spaces, compatible with the penalization operator L as well as the 
diffusion operator. Some technical preliminaries devoted to those spaces are proved in Paragraph 3.3: 
in particular in Paragraph 3.3.2 we prove the continuity of the bilinear operator Ql in those function 
spaces. Finally the last part of this chapter is devoted to an additional smoothing property on the 
divergence. 



3.1. Definition of suitable functional spaces 

By construction the operators A^ and Ql appearing in the limiting filtered system {SWq) are defined 
in terms of the projections (Ilx)i\^e on the eigenspaces of L. In particular, they are not expected to 
satisfy "good" commutation properties with the usual derivation V. Therefore in order to establish 
a priori estimates on the solutions to {SWq) we have to introduce some weighted Sobolev spaces 
associated with some derivation-like operator which acts separately on each eigenmode of L. 

Let us therefore introduce the following norms. We will write as previously ^n,k,j for the projection on 
the eigenmode ^n.fc.j of L and 11^ for the projection on the eigenspace associated with the eigenvalue iX 
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of L. Finally we define 

5 = Nx Z X {-1,0,1}. 

Definition 3.1.1. — Let s > be a given real number. We define the space iJ£ as the subspace 
of (L^(R X T))'^ given by the following norm: 

\in,k,j)€S 

Due to the definition of the eigenvectors of L seen in the previous chapter, one can prove the following 
proposition. 

Proposition 3.1.2. — Let s > be given. Then one has the following property: 

V<PeHl mni - Wild- A + [3^xly/H\\L.^^,T,y 

In particular, H'l is continuously embedded in LL^IH x T), and for all compact subsets fl ofRxT, 
HqI^I) is continuously embedded in iJ£. 

Moreover for all $ G Hf^ n (KerL)-*-, we have 

mHl - J2 l|nA<J'll?/MRxT) 

\iAe6\{0} 

Finally if ^ belongs to KU P, as defined in Definition 2.2.4 P^-S^ 16, then 

\iAee\{o} 
Proof. — Let us first prove the first equivalence: let $ € -ff£ be given. Then we have 

\\{ld-A + f3^xlr/H\\l, = \\ J2 (Id-A + /32x2)^/2n„,fe,,$| 

By the identity 

-^P'; + iSxli^n ^ I3{2n + l)lbn 

the orthogonality of the family {ipn)ne'N a-nd the explicit formulas (2.2.6), (2.2.7) and (2.2.9) for '^n,k,j, 
we infer that, for all integers a, 

||(Id - A + p^xlr^ii^^k,! - (1 + n + er^n^kjL^ <C{l + n + Pf-^ 

which implies in particular that 

(Id -A + /3a;2)-*„.fe.,|(Id - A + ffxlY'^n.k.jA <C{l + n + k^y-\ 

' '■' ' '■' /L2(RxT) 

On the other hand, 

((Id - A + pxlY-^n,k,,m - A + pxlY^n^^k',j^) ^ = if n ^ n* or fc ^ r , 

so we find that 

!|(Id -A + /3x?)-<i>||i.(j,,T)- EE((Id -^ + /^^?)''n„,fc,,'&|(Id -A + /3a;2)-n„,fe,,.$)^, 

n,k,j j" 
n,k,j 



L2 
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We then obtain the first equivalence for all s > by interpolation. 
Then, from the inequality 

V<i>eHl, imH'iR.T) < C\\{ld - A + pxly/'^i^^^^^T,^ 
along with the fact that and for all $ e C°°(R x T) supported in [—R, R] x T, 

||(Id -A + f3xlr/H\\mnxT)<C{l + Ry/^nH^^nxT) 
we get the embeddings H^{n) C HI C H-'iK x T) for all fi CC R x T. 

The second result of the proposition is easy, using Proposition 2.2.5 page 17: 

V*Ae6\{0}, ^ Y. \\^n.kMHHnxT)<mxnHHnxT)<Cs ^ l|n„A.$llH=(RxT), 

T{n,k,j) — X T{n.k.j)—X 

and recalling that by Proposition 2.2.3 page 12, we have 

-^(1 + n + ky/^ < ||*„,fc,,||H=(RxT) < Cs{l + n + ky/\ 

Finally the last result, concerning Kelvin and Poincarc modes is simply due to Lemma 2.2.1 and 
Proposition 2.2.3 . 

The proposition is proved. D 

Remark 3.1.3. The Hf^ estimates are both regularity and decay estimates. In particular, the 

embedding H'1 C L^ (R x T) is compact, and we have the following equality 



fl ffi = 5(R X T) 



s>0 



Note that these spaces are also used by Dutrifoy and Majda [7] to study the uniform wellposedness of 
a non viscous version of {SWe). 



3.2. Statement of the wellposedness result 

The main results of this chapter are the following two theorems. We have written n_L$ for the 
projection of $ onto (KerL)-*". In the next theorem, we state the global existence of weak solutions 
and the local in time existence (and uniqueness) of strong solutions. 

Theorem 2 (Wellposedness results for all /3). — There is a constant C such that the following 
results hold. Let $° G L'^(R x T;R^) be given. Then 

• there exists a global weak solution $ G i°°(R+;L^(R x T)) to (SWq), such that Il±^ belongs to the 
space L^(R+; Hj^), and which satisfies for every t > the energy estimate 



^ll*Wlli^(RxT) + ^^* llv(no$)'(t')||i.(RxT)rfi' + ^l^ llv(n^<i>)(i')lli.(RxT)'i^' < ^ll*"lli^(R 



xT) 



• if we further assume that IIo^'^ belongs to Hf^ for s > 0, then Hq^ (which is unique) belongs 
toL-(R+;i/i). 
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• if Hq^^ belongs to _L^(R x T) and 11^$'^ belongs to H^ , then there exists a maximal time inter- 
val [0,T*[, with T* = +00 under the smallness assumption 

l|no$°||L^(RxT) + lin^^^llffV^ < c-'iy, 

such that <^ is the unique solution to (SWo), anrfn^$ belongs to L^^{[0,T*[, Hl^^)r\Lf^^{[0,T*[, hI^^). 

• if n^$o belongs to HI for some 1/2 < s < 1, then n^$ belongs to L^J[0,T*[,HI) D 
LU[0,T*IHI+'). 

The previous theorem is much improved if a countable set of values for /3 is removed. 

Theorem 3 (Wellposedness results for generic /?). — There is a constant C and a countable 
subset M of R+ such that for any (3 e R+ \ Af, the following result holds. Let <i>° e L^(R x T; R"^) 
be given. Then (SWo) is globally wellposed, in the sense that there is a unique, global solution $ 
in L°°(R+; _L^(R x T)) such that n^$ belongs to the space L^(R+; H\), and which satisfies the energy 
inequality of Theorem 2. 



• if we further assume that 11^$'^ belongs to H^, for < s < 1, then Il±^ belongs to Lj^^(R+, i/£) n 

LUn+,Hi+'). 



Remark 3.2.1. These results are based on a precise study of the structure of (SWo), o,i^d in par- 

ticular of the ageostrophic part of that equation, meaning its projection onto (KerL)-*-. One can prove 
in particular that the ageostrophic part of (SWq) is in fact fully parabolic. That should be compared to 
the case of the incompressible limit of the compressible Navier-Stokes equations, where again the limit 
system is parabolic, contrary to the original compressible system (see [5], [8], [25]). Note however 
that (SWo) actually satisfies the same type of trilinear estimates as the three-dimensional incompress- 
ible Navier-Stokes system, which accounts for the fact that in Theorem 2 unique solutions are only 
obtained for a short life span (despite the fact that the space variable runs in the two dimensional 
domain R x T). In the case of Theorem 3, we use the study of resonances of the previous chapter 
which shows that the limit system is linear, except for its projection onto Kelvin modes; but Kelvin 
modes are essentially one- dimensional so energy estimates are much improved compared to the case of 
Theorem 2, and that is why global wellposedness is true in L^ , for arbitrarily large initial data. 

The rest of this chapter is devoted to the proof of those theorems. Some preliminary results are proved 
in Section 3.3 below, namely the fact that the ageostrophic part of the limit system is parabolic, along 
with trilinear estimates. In Section 3.4 we prove Theorem 2, whereas the proof of Theorem 3 can be 
found in Section 3.5. The last section will be devoted to an additional regularity result, giving an 
estimate of the divergence of $' in both cases, which will be useful in the next chapter. 



3.3. Preliminary results 

Let us prove some results that will be used throughout this chapter: in Section 3.3.1 below, we prove 
that the limit system, projected onto (KerL)^ is parabolic. In Section 3.3.2 we prove crucial trilinear 
estimates. 



3.3. PRELIMINARY RESULTS 33 

3.3.1. Parabolicity of the ageostrophic limit equation. — In this section wc are going to 
prove that the projection of the hmit system onto (KerL)^ is parabolic. To obtain that result, the 
important remark is that, for each eigenmode of L, the first and third components of the eigenvectors 
(corresponding to rj and U2) have very similar behaviours, and thus controlling the regularity of the 
last two components is sufficient to have an estimate on Il±^ in Hj^. A result of quasi-orthogonality 
in H^(H X T) of the nonzero eigenmodes of L leads indeed to the following result. 

Lemma 3.3.1. — Let s > be given. There is a constant Cs such that for any $ G (KciL)^, we 
have 

\mis+.<c,m-A'^^)Hi, 

meaning in particular that the projection of the system (SWo) onto (KerL)^ is fully parabolic. 



Proof — The proof of that result consists in using the structure of the eigenmodes to prove that the 
diffusion - acting a priori only on the velocity field - has also a smoothing effect on the pressure, and 
more precisely that 

(3.3.1) V(n,fc,j)e5, l|n„,fc,,<i>||H.(RxT) < C'||(n„,fc,,<i>)'||ff.(RxT), 

where Tl„^k,j denotes the projection on the eigenmode '^n,k,j of L (with the notation of the previous 
chapter) and C" is a nonnegative constant (independent of n, k and j). By formulas (2.2.6), (2.2.7) 
and (2.2.9) we deduce that for any integer s, we have 

|!92(*„,fej)o||L2(RxT) == l|t^2(*n,fc,i)2||L2(RxT) 

using the orthogonality of ipn-i and V'n+i, and that 

7^l|t^l(*n,fej)2||L2(RxT) < l|5l(^n,fcj)ollL2(RxT) < C'll^^l (^n,fc,j)2 || L2(rxT) ■ 

This implies in particular (3.3.1). 

By Remark 2.2.6, page 18 we then deduce that for all <I> € (KerL)^ 

($|-Ai^<i>)i.(RxT)- E (nA<i'|-Ai^(n;,$))^.(j,^T)^ J2 ll(nA$)'|ll,,(RxT) 

iAee\{o} iAee\{o} 

(n,fc,i)es* 

- CC' ^ llnn,fcj*ll|i(RxT) 



recalling that 



1(1 + n + ky/' < ii*„.fc,,ii^.(RxT) < c^(i + " + k'r^'- 



We therefore obtain the first inequality using Proposition 3.1.2. 
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In a similar way, by Proposition 2.2.3, page 12, we have for all $ e (KeiL)-^ 

($1 - a'^$)h£ = Y. (1 + " + fc')'(n,.fc,,$| - A^$)i2(K,,x) 

(n,fcj)eS* 

{n,k,j)£S* T(n* ,k* ,j*)—T(n.k.j) 

(«,fcj)es* 

T(n,fc,j)#±^/^ 

where i?^ is the contribution of the modes ±\//3/2, defined by 



/c=±^/3/2 

E (i + fr 

fe=±A//372 



0j*)e(o,signfc)2 
(ij*)e(o,signfe)2 



ffi 



Using (3.3.1) and Proposition 2.2.5 page 17 leads then to the expected estimate. 



D 



Remark 3.3.2. — Note that these inequalities indicate in particular that the notion of homogeneous 
or inhomogeneous spaces does not make sense for these weighted Soholev spaces. 

Moreover we recall that there is no analogue of Proposition 2.2.5 in the case of geostrophic modes, so 
that the geostrophic equation does not have that ellipticity property. 

3.3.2. Derivation of the trilinear estimates. — An important step in the proof of Theorems 2 
and 3 consists in establishing some control on the nonlinear term arising in (SWq) in terms of the 
Sobolev norms introduced in Section 3.1. Such estimates are obtained using classical para-differential 
methods. Obviously a more general statement could be written, at the price of more technicalities. In 
order to keep the proof as simple as possible we choose to state only those estimates that will be used 
in the following. 



Proposition 3.3.3. — Denote by Ql the limit nonlinear operator defined by (2.3.6), and let a be 
any real number greater than 3/2. Then the following trilinear estimates hold : 



\u^<i>i]i;m^nli\^,^^ 



+ q|n^$||^i||n^$*||^il|n^<i>,|i;^/f||n^$,||^/'(j,^T) 



($,|gi(<i>,$*)) 






\^o'^*\\LHRxT)m±nm 



L2(RxT) 



+ Ca||ni$,||i2(RxT) l|no$||L2(RxT)l|n±$*||ffi 



|no$*||L2(RxT)l|n±$| 



Hi 
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and for all s < 1 
(<i>,|Qi($,$*))^. 



<c ||n^a>| 



H-: 



ini$* 



\m 



(l|no*| 



in, $* 



H- 



in,*! 



m 



ini$ 



-L^*llH? 



^a + . vil^^U*lli^(R.xT)l|ii±'i' l|H£ + l|no$*||L2(R,xT)l|n_L$||H£ 



,3/2||ni$*llt,.||n!$ 



* II 1 — s 



|n_L$* 



Iff- 



in,*! 



h; 



ni$|| 



h; 



+c||nx$*|| 

< c ( ||nj_<i>||^3/2||ii_L-i' 11^ 

+qinx$*||H.+i (||no$||L2(RxT)l|ni$1lH£ + mo<^*\\mRxT)m±nHi) ■ 

Remark 3.3.4- — l.The estimates presented in that proposition are exactly the analogue of the usual 
trilinear estimate for the three-dimensional Navier-Stokes equations. For instance in three space di- 
mensions one has 



I** 



"Sll-S, ,J|V$*||7-2fR31 + ||$*|| . 3, , l|V$||r2|-R3 



< C||$J| . 1 

-II *IIh7(R2) 



l*llHi(R2)ll^'^*lli^(I^^) 



l$1l^i(^.)l|V$|k2(H2) 



whereas in two space dimensions one would expect 

The reason for the loss of one half derivative compared to the usual two dimensional case is linked to 
the fact that differentiation with respect to xi corresponds to a multiplication by ^/n instead of n. 

2. The restriction s <\ is due do the particular structure of the nonlinear term, in particular to the 
coupling between Rossby modes. For the Rossby mode associated to the eigenvalue iX the regularity 
is indeed measured by 1/A. Therefore the condition of resonance X = fi + fi (which is of course not 
equivalent to A"^ = /i^^ + ft^^ ) does not allow to distribute the derivatives as in the usual paradif- 
ferential calculus. Note nevertheless that the computation (2.4.9) page 27 allows actually to distribute 
one derivative and to obtain a trilinear estimate of the form 



($|Ql($,$))^. <c!|ni$||^a+i||ni$||Hj (!|ni*!|^3/2 + ||no$|U2(RxT 



for all s <2. 



Proof. — The method used to establish these estimates is rather standard : we decompose each 
vector on the eigenmodes of L, then compute each elementary trilinear term, and finally determine 
summability conditions. The fundamental result we will use to estimate the sums is the following 
(see [21]) 

VweF(Nx Zx {-1,0, 1}), Vu;e^'''°°(Nx Zx {-1,0,1}), v*wer(TSix Zx {-1,0,1}) 

(3.3.2) . Ill 

with p,q,r g]1, +oof and — = — I 1, 

r p q 

where the convolution is to be understood in k and n, coupled with the classical result 

(3.3.3) ((l + n + fc2)-i)e^3/2.°°(NxZx {-1,0,1}). 
In the sequel we will use the following notation 

S* = N X Z X {-1, 0, 1} and S** = N x Z x {-1, 0, 1} \ N x {0} x {0}. 

We have by definition of the space _ff£ , 

($*|gL(*,$*))H£ = Y. (i + «* + fc*)'(n„..fc..j.$*|n„.,fc.,,.Qi($,$*))^^. 

(n.,fc.,J.)eS 
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We can then write 

(<i>*|Ql(<i>,$*))h£- E 






E {l + n, + kly[nn,,k,,3,<^*\Un,.k,.j.Q{Ilx<5,Ux,<^*) 



L2 



Let us start by estimating the purely ageostrophic part of Qj^, denoted Qj^ and defined by 



(3.3.4) 
We have 



Qi($,$* 



J2 n;,Q(n^$,n^.$* 






Hf 



E (1 



fc?)Mn„.,fe.j.$, 



n 



n* ,A;* j' 



,Q(n„,fcj$,n„.,fc-j***) 






L2 



where the eigenvalues iT(n, k,j), iT{n^, fc*,.j*) and iT{n*, k* ,j*) run over 6 \ {0}. 

Thus using the regularity estimates on the eigenvectors (^n,fcj) of L stated in Proposition 2.2.3, 
page 12, we get (writing to simplify r for T{n,k,j), and similarly r, ~ T{n^,k^, j^) and r* — 
T{n*,k\f)) 



E {i + n,. + kly (n„^,fc.j.** n„.,fe.j.g(n„,fc,i$,n„.,fc.j.$*) 



fc*=fc+fc* 



L2 



< C E (! + "* + fc*)1(^*)„.,..,,JI^«A.II'P:%.*.,. I((" + fc')'^' + in* + (r)2)V2), 

fc*=fc+fc* 

iT,iT,,ix*ee\{o} 

where ^Pn,k,j is defined as in the proof of Proposition 2.2.3 by 

For the sake of clearness, we will simplify (abusively) the notations as follows : we will denote re- 
spectively by ip, (p» and ip* the coefficients Pn,k,j, {'^*)n,,k,.j, and (/3*, j,* -., and by ^ the sum over 
(n, k,j){n^, k^,j^,){n* , k*,j*) G (S*)^ satisfying the following constraints 

T(n,, k^,jf) — T{n, k,j) + T(n*, k* ,j*) and fc* = fc + fc*. 

It is fundamental for the following estimates to notice that those constraints in fact imply that when k 
and k* as well as j,j* and j, are fixed, then the condition T{n,,,k + k*,j^) ~ T{n, k,j) + T{n* ,k*,j*) 
implies that a given n and n* constrain the value of n,. Indeed we recall that according to Lemma 2.2.1, 
page 9, there is only one value of n, associated with one value of T{ni,, k.,,jt:) ^ 0. Note however that 
contrary to the usual case when there is an actual convolution (as is the case for the Fourier variable k 
here), we have no ovbious estimate on n*, as a function of n and n* . So the usual methods of 
distribution of derivatives cannot be fully used here, as derivatives in the xi direction, acting on $, 
cannot be traded for derivatives on $ or <I>*. 
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(i) If s = , by the Cauchy-Schwarz inequality, we obtain 



1/2 



1/2 



By (3.3.2) and (3.3.3), we therefore get 

^|^,||^||^*|(n + fc2)i/2 

x[\\{n + k'y/'^\U.^S')m + n. + k^^)-'/%2.^^s)y.^^^^^^^ 
and a similar estimate for the symmetric term. 

By definition, 

(3.3.5) |l(l + n + fc2)i/2^||,,(5.)<qin^$|l^,. 

Plugging this estimate in the previous inequality leads to 



1/2 
1/2 



(3.3.6) 



(i>*IOl(*, 



$* 



L2(RxT) 



1 3/4 1 



1/4 



3/4, 



1 3/4 



<c\\n^^4'^m±^4'i!^\^.T)m^n%m±<^*\\% 



X ( ||n^<f>1ll/f ||n^$|li{{R,T) + lin^^lll/f lin^**"'/' 



Iff 



Iff 



l2(rxt; 



1/3 



1/3 



1/3 



(ii) If 3 — . another way to estimate the L^ scalar product is as follows 

^|^,|M|^*|(n + fc2)i/2 

< (e ((1 + - + k')'/'\v\y^' ((1 + -* + k'.)'^'\v.\f' (1 + -* + (krr'^'^ 

x(^^((l + n + fc2)3/4|^|)'/'((l + „* + (r)2)V4|^*|)'/'(l + „^+^2)-3/4^ 

X (E(i + " + ^')"'^' ((1 + "* + kiy^'\^*\y^' ((1 + ^* + (fc*)')'/'i^i)'^') 

and 

^|^,||^||^*|(n* + (r)2)V2 

< (5:(n* + (r )2)V2|^*| ((1 + n + fc2)V2|^|) '/' ((1 + n. + k^'^^^f' (1 + n. + fc^)-^^^ 

X ('^(n* + (r)2)l/2|^*| ((1 + n + fc2)l/2|^|)'/' (^(1 + „^ + fc2)l/4|^^|^^/'^ (^ ^ ^ ^ fc2)-l/4\ ' 

By (3.3.2) and (3.3.3), we therefore get 

^ |^,||^||^*|(n + fc2)l/2 < 11(1 + „ + P)3/4^||^,(^. J(l + „^ + fc2)l/4^^||^^^^^^ 

x\\{l + n* + {k*)Y^^V*\\i%S')m + n + k')-'/%..^^S') 



1/2 
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and 






So we find 
(3.3.7) 

(*.|g.($,$*))^^ 



(RxT) 



<c||n^$||^3/.||n^$,!|^/f||n^$,||V2^^^^l|n^<j,*l|V2||n^<j,*l|V^^^ 



1/2, 



1/2 



c||n^$||^,||n^co*||^,||niCO,j|V;=||n^$,l|v^^^^^^. 



< 



1/2 



(iii) If s = and we have additional regularity on $ and $*, tlien again one can write a different 
estimate. We can write indeed 

(3.3.8) Y. I^*ll^ll^1("* + ik*f)"^ < ||^*||,^(s*)ll(l + n* + (r)2)"/Vl|,2(s.) 

X 11(1 + n* + (r )2)l/2-«/2||^^^^. J(l + „ + fc2)^/V||,2(s.)||(l +n + k')-"^^,.^s'), 

which, coupled with the similar estimate for the symmetric term, gives the expected result. 

(iv) If s < 1 , we have by Holder's inequality 

^(l + n,+fc2)^|^,||^||^*|(n + fc2)V2 

^(1 + n, + k^r\^,\ ((1 + n + k^f/%\y^' ((1 + n* + (F )2)i/2|^*|) '^^ (i + „ + fc2)-i/2\ 
X f^(l + n, + fc2)-|<^^|((l+„ + fc2)3/4|^|y (^(l+n* + (r)2)i/2|<^*|^ (l + n* + (r)2)-M 
from which we deduce 

^(1 + n* + k^,r\^,M^*\{n + fc2)l/2 < Cj|(l + „^ + fc2)s^^||^^^^^^||(^ ^ ^ ^ fc2)3/4^||^^^^^^ 

x||(l+n* + (r)2)V2^*|j^,(^.j||(l+„+fc2)-3/4|j^^^^^^^^ 

and a similar estimate for the symmetric term. 
Therefore, we get 



1/2 



**IOl($,$* 



Hf 



<c\\ii^<i>4Hrm±nH-/4^^'^*\\Hi + m±'^*\\H-/4^^nHi): 



and we conclude by interpolation 
(3.3.9) 



$*|Ql($,** 



Hf 



< c ||ni$ 



±>i'||fr3/2 



H' 



|ni$*||?,J|ni$* 



ff-' 



ini$ 



±^ Wh^^ 



ini$||?,J|ni$ 



Hf\\'-'-±^\\H- + ^ 



X ||n^$,||^-/||n^$,||^ 



• Proposition 2.4.5 shows that 

Ql($, $*) = OL(no$, n^**) + Ql{ti±<^, Ho**) + Ol($, $*). 

In order to end the proof of the proposition, it remains therefore to estimate the terms coupling the 
geostrophic and ageostrophic parts. We start by noticing that the constraint T{n,k,j) — T(n',k,j') 
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implies by Proposition 2.2.1 page 9 that necessarily n — n' and j — j' , except if n = 0. But according to 
Remark 2.2.2 that case corresponds to two different values of j for the same eigenvalue only if 2fc^ = (3, 
in which case the multiple root is fc — ±-\//3/2. That means that one can write 



($*|QL(no$,ni<i>*))^e |< 






xTl 



E 
E 



(1 + ^y (no,,,/i7i,o'J'*l^(no$,n„ /— ^,.$ 



2' 



(1 + '-^y ino,v^./*IQ(no*,n„, /^^o** 



L2(RxT) 
L2(RxT) 



Integrating by parts when the derivative acts on IIo'I', we get 



(n„,fc,,$,|Q(no$,n„,fc,,.$*)) 



L2(RxT) 



<C|((^,)„,fc,,||(^*)„,fc,,.|(" + fc')'^'l|no$||L2(RxT) 



By the Cauchy-Schwarz inequality, we then get 

\{<P.\QL{no<i>,n^^*))Hi\ 

< q|no<I'||L2(RxT)||(l + n, + kly^\4,2^s')\\{l + n* + (fc*)2)(-+i)/V*|l,2(s.) 
Remark that the derivatives can be distributed either on $, or on $*. 
We finally deduce that 

'$*|gL(no$,$*) 



(3.3.10) 






Note that this term is not zero as in the case of the usual Sobolev spaces, because the spectrum of L 
is not symmetric with respect to 0. 

One should also remark that in (3.3.10), no derivative acts on any vector field in KerL. This can seem 
somewhat surprising, but is due to the very strong constraint induced by the resonance: instead of a 
summation over three types of indexes (namely (n, fc), (n*, fc*) and (n^, fc*)), one only sums over n. 



Combining (3.3.6), (3.3.7) and (3.3.8) with (3.3.10) (with s = 0) provides the first estimate of the 
proposition, while (3.3.9) and (3.3.10) (with s < 1) give the second one. 



The proposition is proved. 
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3.4. Proof of Theorem 2 



The proof of Theorem 2 is divided into four steps. In Paragraph 3.4.1 is proved the existence of weak 
solutions, and the propagation of regularity of the geostrophic part is proved in Paragraph 3.4.2. The 
construction of strong solutions is performed in Paragraph 3.4.3, while the propagation of regularity 
of the ageostrophic part is proved in Paragraph 3.4.4. 
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3.4.1. Weak solutions. — In this section we are going to prove tlic existence of weak solutions to 
the hniit fihered system (SWo). We foUow the hnes of the classical proof of the Leray theorem, stating 
the existence of weak solutions to the Navier-Stokes equations. 

• Definition of the approximation scheme 

Denote by Kn the truncation operator defined by 

(3.4.1) Kn= J2 n„,fcj. 

Clearly the operator K]yQ]^{K]y^, K^^) — vKjq/^'j^Kjq^ is continuous on L^(R x T) (with a norm 
depending on N). Therefore, we deduce from the Cauchy-Lipschitz theorem that there exists a unique 
maximal solution $(jv) e C{[Q,Tn[,L?'{^ x T)) to 



(3.4.2) ' 

$(Ar)(t = 0) = i^W$°. 

Note that the uniqueness implies in particular that '^im) ^ Kjq'^ijsiy 



Now let us write an energy estimate on (3.4.2). The quadratic form being skew-symmetric in L^ (this 
can easily be seen by its definition as the limit of the filtered quadratic form in (2.3.3) page 18) we 
find that 

2^ll*(JV)(*)lli^(RxT) - '^(AL$(iV)|$(Ar))L2(RxT) = 0. 

Applying Lemma 3.3.1 implies that 

-(A'j^nx*(Ar)|ni$(iV))L2(RxT) > C'o"^ l|n±$(Ar) |1 hI(rxT) 

so we infer by Gronwall's lemma that 

ll*(JV)(^)lli2(RxT) + ^^ j \\^i^a'^iN)Yit')\\h^n>cT)dt' + 2^J' ||n^<I'(Ar)(i')llffi(RxT)'^^' 
(3.4.3) < ||i^iV$°|li.(RxT) < ll<I'°lli^(RxT), 

SO that the approximate solution is defined globally, i.e., Tj^ — +00. Moreover the proof of Lemma 3.3.1 
also implies that 

n^$(w) is bounded in L^(R+; H^). 
• Existence of a weak solution 

We will only sketch the proof of the existence of a weak solution, as it is very similar to the case of 
the 3D incompressible Navier-Stokes equations. By (3.4.3) we deduce that 

(($(Ar))o) is uniformly bounded in L'^(R+,L^(R x T)) 
($(^)) is uniformly bounded in i°°(R+, L2(R x T)) n Lj*„^(R+, //^(R x T)) and 
n^$(Ar) is uniformly bounded in L°°{K+,L^{K x T)) n L^{Il+,Hl). 
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For any h > 0, denote by (5^(f>(^)(i,a;) — $(Ar)(i + h,x) — $(7v)(i,a;)- Then 

\\^h^{N){t)\\L^([a,T]x-RxT) 

T / pt+h \ 

U/i$(7V)(t) / 9t$(7v)(s)ds| dt 

\ * /L2(RxT) 

ft 

('5h$(Ar)(i)|gL(*(Ar)(s),$(Ar)(s)))^, dsdi 

"'t ^ ^ 

By Proposition 3.3.3 and the positivity of — A^, we deduce that 



\^h^{N){t)\\L^([a,T]xRxT) 







3/4||TT,X.A,.„^+Mll/4 llTT,^,,.,^„MlV4 



<c / / ||n^<5,$(^)(t)||^ri|n^<5.*w(i)lli4(KxT)lin±$wWli;;: 



X \\n±<P^N)is)\\%\-^^T,.^dsdt 



j-T pt+h 

+ C / ||(5,,n_L*(Ar)(t)||L2||no$(Ar)(s)|U2(R,xx)l|n_L$(Ar)(s)||^irfsdi 

Therefore, using the uniform L°°(R+,L^(R x T)) bounds on ^(n) and dh^(N)j and the uniform 
L^(R+,iJ]^) bounds on n^$(7v) and U±6h^(N) coming from the energy estimate, we get by Holder's 
inequaUty 

and thus 

||<5?i*(Ar)(0llL2([O,T]xRxT) ^ Crh^/^ . 

By interpolation, on gets therefore that (up to extraction) 

no$(Ar) ^ no$ weakly in L2(R+, Lj^„^(R x T)) 

<i>'(^) ^ $' strongly in LL(R+,iL(RxT)) 

n±$(Ar) ^ n_L$ strongly in LL(R.+ , ^^(R x T)). 

Note that, because of Remark 3.1.3, the last convergence is actually global in space. We are then 
able, as in the usual case of the 3D Navier-Stokes equations, to take limits in the weak formulation 
of (3.4.2), which proves that $ is a weak solution to {SWq). 

• Strong-weak uniqueness 

In general such a weak solution is not unique and the Cauchy problem is not well-posed in _L^(R x T). 
Nevertheless we have the following strong-weak uniqueness principle. 

Proposition 3.4- 1- There is a positive constant C and a nondecreasing, positive Junction C(t) 

such that the following holds. Let $ and $* be two weak solutions to [SWo) with respective initial data 
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^^ and ^1, satisfying the energy estimate. Assume that there exists some T > such that 11^$ belongs 
to L°°{[0,T], hI^'^) n L^i[0,T], hI^^). Then for all te [0,T], the function S<P = <P^ - <P satisfies 

wmmh < \\mo)\\h c^v{cmi + iia>o|ii.) + ^(i + m^nl^^^^^.y^Hi/^^) I m^^mi^^.^dt' 

In particular, $* = $ on [0, T] x R x T if $° == (f>°. 

Proof. — In order to establish the stabihty inequaUty we start by writing (formally) the equation on 
(5$ = $* - $ 

(3.4.4) dt6<i> + gL((5$, 6^) + 2Ql{6<P, $) - j/A^(5$ = 0. 

Proposition 3.3.3 implies that 

|(,5$|Ql(j$,$))^. |<q|ni$ll v2||ni<5$||Hi||ni<5$|U2 
+ c||ni,5$|U2(||no$|U2||ni,5$||ffi +||no(5$|U2||ni$||^i). 

We then deduce (using the same argument as in the construction of a weak solution page 40 for the A'^^ 
term), that 

ll^*Wlli2(RxT) - ll'5<i>°lli2(R,T) + 2i^/ ||no(<5<i>)'(t')||L2(RxT)dt' + ^y' \\n^5Hf)\\l,^dt' 

< a J^ (||n^<i>(t')ll^3/. + ||n^$(t')||4^,) \\u^5^t')\\l,dt' 
+ c,J (||no$(t')|li2||ni<5$(t')||i2 + ||nx$(t')||?,i||no(5$|||,)dt', 

using the embedding Hj^ C L^. Gronwall's lemma yields 

mml. < mYL^e^p^c^J (i + \\iioHt)\\1,^^,^^ + ||n^<f>(r)r^, + ||n^a>(r)||^3/.)rfT), 

and the conclusion comes from the fact that n^$ belongs to L*{[0,T], Hj^) by interpolation be- 
tween L°°{[0,T],Hj^ ) and L^{[0,T],Hj^ ), along with the energy estimate on $. The proposition is 
proved. D 

3.4.2. Propagation of the geostrophic regularity. — The following regularity result for the 
geostrophic equation is inspired by the Weyl-Hormander symbolic calculus, even if that theory does 
not seem to be appliable directly due to the possible singularity at xi — 0. 

Using the formula giving XIq in Proposition 2.1.2, we first see that the geostrophic equation 

dt^ - j/HoA'no* = 
can be brought back (at least formally) to the scalar equation 

dtU2 - vD{DD'^ + Id)-^D^Au2 ^ 0, 



where we recall that D is the differential operator defined hy D- = di [ — — . Then by a simple 

\l3xi ' 

change of variables this scalar equation becomes 

dt^p — vAip — 0, 
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where A is sonic self-adjoint scalar pseudo-differential operator (possibly singular at xi 
principal symbol of which is given by 



0), the 



a(xi,$i 



e^ 



neglecting the possible singularity at xi = 0. Then, in order to propagate Sobolev regularity on tp, we 
should have to control some commutator of the type [A, V^], which is not so easy because A cannot 
be written simply in terms of the usual derivatives V. 

In order to find a convenient way to measure the regularity, we therefore use formally the results of 
symbolic calculus. Note that the Weyl-Hormander theory is used here just to guide intuition, the 
result of propagation being actually proved by explicit computations. In order to determine the class 
of operators A should belong to, we have first to characterize the metric. Computing the partial 
derivatives of a with respect to xi and ^i 



dxia 



(a;i,6 



2/32x1 



%a 4 

a 4i 



26 



shows that the Hormander metric to be considered is the one associated to the harmonic oscillator 

„2 



g{dxi,d£,i) = /3 



2 '^•^1 



dei 



1 



x^ 



1 1+e^ 

Then it is natural to measure the regularity by powers of the harmonic oscillator, and therefore to 
study the propagation equation 

dt{-dl,, + (3'xly^ i^A{-dl,^ + p^xfr^ = i^[(-dl,^ + p'xfy^A]^ 

The fundamental result of the Weyl-Hormander theory states the following: if A is a pseudo-differential 
operator (meaning in particular that there is no singularity at xi =0), the commutator occuring in the 
right-hand side of the previous equation is a pseudo-differential operator of lower order (for the metric 
g), meaning that we expect || [{-dl^^^ + P^ xfy , A]ip\\ ^2 (^^^ to be controlled by |1 (-9^^^^ + P^ xfy ip\\ l^ (r) 
and -{{-dl,^+/3^xly^\A{-dl,^+f3^xly^). 

Nevertheless, as we are not able to prove in a simple way that there is no singularity at xi = 0, we 
shall not use the general theory of pseudo-differential operators and will proceed instead using explicit 
computations. We have seen in the previous chapter that the family (^ri,o.o)neN defined by 

/ 



*n,o,o(a;i) 



/ (n+1) 
2{2n+l' 



■■)pn-l{xi) 



(n + 1) 



V 



■1pn-l{xi) 



2(2n+l 



2(2n-hl) 



-1pn+l{xi) 



\ 



V'„+i(a;i) 



/ 



2(2n+l) 
constitutes an Hermitian basis of KerL, and that, due to the properties of the Hermite functions, 

Vn e N, \\{-dl^x^ + iS^xly^n^ofiW^R) ^ (1 + ")' ll*".o,o|li2(R) • 
Therefore it is natural to study the propagation of the i/£ norm of IIo^, recalling that 



|no*lll,j = E(i + ")1^„l' 



riGN 



where we have defined 






(*„,o,o|no$). 
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In the following we will denote by Ng the operator defined by 

Vn e N, 7V,no$ = ^(1 + n)V^*„,o,o, 
so that 

l|A^«no$||L2(R) = ||no$||H£- 

We have 

• We start by computing noA'^„^o,o- From (2.2.1) page 9 we deduce that 



/ 



dx^'i'nfifiixi) 







and 



^a;ia;i*n,0,o(^l^ 



V Y 4(2n+l) 

/ 

2^2(2n + l) 



/3(n + l)(n-l) , , , //3n(n + l), , , //3n(n + 2), , . 



2n+l 



4(2n+l) 



\/{n + l)(n - l)(n - 2)'(/;„_3(xi) - (3n - l)\/riTTV'n-i(xi) 



y +(3n + 4)V^V'n+i(a;i) - ^n(n + 2)(n + 3)'(/'„+3(xi) y 

with the usual convention that f/'n ^ for n < 0. Therefore, using the orthogonality of the Hermite 
functions in _L^(R), we get 

with 



a. 



(-4)^_P / (n-4)(n-2)(n-l)(n+l) 



(3.4.5) 



4 y (2n+l)(2n-7) 

a(-^)^^(4n-2),S"-^;i"+^\, 
4^ ^Y (2n+l)(2n-3) 

i-n-, /3 6n^ + 6n — 1 

4 (2n + 1) 

«(2)^^(4„ + 6) 



n(ri + 3) 



a. 



(4) 



(2n+l)(2n + 5) 

13 /n(n + 2)(n + 3)(n + 5) 
'ly (2n+l)(2n + 9) 



• From the previous computation we deduce that 

[AT,, noA']*„,o,o = ((n - 3)« - (n + l)«)al-^)*„_4,o,o + ((n - 1)^ - (n + l)^)at'^*„_2,o.o 
+ ((n + 2)^ - (n + l)^)ai'^*„+2,o,o + ((n + 5)* - {n + l)*)al^)*„+4,o,o 
Thus, using the definition (3.4.5) of the coefficients a, we get 

\a„\<C{n + l) 
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and 

||[iV„noA']*„,o,o|li2(R) < ^(71 + 1)" < C,||A^,*„,o.o||l2(r). 
Because of the quasi-orthogonality of ([7Vs,noA']^„^o,o)neN, we have actuahy the more general com- 



mutator estimate 



|[iv„noA']no$|li.(R) 

= 5] |((n + 1)^ -in- 3)^)ai-lV„+4 + ((" + 1)' " (" " ^y)^t'^Z 



+ {{n + 1)^ -{n + 3r)<U^_, + {{n + 1)^ - [n + S)^)^^!^^^ 



.n+2 

2 



< c, Y,{n + 1)2- (w,.^,r + \^„^,v + i^„_,r + 1£ 



.n+J + \^n+2\' + \^n-2\" + \^nJ' 



n 

which can be rewritten 

(3.4.6) |l[iv„noA']no<J>|l^.(R) < cWNMonmu)- 

Note that, due to the particular choice of the operator Ng, there is some additional cancellation, 
meaning that the commutator [A^sjIIoA'] which is expected to be a pseudodifferential operator of 
order (2s + 1) is actually of order 2s. 

• It is now very easy to propagate regularity using Gronwall's lemma. We recall that 

9tiv,no$ - z.noA'iv,no$ = i.[a^„ noA']no$, 

from which we deduce that 

\\Ns^^Ht)\\hiTL) + ^ I ||v(7v,no$)'(T)||i.(j,)dr<||iv,no$"||i.(j,) + c, / \\N,n^^^{T)\\l.^^^dT 

Jo Jo 

and finally 

(3.4.7) ||7V,no$(i)||i.(R) + i^ I ||V(7V,no$)'(r)||i.(i,)dr < WN^H^-^Yl^^^) exp(at). 

Jo 

This concludes the proof. 

3.4.3. Local strong solutions. — In this section we are going to prove the existence of unique, 
strong solutions for smooth enough initial data. As in the case of weak solutions discussed in Sec- 
tion 3.4.1 above, we will not write the full proof, but detail the estimates enabling one to use the usual 
Fujita-Kato theory of strong solutions to the 3D Navier-Stokes equations (see [4] for instance). 



• Global existence of strong solutions for small data 

We prove here that under a suitable smallness assumption there exists a (unique) global strong solution 
to (5Wo) such that n_L$ belongs to i°°(R+, iJ^^^) n L'^{B + , H^^"^). 

As previously we start from the solutions '^im) of the approximation scheme (3.4.2). We have of course 

^||no<I'(Ar)(i)||i.(j,,T) + 2^l|V(no<i>(Ar))'(i)||i.(j,,T) < « 
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and, by Proposition 3.3.3 
dt 



^M^^iN){t)\\%u. - 1v (ni<I>(Ar)|A^nx$(iV))^i/2 (t) 



< -2(ni$(jv)|niQL(*(Ar),$(Ar)))^i/2 (t) 

3/2 1/2 

By Lemma 3.3.1 and the obvious embedding Hj^ C i^^ , that inequahty can be written 

4lin^$(^)(t)ll^,,, + 2-^||niCi>(^)(t)||^3,, 

(3.4.8) 1/2 



< C||n_L$(W) 11^3/2 (^||n_L$(Ar) 11^1/2 + ||no$(7V)||L2(RxT)j(t)- 

As usual we notice that this inequahty is useful only if ||n^$(^-)(t)|| 1/2 and ||no$(jv)||L2(R^xT) s-re 
small, which is a typical phenomena of global results under a smallness condition. 



/ \ 2~ 

/ V 



Define 

DN = \i^ R+ /Vt' < i, ||no<i>(Ar)(t')lli2(RxT) + l|n±<I'(Ar)(t')ll^i/2 < (^^^ ^ 

where C is the constant appearing in (3.4.8), and let us impose the following smallness assumption on 
the initial data: 

l|no$^)lli2(RxT) + lin^<J'^)ll^v2 < (^^ 

Then clearly Djy is not empty. By construction, no$(Ar) belongs to C(R+,L2(R x T)) and nx(f>(^-) 
belongs to C(R+, Rj^ ) thus D^ is a closed set. 

Denote by Tat = max I? at. If T^v < +00, then 



|no$(„)(rAr)||i2(R,T) + l|n±<f>w(rAr)||^^,/, < (^^^ 



1/2- 

and we deduce from (3.4.8) that 
d 

which is in contradiction with the definition of Tn = maxD^r. Therefore Tj^ = +00. 



Then we deduce immediately that for all t E R+ 

I|n^<i'w(t)ll^v2 + j^ / ||n^<f>(Ar)(t')ll^3/2rft' < ||n^$?w)||^v2. 

i '-^1/2 Jo ^ ^ 

Up to the extraction of a subsequence that converges to a Leray solution $ of {SWq), the previous 
estimate implies that 

The strong-weak stability principle established in the previous section provides then the uniqueness of 
such a solution. 

• Local existence of strong solutions 
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Let us now consider the case of large data. The idea (see for instance [4]) is to spUt <I>(7v) in two parts 
as follows 



*(JV) = '^w + '^In) 



where ^1}^) is the unique solution to 



^t*(Af) ~ "^L-^iN) 



(3.4.9) 



•I'wW 



Y^ n„,fcjni$('jv) + no$(Ar), 

(«,fcj)es 



and A > is a truncation parameter to be determined (independent of N). Using Proposition 3.3.1, 
it is easy to check that 

||ni<i>^)||ioo(R+^i2(RxX)) < ^11$ ||l2(RxT), 

By (3.4.2) and (3.4.9) we deduce the equation satisfied by ^J^n : 

-JArgL($fjv),*fAr)) + 2JAr( 
/ 



9t^tN) - '"^'l^In) + '/jvQL($fjv), ^fAT)) + 2JArgi($>^), $<^)) = - JArQi($<^) , $<^)) 



(3.4.10) 



<J'm(0) = $?Ar)- 



(N) 



E 



n 



ni$?^,, + no$ 






n,k,]l-l-±'i'{N) 



1/2 



O'i'(Af) 



/ 



We are going to show that $7^-, remains small in Hj^ on a time interval which does not depend on N 

-I /o 

Let us write an energy inequality in iJ^ on (3.4.10): we have 

^^II*WWII'^V2 + ^ll*fAr)(t)ll^3/. = -(JivgL(<I>f^),$fA.))|$fA.))^i/2(t) 

1/2 

+ 2(JArgL($>^),$<^))|$>^))^V2(i) - (JArgL($fjv)'<J'w)l*w)<^W' 

and Proposition 3.3.3 yields 

+ C||$f^)(i)||^i/2||no*(Ar)(t)|U2(RxT)||ni$^Ar)(t)llH^/2. 

Now consider the set 



D 



N 



teR+/yt'<t, ||$>v)(*')llffV2 < 



(where C is the constant appearing in the right-hand side of the previous inequality) and T^ — sup-Djv- 
We are going to prove that there exists T > such that 

VA^eN*, Tn>T. 

We notice that if A is chosen large enough (independently of N), then Dn is not empty: we can indeed 
choose A so that 

(3.4.11) ||cO>^^(0)||^,/, <^^,^, ViVeN. 



4CCi 



/2 



48 CHAPTER 3. THE ENVELOPE EQUATIONS 



As long as t < Tjv, we can write the previous inequality in the following way: 



/2 

+ 7ll*fAr)WllLv^(l|n±$^Ar)WII^3/.|| + ||no$(Ar)(i)||i.(j,,T: 

+ c||$fjv)(t)||^i/2||no$(w)(t)l|L2(RxT)l|nx$fjv)WllH^/2 

1 ^ - .4, 



So we get 



where 

F{A, $", ^) = C (^1 + ^ j (1 + A-)(l + \m\hii..T) 
Gronwall's lemma enables us to infer that for all t < Tn, 

Since A is chosen so that (3.4.11) is satisfied, it suffices now to choose T in such a way that 

exp {TF{A, $0, i.)) < 2, exp {TF{A, $«, i^)) - 1 < ^ [c^^ 
so that for any t <T, 

hence necessarily 

V7V £ N*, Tat > T. 

• Gathering those results, we infer that any limiting point of ^fj^\ + ^nv-, (in particular any Leray 
solution to {SWq)) satisfies 

no$ e L°°(R+, L^(R X T)), (no$)' e l2(r+, i/i(R x t)) 
nx$eL-([o,T],i/y')nL2([o,T],i/^/'). 

The weak-strong stability principle gives then the uniqueness of such a solution on [0,T]. 

3.4.4. Propagation of the ageostrophic regularity. — In order to construct a strong approxi- 
mation of the filtered solution to the Saint- Venant system in the next chapter, we will actually need 
further regularity on the solution of the limit filtered system, which is obtained in a very standard 
way from the trilinear estimates stated in Proposition 3.3.3. So suppose that XIj^^" belongs to iJ£, 
with 1/2 < s < 1, and consider as previously the sequence ($(Ar)) of approximate solutions to (SWq) 
defined by (3.4.2). From Lemma 3.3.1 and the last estimate in Proposition 3.3.3, we deduce that for 
all s < 1, 

'^"" ^ '-12 , 2i/„„ ^ ,,,,,2 



-||ni<j>(Ar)(t)|l^. + — ||n^$(jv)(i)||^.+i 

< C||nx$(Ar)(t)||//= ||ni$(Ar)(i)||^. + i (||nx$(Ar) (t) ||^3/2 + \\nQ<^iN)m LHUxT] 

< ^l|n±$(jv)(t)||^.+i + ^\\U^-i^N)m%i (l|n±*(Ar)(t)llHj/2 + ||no<I'(Ar)(t)|U2(a^T)^' 
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and thus by Gronwall's lemma 

^ L-s Jo ^ 

Taking limits as A^ ^ oo in the previous inequality shows that 

n^$eL,-([o,T*[,i/£)nLL([o,T*[,H^+i), 

and that 

l|ni*(i)llH£ + 7^ f m^Ht')\\l.^,dt' 

^ Cs Jo "^ 

<||ni$"|||.cxp('^y" (||ni$(t')||^3/2 + ||no$(f')llL^(RxT))'dt'') , 

which proves the propagation of regularity result, and completes the proof of Theorem 2. D 

3.5. Proof of Theorem 3 

In this section we shall prove Theorem 3: Paragraph 3.5.1 is devoted to the global wellposedness result 
in L^, while the propagation of regularity results are given in Paragraphs 3.5.2 and 3.6. 

3.5.1. Global wellposedness. — In this section we shall prove the first part of Theorem 3, namely 
the fact that except for a countable number of f3, the limit system is globally wellposed in L^. This 
turns out to be an easy matter in view of the resonance results obtained in Section 2.4 in the previous 
chapter. 

Indeed Proposition 2.4.1, page 19 indicates that except for a countable set of values for (3, the limit 
system reduces to the following: 

dtHo^ - z/noA'^$ = 0, 

dtUR<^ + 20i(no$, Uii<P) - i^UrA'l<p = 0, 

dtUM^ + 2Ql{I1o<P, Hm*) - /^HmA'^* = 0, 

dtiip<P + 20L(no$, np$) - z/npA'^$ = o. 

So the limit system is a linear equation on all modes but Kelvin modes; Kelvin modes being essentially 
one-dimensional, it will be easy to prove the wellposedness of the system. In fact the only point to be 
proved is the uniqueness of the solution, since existence was proved in the previous section. Uniqueness 
will be immediate in the case of all non-Kelvin modes so let us concentrate on the equation on 11^^ <I>. 
Let $ and <I>* be two solutions, and define (5$ = $* — $. Then 

dtiiK6^ + 20i(no(5$, 0^,5$) + 20i(no$, 0^,5$) + 20i(no<5$, n^*) + Ql{tIkS^, n^,5$) 

+2QL(nx$,nx5$) - lyllKA'^S^P = 0. 

Now let us write an energy estimate in L^, in the spirit of the computations of Section 3.4.1 above. We 
note that in the case of Hk, the decomposition on eigenmodes of L simply corresponds to the Fourier 
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decomposition, so that 

(nK<5<j'|QL(nK$,nK<5<i'))i.(R^T) - (nx<5<J>|0(nK$,nK<5$))^.(i,^T)- 

Moreover in that case, the usual Sobolev spaces H^ coincide with the 7J£ spaces, since n = 0. Using 
the fact that Ho projects onto a;2-independent functions and that the dependence in xi of HkS^ is 
that of the Gaussian i/jq; it is easy to see that 



and that 

(nK<5a>|0(nK$,nK<5a>))^.(i,,T) 



< ||nx5$||Hi/RxT)l|no$||L2(R)||nK(5$||i^/R.i2/T)), 



(R)l 



= (R;L2(T)), 



< ||nK<5$||ffi(RxT)||nK$||L-(RxT)||nK(5*||L2/R,xT 



(RxT)l 



2(RxT) 



C. 



- ^lin^'^*llHi(RxT) + -linK$|lio.(RxT)linK'5$||i2(RxT)- 



We recall indeed that 



IIk • = V! (*o.fe,o I • )L2(RxT)*o.fc,o where *o,fc,o(a;i, 2:2) = ^==e*''' 



fcez* 



^ M^i) ^ 



(^ V'o(a;i) J 



One sees easily that |jnif<i>||^oo(RxT) < C'||n/i-<i>||^i(RxT)i so finally an energy estimate (coupled with 
a Gronwall lemma) gives 

Xexpf^f \\^K^t')\\j,^ dt'Y 

Since IIkS^{0) — 0, uniqueness follows from the energy bound on 11^^$. Note that we have recovered 
here the usual, two-dimensional Navier-Stokes type estimates since in the case of purely Kelvin modes, 
the quadratic form and the spaces involved are the same as in the Navier-Stokes case. In fact Kelvin 
modes are even one-dimensional (up to a multiplication by ^'0(2:1)) so it is possible to improve those 
estimates (that will be done in the last section of this chapter). 



3.5.2. Propagation of regularity. — Let us now prove that if the initial data IIj^^" is in i/£ 
with s e [0,1], then that regularity is propagated to n^<i>. Recalling the special form of the limit 
system for almost all /3, the only equation we need to study is the one on Uk^- Indeed denoting 

* = (nfl + np + nM)$, 

we have 

9t* + 2QL(no$, *) - lyA'^^ = 0, 
and Proposition 3.3.3 gives directly 
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Now let us turn to the Kelvin modes. In that case we simply use again the fact that H'' and _ff£ spaces 
coincide in the case of Hk- We have, using Proposition 3.3.3 again, 

l|nK<f>(t)||l,.(RxT) + ^/ l|nK$(i')llH= + i(RxT) < / (IiKmi^K'^,^K^))HHnxT)(t')dt 

+ c f ||no$(i')||LH|nK<J'(t')llHnRxT)l|nK<J'(t')llH-+MRxT)di'. 

Jo 
Two-dimensional product rules give 



< ||nK$||if=+i(RxT)l|nK$ ® UkHh'Crxt) 

nK$|| 1, 

^ "ff2(RxT) 



so that finally 



< l|nK<i>|||.+i(RxT)linK*llH=(RxT)l 

< ^mKnWiuxT) + 7!inx<i>||l,.(RxT)linx$||S^^^^^ 



l|n/.$(t)||?,.(RxT) + ^£l|nK$(i')lll^=+i(RxT)<^^ \\TioHt')\\l4iiKHt')\\l.^^,^^ 



dt' 



^^VK<i>(t')ll?.=(RxT)linK<i'(t')ll^(^^^)< 



and the result follows from Gronwall's lemma and the energy estimate on ni^:<I> 
That concludes the proof of Theorem 3. 



n 



3.6. A regularity result for the divergence 

In this section we are going to prove an additional regularity result for the system [SWq), which will 
be useful to study the strong asymptotics of the rotating shallow-water system in the next chapter. 

Proposition 3.6.1. - Let $" belong to ^^(R x T). 

For all 13 € R+, if Il±<^° belongs to H^^^ and (Up + 0^)$" belongs to Hf; for a > 3/2, then the 
solution $ of (SWo) with initial data $*^ defined on [0, T * [ satisfies for all t G [0, T*[ 

f ||V-$'(t')||L~(RxT)di'<+«3, 

Jo 
where we recall that $' denotes the two last components of ^. 

Furthermore the regularity assumption on the initial data can be relaxed for all but a countable number 
of (3. Indeed, for all j3 G R+ \ M where J\f is the countable subset of R+ defined in Theorem 3, if 
(lip -|- nx)<I>° belongs to iJ£ for a > 1/2, then the solution $ of {SWq) with initial data $° satisfies 
for all t G R+ 



/ ||V-$'(t')||L~(RxT)di'<+«^- 

Jo 



Proof. — Such a result is established by decoupling the equations on the various parts of $, and proving 
that the regularity is propagated for the Poincare and Kelvin modes, while a smoothing property on 
the divergence holds for the nonoscillating, Rossby and mixed Rossby-Poincare modes. 
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Let US decompose $ on the supplementary subsets Kcri, R, M, P and K 

$ == Eo* + I\-r(^ + IIm* + np$ + 11;^$, 
and estimate each projection separately. 

• For nonoscillating, Rossby and mixed modes the smoothness of the divergence is not due to a 
propagation result but to a stationnary property of the eigenvectors. 

Clearly, by definition of KcrL, V • (Ho*)' = 0. 

Let us consider the Rossby and mixed Rossby-Poincarc modes. By definition of the Rossby modes we 
deduce the following relation 

V$ e i?, V • $' = ^ V • $';, = ^ iA($A)o 
with the notation $a = H^^. It follows that 

w^-'^'wi^ - II E '^("^ML 

< c J2 IIM$a)oII?^2(r,t)> 

by Proposition 3.1.2. But, as Rossby waves correspond to j — 0, we have 

||a($a)o!I1/2(r.xt) < c'lAp 2^ ll(n„,fc,o$)o||?/2(R.xT)' 

using Remark 2.2.6. Recalling the explicit form of ('^n,k,j)o, we see that 

But for Rossby modes, the following asymptotics hold as |fc| or n goes to infinity: 

^^"("'^'°)^F+/3^2n+l)- 
So we infer that as |A:| or n goes to infinity, 

|A| ||(n„,fc,o^)o|lH2(R,xT) < C'|l(^„,fc,o$)ollHl(RxT)• 
Finally we infer that 

llV-^'ll^,. < C J2 I1M$a)o|1L2(hxT) 
iAeSi? 

- ^ Z2 ll(nn,fe,0'S)o||Hi(RxT) 

{n,k,a)eeR 

By the embedding of H^(R x T) into i°°(R x T) we conclude that V • (Bt?*)' belongs to the 
space L^([0, r];i°°(T x R)). The same result can easily be extended to the mixed Poincare-Rossby 
modes (it is in fact easier since n = in that case) and we obtain 

l|nM<I'||L2([0,T],_H"i) < C't, ||V • (Bm$)'||l2([o,T],L~(RxT)) < Ct ■ 

Finally we deduce that 

II V • ((IIo + IIr + Bm)$)'||l2([o,t],l°°(rxt)) < Ct ■ 
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• In order to establish a similar estimate for the Poincarc and Kelvin modes, we prove that the equation 
governing these modes propagates the HI regularity without restriction on s. Indeed we have seen in 
Remark 3.3.4 that the restriction to s < 1 for the propagation of regularity stated in Theorem 2 is due 
to the coupling between Rossby modes. 

The idea here is to study the propagation of the following norm on (KerL)^ 

which controls the i/£ regularity of the Poincare and Kelvin modes only, due to the following easy 
estimate (see Proposition 3.1.2): 



(3.6.1) 



C-^mK + Hp^Hf < lln^^lls < C'iWillK + IlpmiHf + ||n^$||L2(RxT 



(RxT)j- 



This norm is convenient to deal with the condition of resonance occuring in the nonlinear term of 

(SWo). 

Similar arguments as in Proposition 3.3.3 allow to write the following trilinear estimate: 

(3.6.2) |(<i>|Qi($, <&)),! < c,||ni$||^(||nia>||,+i + \\n^^H^J{\\u^^\^:^r^ + ||no$|U2(R^T))- 

With the same notation as in the proof of Proposition 3.3.3, we have indeed 



{^,\Ql{<^,'^*))= 



def 



Yl (1 + ^*)' (n„.,fe.j.** n„^,fc.j.Q(n„,fej$,n„.,fe.j.$ 



T^—T + T* 

■i-r*,ix*ee\{0} 



L2 



< Cs ^(1 + r!y^'{{l + T^yl^ + (1 + (r*)2)^/2) |^.||^||^*|((n + k^^l^ + [n* + {k*fY"\ 
from which we deduce that 

|($|gL($, $)). < a||nx$||,(||ni$||,+i + nn^^Hai )||ni$||^3/2. 

In the same way, for the geostrophic part, we have 

($|gL(no$,ni$)). < c,||ni$||,(||ni$||,+i + \\ii^^\HimMLmi^Ty 

Note that the H^ norm appearing above is used to control the gradient of the Rossby and mixed 
modes. 

Using (3.6.2) and Gronwall's lemma, we get the following propagation result 

m^HMl + ^ £ (||n^<i>(t')ll?+i + mi^Ht')\\l^A' 

<||nx$"||2cxp(^^^*(||n^ci>(i')||^3/2 + ||no$(OllL^(RxT)) dt'] . 

(see Paragraph 3.4.4 for the detailed proof), and therefore by (3.6.1) we get 

UUk + i^Pmm^ + ^ ^* ll(nx + np)$(i')|p^...rfi' 

<(||(nK + np)a>0||2,.+||$0||2,(^^^^)cxp('^ / (||n^$(t')||^3/. + ||no$(t')||L^(RxT)) dt'] . 
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By Sobolev embeddings we then deduce the expected control on the divergence 

II V • {{Uk + np)$)'|U2([o^T],L~(RxT)) < Ct 

for all T <T* where T* is the lifespan of the strong solution $. 
That concludes the proof of the proposition in the case of general (3. 

• By Proposition 2.4.1 page 19, we recall that, for all (3 G R+ \ A/", the only possible resonances are 
Kelvin resonances. 

Let us now consider the equation governing the Poincarc modes which can be seen as a linear parabolic 
equation whose coefficients depend on Ho^. Therefore it is very easy to propagate regularity once one 
has noticed that the multiplication by |A| for the Poincarc mode n^'I' is "equivalent" to a derivation. 

Introduce as previously the notation 

{n,k,j)eS 

so that 

{n,k,j)eSp 

where 

S-p = N* X Z X {-1, 1} U {0} X Z+ X {1} U {0} X Z; X {-1} U x x {-1, 1}. 

We can use Proposition 2.2.3 page 12 to deduce that for each (n, fc,j) in Sp the equation govern- 
ing (fn.k.j can be decoupled (recall that Hq^ only depends on xi): 

dt^n,k,j — '^^n,k,j{^n,k,j\^''^n,k,j)L^{IlxT) = —'^fn.k.ji'^ n,k,j\Q{'^ n,k,j ,^0^)) L^ (TLxT) 

which can be rewritten 

dt (</?„,fcjexp(-i/t(*„^fej|A'*„^fej)i2(RxT))) 

= -2<y5„,fcj(*„,fcj|<9(*n,fej,no$))£2(R.xT)exp(-I^t(*„^fej|A'*„_fcj)i2(R,xT)) • 

By Gronwall's lemma and the estimates 

|(*n,fc,j|(3(*„,fcj,no$))L2(RxT)| <Ci(n + fc2)l/2, 
-(*„,fcj|A'*„,fcj)L2(RxT) > C2{n + k^), 

we then deduce that there exists a nonnegative constant d, (depending only on ly) such that, 
(3.6.3) y{n,k,j)eSp, \ipn.k.jit)\ <\^n,k,j [0)1 cM-C An + k^)t). 

We have 

liv • (np$)'(i)|i^»(j,,T) < E \'Pr.,kAt)\ liv • (*„,fc,,)'(i)IL»(RxT) 

{n,k,j)eSp 

<C Y. \fnMAt)\{n + ky^^ 

(n,k,j)eSp 

since {^n.k.j) is uniformly bounded in L°°(R x T). Thus, by (3.6.3), 

||V • (np$)'(t)||^^(j,^T) <C E l^".fc.j(0)l cxp(-a(" + k^mn + fc2)i/2. 

{n,k,j)eSp 



(n + fc2)-i-" 
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Integrating with respect to time leads then to 

{n,k,j)eSp 

<c'J J2 \vn,kAorin+k'r] E (^ 

\(n.k.j)eSp J \in,k,j)eSp 

from which we deduce that for a > 1/2, 

\\^-i^P'^y\\LHlO,T],L^iILxT)) < C\\Iip<^"\\H2 

where C depends only on v and a. 

It remains then to establish the propagation of regularity for the Kelvin part of the equation, which is 
nonlinear and has no smoothing effect for the divergence as the Rossby part. 

The crucial point here is to recall as above that this equation is actually one-dimensional (modulo a 
smooth function with respect to xi). The propagation of regularity result proved in Paragraph 3.5.2 
implies that as soon as the initial data is in 7J£ with < a < 1, then the solution lies in L^(R+; H^^^)- 
In particular V • (n^^)' lies in L2(R+;7?£). So if a > 1/2, using the fact that if"(T) is embedded 
in L°°{T), the result follows directly. 

Proposition 3.6.1 is proved. D 



CHAPTER 4 



CONVERGENCE RESULTS 



The aim of this chapter is to study the asymptotics of the rotating shallow-water system (1.4.3) 
presented page 5. In particular wc will see that the system (2.3.5) obtained formally in Section 2.3, 
page 19, is indeed the limit system, after application of the filtering operator exp(— i_L/e). In order to 
simplify the presentation let us recall here the two main systems we will be considering in this chapter, 
namely the shallow-water system 



{SW, 



dtTj -\ — V- ( (1 + erj)u\ 
dt ((1 + £ri)u\ + V • f (1 + £ri)u (^ u\ + —(1 + eri)u-^ + -(1 + er])Vri - vAu = 



'7|t=o = ff\ "|t=o = w", 



and the limit system 

where A'^ and Ql denote the linear and symmetric bilinear operator defined by (2.3.6) page 19. 
We also recall the formal equivalent form of {SWe), 

dt{r], u) + -L{i], u) + Q ((77, u), (77, uj) - i^A'{r], u) ^ R 

e 

{v,u)\t=o = (»?",w°), 
where Q and A' are defined by (2.3.4) page 18 and 

i? = (0, -t/-^^Au). 

The study of the asymptotics of (S'Wg) will be achieved through three different methods, which provide 
three different types of results. In Section 4.1 we describe the weak limit of the weak solutions to {SWe) 
as e goes to zero, which is proved to satisfy the geostrophic equation studied in the previous chapter, 
i.e., the projection of (5Wo) onto KerL. The statement is given in Theorem 4 below. Then for smooth 
enough initial data, we prove in Section 4.2 the strong convergence of the filtered sequence of solutions 
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towards the unique solution of [SWq). The precise statement depends on the setting, as in Chapter 3: 
for all /? > we are only able to prove results locally in time (globally for small data) whereas if a 
countable set of values of /? is removed, then the convergence is strong for all times, and the smoothness 
assumptions on the initial data are less restrictive (see Theorems 5 and 6). Finally in Section 4.3 we 
propose an intermediate study between those two asymptotic results, by considering the asymptotic 
behaviour of the filtered sequence e~*^'^(?7e, u^), where (77^, Ug) is a weak solution to (SWe). We prove 
a strong convergence result towards a weak solution to (2.3.5), where unfortunately due to the lack of 
compactness of rjg in space, a defect measure remains (see Theorem 7). In order to circumvent that 
difficulty we propose an alternate system to the Saint- Venant equations {SW^), where capillarity effects 
are included. Technically the effect of capillarity is to have a uniform control on erj^ in strong enough 
norms so as to obtain an evolution equation for Ug. A strong convergence result for e~*^''^(?7e, Ug) is 
established in that new setting, see Theorem 9. 

In this chapter, many results and notation of the previous chapters will be used. However precise ref- 
erences will be made each time, so that this chapter can be read independently of the others (assuming 
the results of course) . 



4.1. Weak convergence of weak solutions 

The first aim of the chapter is to describe the weak limit (77, u) of (ry^, u^) as s goes to zero. 
Theorem 4 (Weak convergence). — Let (rj'^,u'^) E L^(R x T) and {rf^^u^D he such that 
^^^^^ \jm' + {l+e4)\ul\-)dx<£\ 

(r/^,u°)^ (?7°,u") ini2(RxT). 

For aU e > 0, denote by {-q^^u^) a solution of (SWe) with initial data {rf^,u^), as constructed in 
Corollary 1.4.1 page 6. Then {rj^^u^) converges weakly in L^^^(R+ x R x T) to the solution {ri,u) e 
L°°(R+, L^(R)), with u also belonging to L^(R+, i7^(R)), of the following linear equation (given in 
weak formulation) 

(4.1.2) ui=0, /3a;iU2 + ai?7 = 0, 
and for ah ir]*,u*) £ L"^ x i?^(R) satisfying (4.1.2) 

(4.1.3) / (7777* + U2U*2){t, x)dx + iy / Vm2 • VM2(t', x) dx dt' = / (r/"??* + u^m;)(x) dx. 

Remark 4. 1.1. — • Theorem 4 shows that the system satisfied by the weak limits of rjg and u^ is lin- 
ear. There is therefore no convective term in the mean flow: system (4-.1.2, 4-1-3) actually corresponds 
to the projection of (SWq) onto KciL: as seen in Section 2.4-2, that projection can indeed be formally 
written 

dt{v,0,u2) - iyllo{0,0, Au2) ^ 0, 

(r/,u)(t)=no(r/,u)(t)Vi>0, 
(r,,M)|t=o = no(7?°,u°). 

• Note that {r]'^,u'^) do not necessarily satisfy the constraints (4-1.2), so in general {ri,u)\t=o is not 
,0 „o^ 



equal to (rf , u ) 
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• The study of the waves induced by L, in Chapter 3, revealed the presence of trapped equatorial waves, 
which however do not appear in the mean flow described by Equation (4-1-3): no constructive inter- 
ferences take place in the limiting process, in other words the fast oscillating modes decouple from the 
mean flow, without creating any additional term in the limit system (that feature was already observed 
in [10] in the case of inhomogeneous rotating fluid equations, modelling the ocean or the atmosphere 
at midlatitudes) . This will be obtained by a compensated compactness argument in Section 4-T4-- 

4.1.1. Constraints on the weak limit. — We recall (see Chapter 1) that the uniform energy 
bound on (jjg, u^) implies the existence of a weak limit {r],u). In this paragraph we are going to prove 
that the weak limit belongs to KerL. 

Proposition 4-1-2. — Let (77°, u°) e L^(R x T). Denote by (?7e,Ue)e>o a family of solutions 
of (SWe), and by {ri,u) any of its limit points. Then, {r],u) G L°°(R+, L^(R)) belongs to KeiL, and 
in particular satisfies the constraints 

(4.1.4) ui = 0, I3xiu2 + diT] ^ 0. 

Proof — Let X; V' '= I'(R+ x R x T) be any test functions. Multiplying the conservation of mass 
in (SWe) by ex and integrating with respect to all variables leads to 

{erj^dtx + (1 + £%)«£ • Vx) dxdt = 0. 

Because of the bounds coming from the energy estimate (1.4.5), we can take limits in the previous 

identity as e goes to to get 

r r 

u ■ Vxdxdt = 0. 



Similarly, multiplying the conservation of momentum by £-0 and integrating with respect to all variables 
leads to 

(e{l + eri^)ugdt-)p + e{l + eri^)us-{ug-V)-)p + (3xi{l+erig)u^-'il;^ + {l + -rii;)rigVip + iyUg-A'iJj]dxdt = 0. 

Once again the bounds coming from the energy estimate (1.4.5) enable us to take the limit as e goes 
to 0, and find that 

(ryV • -0 + l3xiu ■ ^^)dxdt = 0. 

It follows that {r]{t),u{t)) belongs to KerL for almost all t G R^, and we conclude by Proposition 2.1.1 
page 7 that (ri,u) does not depend on X2 and satisfies the constraints (4.1.4). D 

To go further in the description of the weak limit {r],u), we have to isolate the fast oscillations generated 
by the singular perturbation L, which produce "big" terms in [SWe), but converge weakly to 0. 

Therefore, a natural idea consists in introducing the following decomposition 

(»7e,Ue) = no(%,U£) +n^(77£,Ue), 

where Ho is the L^ orthogonal projection onto Keri and IIj^ the L^ orthogonal projection onto (KcrL)^. 

The idea to get the mean motion is then to apply Ho to [SWe) '■ since L is a skew-symmetric operator, 
we have IIoL = and we expect 9tno(%,Ue) to be uniformly bounded in some distribution space. 
The difficulty comes from the fact that one has no uniform spatial regularity on rjg. That is why we 
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will actually consider the weak form of the evolution equations; the point is then to take limits in the 
nonlinear terms. 

4.1.2. Rough description of the oscillations. — The analysis of the nonlinear terms lies essen- 
tially on the structure of the oscillations. A rough description of those fast oscillations will be enough 
to prove that they do not produce any constructive interference, and therefore do not appear in the 
equation governing the mean (geostrophic) motion. The much more precise description given in Chap- 
ter 2 will not be used in this section, but will be necessary to discuss the strong asymptotic behaviour 
of the solutions in the next sections. 

In the following statement we have considered a regularization kernel defined as follows: let k be a 
function of C^(R^, R+) such that k{x) = if |x| > 1 and / ndx = 1. Then for any 5 > we define Kg 

by 

Ks{x) = S^ k{5^ x). 

Proposition 4-1-3. — Let (ry",^^) G i^(Rx T) and (77°, u") satisfy assumptions (4-1.1), and denote 
^y (('7ej "e))e>o a family of solutions of (SW^) with respective initial data (77°, u°). 

Then t]^ — ks * Vs o,nd m^ — Kg -k ((1 + eri^)u^) — u^ -\- e{ri^u^y satisfy, for all T > 0, the uniform 
convergences for a/l J7 CC R x T 

||?7e — ?7e||L°°(R+.H=(o)) ^0 as S ^ uniformly in e > 0, for all s < 0, 

(4.1.5) \\u^ — M£||L2Qo,T];ffs(o)) ^0 as 6 ^ uniformly in e > 0, for all s < 1, 

||?7£Ue - {rieU^Y \\l2 (^[QT];H'' (n)) -^ as 6 ^ uniformly in e > 0, for all s < 0, 

as well as the approximate wave equations 

edtvt + V- mf = 0, 

edtml+(3xi{ml)^ + \/r]f^esl + Sa^^, 

denoting by s^ and a^ some quantities satisfying, for all T > 0, 

SUpSUp||o-f||i2([o,T];Hi(R-xT)) < +00, 

(4.1.7) '>"^>° , 

V(5 > 0, sup ||4IIli([0,T];H1(RxT)) < oo- 

£>0 

In particular the approximate vorticity w* = V^ • m^ satisfies 

(4.1.8) edtiiu', - Pxi4) + f3ml, = egf + Sp^ 
with, for all T > 0, 

(4.1.9) sup \\pt\\L^([o.T]:L^(RxT)) < +00 and \/5 > 0, sup |l(7e|lLi([o.T]:L2(RxT)) < +00. 
e>0 e>0 

i5>0 

Proof — We proceed in two steps, first stating the wave equations for (775, mg), then introducing the 
regularization (77^,771^). 

• The first step consists in establishing some bounds for 

edt{i]e,nie) + L(T]^,m^). 

We have 

edtrie -f V- TTie = 0, 
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and 

edtirie + fixiin^ + '^rj^ — ere, 
with 

Tg = —V • {rrie (Xi We) — Ve'^Ve + vAu^. 

Let us now find a bound for r^. It is made of three contributions. The easiest to handle is Au^. 
Indeed u^ is bounded in L'^(R+,H^(R x T)), so Au^ is bounded in L'^(R+,H-^(R x T)). 

Next let us consider the nonlinear terms 

—V • (nie ® Ue) — rjeVrje — —V • (rrie ® u^) ~ -Vry^. 

By the energy bound (1.4.5) we infer that they are bounded in _L°°(R+, T4^^^'^(R x T)). 

Therefore in particular 

(4-1-10) lke|lL2([o,T];H-V2(RxT)) <Ct- 

• Now let us proceed to the regularization. We recall that the energy inequality (1.4.5) provides the 
following uniform bounds 



ll'7e||L~(R+,L2(RxT)) < C, 



"ellL2(R+^Hi(RxT)) - ^■ 



In particular we have 



1^.."^ 



%llL2(RxT) < llvl+^"£lli2(RxT) +C'^ll%llL2(RxT)l|w£||L2(RxT)ll^i£|liJi(RxT)' 

thus by the Cauchy-Schwarz inequality 
1, 



and finally 

(4.1.11) ||«£|1l2([o,t];H1(R-xT)) < C't- 
We also deduce from the usual product laws that 

(4.1.12) heUe\\L^([o,T];H-(RxT)) < Cs for all S < 0. 

The convergences (4.1.5) are then obtained by the Rellich Kondrachov theorem. 

• By convolution we get (with obvious notation) 

edtvt + V- uf = 0, 

and 

edtml + /3a;i(mf)^ + Vryf = erf + /3a;i(mf)^ - (/3a;imf)^. 
Notice that 

ximf (x) - {xime)\x) = / Ks{y)me{x - y){xi - {xi - j/i)) dy 



yiKs{y)me{x-y)dy. 
That implies that 
(4.1.13) Ximl{x)-{xim,)\x)^5K^p{y)^me, 
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where Kg (y) = d^'^K^^'{6^^x) and k'^-'^'(x) — xik{x). We therefore infer that 

/3a;i(mf)^(x) - /3(xim^)^(x) = fof (x) + e/3,54'' * (%Me), 
where for all T > 0, 

Supsup ||cr£||L2([0.T];Hi(RxT)) < +00- 
<5>0 £>0 

It remains then to control 

4 - rf + /3(5k(i) * {fi,u,) . 

By (4.1.10) and (4.1.12) we get 

V(5 > 0, sup ||sf ||li([o,t];H1(RxT)) < oo. 

£>0 

• Taking the vorticity in the second equation of (4.1.6) leads then to 

dtcvf + V^ • (/3a;i(mf)^ + Vryf) = eV^ • sf + 5\/^ ■ 4, 
from which we deduce that 

atwf + f3xiV ■ mf + /3mf 1 = eV-^ • sf + JV-^ • erf. 
Combining this last equation with the first one in (4.1.6) gives finally 

dtiujf ~ I3xi4) + I3ml^ = eV^ • sf + 5V^ ■ af , 
from which we deduce the estimate (4.1.9) on the remainder. The proposition is proved. D 

4.1.3. Proof of Theorem 4. — We consider an initial data [rp^vP] E L^(R x T) and a fam- 
ily (?7°,u") such that 

i|(K"p + (l + eaKnd.<£0, 

(r/^,u°)^ (?7°,u") ini2(RxT). 

We consider a family ((77^, Ue))e>o of solutions of (SWe) with respective initial data (77°, u") (given by 
Corollary 1.4.1 page 6), and {ri,u) any of its limit points. Finally we consider (77*, u*) in (L^ x H^){'R) 
such that (?7*, u*) belongs to the kernel of L. In particular by Proposition 2.1.1 we know that 

wj = and PxiU2 + d\ff — 0. 

Our aim is to prove that 

(4.1.14) I {-q-q* +U2ul){t,x)dx + v I [ Wu2 ■ Wu*{t' ,x) dx dt' ^ / (//"ry* + m^u*)(x) dx. 



• In order to establish such an identity, the idea is to take limits in the weak form of System {SW^), 
which will require some further regularity on (77*, u*), and then to extend the limiting equality to all 
vector fields (77*, u*) e L^ x i7^(R) n Keri by a density argument. 

Note that the classical regularization method cannot be applied here, since the kernel of L is not stable 
by convolution. In view of the explicit formula (2.1.2) page 8 giving the projector Ho (which is written 
in terms of the singular pseudo-differential operator 9i((/3a;i)~^-)), it is actually natural to consider 



4.1. WEAK CONVERGENCE OF WEAK SOLUTIONS 



63 



the Hcrniite functions introduced in the previous chapter, and we recall (see (2.2.7) page 11) that any 
element of Keri is a linear combination of the following 

/ 



\J]n 7 ^n 



1 



^27r(2n+ 1) 



—^1pn-l{xi) - J -Ipn+lixi) 







— ^■0„_l(xi) - J-1pn+i{xi) 



for n > 1, 



and {rio,uo) 



ipo{xi) * 




We will therefore restrict our attention to these particular vector fields which are smooth and integrable 
against any polynomial in a:;i (recall that 

[3x1 



ipnixi) = exp ( — 1 P„(xi V/3) 

where P„ is the n-th Hermite polynomial), and then conclude by a density argument. 
• Using the conservations of mass and momentum (SWe) it is easy to see that 

(^Tyn + m-e,2Un,2) (t, x) dx + 1/ / VWe,2 ' Vu„,2(i', x) dx dt' 



iVeVn + 'm'^.2'^n,2){x) dx + / (m^ • (ue • Vu„)) (t', x) dxdt'. 



Now we need to take limits as e goes to zero in all four e-dcpendcnt integrals appearing in that 
expression. 

Clearly the three first terms converge to their expected limits, as 

{VeVn + U^.2Un,2) {t, x) dx —> / {iJTJn + U2U«.2) {t, x) dx 

and 

l^ / VU£,2 • Vun,2(i', x) dx dt' ^t V I I VU2 • VUn,2(i', 2;) dx dt' 

for alH > 0, as e goes to zero. 
So the only term we need to worry about is the coupling term 

/ me,2Ue,ldiUn.2{t',x) dxdt'. 

We will prove in the following lemma that it actually converges to : 



lim 



fns,2Ue,idiUn,2{t' tx) dxdt' = 0, 



which is due to the special structure of the oscillations pointed out in Proposition 4.1.3. This result 
clearly ends the proof of Theorem 4, and is proved in the next paragraph. 
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4.1.4. The compensated compactness argument. — Let us prove the following lemma. 
Lemma 4-1-4- '^" With the previous notation, we have locally uniformly in t 



lim 

£^0 



TT^e,2U£^ldiUn^2(t' , x) dxdtt' = 0. 



Proof. — Let us introduce the same regularization as in Proposition 4.1.3, defining 



rj^ — rj^-k KS, u^ — Ugr-k Kg and m^ ~ m^ • ks- 



Then 



TTle,2U£^ldiUn^2(t' , x) dxdt' = 



(4.1.15) 



m^ 2'^e idiUn.2{t' , x) dxdt' 

m^ 2(^6 1 ~ '^e i)3iUn,2(i'j 2;) dxdt' 

m^ 2(^5,1 " ^e i)diUn,2{t\ x) dxdt' 

("^£,2 ^ "^e 2)'^e,ldlUn^2{i' , x) dxdt' . 



• By the energy estimates and the bounds on the Hermite functions given in Proposition 2.2.3 page 12, 
we can prove that the two last integrals converge towards zero as 5 goes to zero uniformly in e. Indeed 
for all q; > there exists some bounded subset Via x T of R x T such that (recalling that n is fixed) 

Then, for < s < 1 and for any s' > 0, 

("^£,2 — fne 2)^e.ldlUn^2it' , x) dxdt' 

< ||to£,2 -'™£,2llL2([o,T];H--- = '(0„xT))ll"e,llU2([0,T]:Hi(RxT))ll9lU„,2|lTVi.==(R) 
+ 2a||me,2|lL2([o^T]:H- =— '(RxT))ll"e,llU2([0,T];Hi(RxT)) 

which goes to zero as a then 6 go to zero, uniformly in e by (1.4.5), (4.1.11) and (4.1.5). 
Similarly, we get, for < s < 1, 

m^ 2 (""£,! ~ ^e i)diUn.2{t' , x) dxdt' 

< II"^£,2|Il2([0,T];H-»(RxT))I|w£,1 ^ W£,l||L2([o,T];ff«(a„xT))||9lU„, 211^^1. =°(RxT) 
+ 2Q:||me,2||L2([o^T];H-»(RxT))l|We,l||L2([o,T]:Hi(RxT)) 

which goes to zero as a then 5 go to zero, uniformly in e by (1.4.5), (4.1.11) and (4.1.5). 

Next we prove that for all (5 > 0, the second integral in the right-hand side of (4.1.15) goes to zero 
as £ goes to 0. We have seen in (4.1.12) that rjeU^ and consequently mg are uniformly bounded in 
the space L^([0,r];iy(R x T)) for s < 0. Therefore, for fixed i5 > 0, {rjeUeY and rn^ are uniformly 
bounded in L2([0,T] x R x T). Then, 

/ '™£,2(^£,i ~ fngi)diUn.2{i' , x) dxdt' 

< £|l'7^g^2lU2([0,T]xRxT)ll('7£'"e,l) || L2([0,T] xRxT) I15l"n,2|l L'=°(RxT) 
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which goes to zero as e ^- for all fixed (5 > 0. 

• So finally we need to consider the first term in the right-hand side of (4.1.15). We are going to prove 
that the limit of that term is zero using Proposition 4.1.4. Integrating by parts, we have, recalling 
that wf = V-^mf, 
t 



m^ 2"^e idiUn.2{t' , x) dxdt' 



t 

' ((aimf_2)mf 1 +771^2(^1 "^ll)) Un,l{^ ,X) dxdt' 
' 

t 

S \™<5 I ™c5 /Y7 ™(5 Q ^S 





t 



{{-ljI + Pxi4)mli + 4{~l3ximl^ + 82^!) + m^^i^ ■ mf )) M„,2(i', x) dxdt' 



^\j J d2 ((™f,i)2 - {ml,f - (ryf )2) u^^2{t', x)dxdt' 
and the last term is zero because 92W„,2 = 0- 
Proposition 4.1.4 now implies that 

edtiJ, - pxi4) + pm^^ - eqi + dp^, 

£dtml2 ~ Pximl^^ + d2r]l = ^^1.2 + ^^ta-^ 

edtVe + V • mf = 0, 

where qf and sf are bounded respectively in ^^([O, T]; L2(R x T)) and L'^{[0,T]:H'^{'R x T)) for 
any T > uniformly in e (by a constant depending on S) , and where pf and a^ are uniformly bounded 
in e and 6, respectively in the spaces L'^{[0,T];L'^{'R x T)) and L'^{[0,T];H^{'R x T)) for any T > 0. 
It follows that 

ct r 

"^e 2'^e ldlUn^2{t' , x) dxdt' 

I I (^^*(^^1^' - ^')' + |(/3^1^e - ^t)4 + |(/32;i?7s - 4)pi) Un,2{t':^) dxdt' 

f i^edtivtmi^) + evisl, + S^ai^) u^Mt' , x)dxdt' 
/o J 

Now we notice that 

/ / i/3xi4 - Uji)pt:Un,2{t' , x) dxdt' <C (^T^/^ ||?7* || ^^ (R,+ .i2(Rxx)) + ll'^fl|L2([0,T];L2(RxT))j 

X 11(1 + a;?)"^'^Mn,2||L==(RxT)||Pel|L2([o,T];L2(RxT)), 

and similarly 

vt'^e,2'Un,2it',x) dxdt' < CT^^^ || T^f || ^oo (R+.i2(R.xT)) ||Mn,2 1| L~ (RxT) Ikg || L2([0,T1;L2(r,xT)) ■ 

So writing 

ll^ellL2([o,T];L2(RxT)) < II V ' U^ ||l2([o,T];L2(RxT)) + ^11 V ' ('7eU^)|| L2([0,T];L2(RxT)) j 



t 
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we infer that 



lim lim — 
lim ( 6 



{13x^4 ~ Ujl)plur,^2{t' ,X) dxdt' 



= 0, and 



r]^a^2'^nMt',x) dxdt' 



= 0, uniformly in e. 



On the other hand, 







(/3a;i77* - uj'^)q'^u„^2{t' , x) dxdt' 



< C" (||?7£lU==(R+;L2(RxT)) + ll^ellL~(R+;L2(RxT))) 



X 11(1 + a;i)^'^U„^2||L~(RxT)||9el|Li([0,T]:L2(RxT)) 



^C" ||?7£||l=c(R+;L2(RxT) 



Ll=°(R+;L2(RxT)) +£|1V • (?7£M£)|1l=o(R+;L2(rxT)) 



X 11(1 + a;i)M„_2||L~(RxT)lkellLi([0,T];L2(RxT)), 



and 



so 



Ve^e 2^ri,2{t\ x) dxdt' 



^ C'||?7g||ioo(R+.i2(Rxx))||Mn,2||L~(RxT)ll'S£||Li([0,T]:L2(RxT) 



lim 



oV/3 



{PxiTjl - w*)g>„,2(i',a;) dxdi' 



= 0, for all (5 > 0, 



lim I £ 



rj^s^2'^n,2{t' tX) dxdt' 



==0, for all (5 >0. 



So we simply need to let e go to zero, then (5, and the result follows. 



D 



4.2. Strong convergence of filtered weak solutions towards a strong solution 



In this paragraph we will prove the following strong convergence theorems. We recall that Ii±_ denotes 
the projection onto (KerL)^, and that the spaces -ff£ were defined and studied in Chapter 3 and 
defined again in this chapter, page 79. 

The first result we will state concerns the case of smooth enough initial data, and requires no restriction 
on fi. 

Theorem 5 (strong convergence for all /3). — Let ^^ ~ {rf' ,vP) belong to L^(R x T), and con- 
sider a family ((ry^ , Ug))e>o such that 

\j {\4? + {l + evl)\ul?)dx<£'^ and 

\ f {\V° -??"!' + (1 + £Ve)K - "°n rfx ^ as e ^ 0. 

For all e > denote by {rj^, u^) a solution of (SWe) with initial data (77", u^). Finally suppose thatll±^'^ 
belongs to H^ and that np$° and Hk^^ belong to iJ£ for some a > 3/2. Then the sequence of 
filtered solutions ($e) to {SWe) defined by 



(4.2.1) 



$. 



-C (?7e,We), 



converges strongly towards $ in L^^^([0, T*[; i^(R x T)), where $ is the unique solution on [0,T* 
of {SWq) constructed in Theorem 2, page 31 . 
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The next theorem requires less assumptions on the initial data; on the other hand one must first remove 
a countable set of values for (3. 

Theorem 6 (strong convergence for generic (3). — There is a countable subset M of R+ such 
that for any (3 G R+ \Af , the following result holds. Let (f>'^ G -L^(R x T; R"^) be given, and consider a 
family ((ry^, u^))e>o such that 

\l {\4? + {l + evl)\ul?)dx<£'^ and 
(4.2.2) '■' 

2 / (l^° - r/"P + (1 + eTi)\ul - uX) dx^Q ase^Q. 

For all e > denote by {rie , Ue) a solution of (SWe) with initial data {r/'^ , u^) . Finally suppose that Hp^'^ 
and Hk^'^ belong to H^ for some a > 1/2. Then the sequence of filtered solutions (<I>e) to (SW^) 
defined by (4-2.1) 



converges strongly towards $ in Lf^^(R.'^;L'^(R. x T)), where $ is the unique solution of {SWq) 
structed in Theorem 3, page 32 . 



con- 



Remark 4-2.1. — • Note that definition (4-. 2.1) of $e does make sense since as stated in Corol- 
lary 1.4.1, one has a L^ bound on Ue- 

• The strong compactness of (<I>e) in _L^^^(R+, _L^(R x T)) cannot be obtained directly using some a 
priori estimates. Indeed we have a priori no uniform regularity on % with respect to the space variable 
X (besides we expect the limiting system to be a mixed hyperbolic-parabolic system). 

• The proof of both convergence results is based on a weak-strong stability property of (SW^). It 
is therefore crucial to be able to construct a smooth approximate solution ^app to $g, writing an 
asymptotic expansion in e whose first term is $. The regularity assumptions on the initial data stated 
in both theorems are precisely that enabling one to guarantee that the limit system has a unique, stable 
solution and propagates regularity. In particular it should be noted that in both cases, the assumptions 
on the initial data imply that V • $' belongs to L^([0, T]; L°°(R x T)) (see Proposition 3.6.1). Once 
the setting is posed so that the limit system does satisfy those properties, the proofs are very much the 
same in both cases. So in the following we will only prove Theorem 6, and leave to the reader the easy 
adaptations in the case of Theorem 5. 

Proof. — As noted in Remark 4.2.1 above, we will only prove Theorem 6 here. 

The idea is, as usual in filtering methods, to start by approximating the solution of the limit system, 
and then to use a weak-strong stability method to conclude. 

So let us consider the solution $ of (S'VFo) constructed in the previous chapter (see Theorem 3 page 32), 
which we truncate in the following way: 

(4.2.3) $Ar = JArn^^ + no^JV, 

where Jm is the spectral truncation defined by 



(4.2.4) J^= ^ n 



A 

JAGSj- 
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with 

©AT = |iT(n,fc, j) e & In < N, \k\ < n\ , 

and nj_ denotes as previously the projection onto (KerL)^. Finally Ho^at solves 

dtUo<i>N - z^HoA'no^AT = 

no$jv|t=o = Yl n„,o,o*"- 

0<ri<Af 

We recall that n„^o.o denotes the projection onto the eigenvector ^„,o,o of KerL. Then for all fixed 
iV € N we have (see Theorem 2) 

(4.2.5) Ho^AT belongs to L°°(R+; iJ£), Vcr > 0. 

Recall that such a result means that Ilo'I'iv is as smooth as needed, and decays as fast as needed 
when xi goes to infinity. 

Moreover by the stability of the limit system (which is linear) we have of course, for all T > 0, 

(4.2.6) lim ||no$Ar-no$|Uoo([o,T];L2(RxT)) =0. 

Note also that for all fixed iV e N, using the smoothness and the decay of the eigenvectors of L, we 
get for any polynomial Q E 'R.[X] 

(4.2.7) g(xi)$Ar e L°°{[0,T];C°^{K x T)) 
We have moreover, for all T > 0, 

(4.2.8) (||ni($ - *jv)||l-([o,t];L2(rxT)) + K^k + np)($ - $Ar)|lLoo([o,T];H£)) ^ as A^ ^ oo, 
and 

(4.2.9) (||ni($ - '^N)\\m[o,nHi) + U^k + np)($ - *Ar)jli2([o,T];H°+i)) -^oasN^^. 

Finally since Jn commutes with A'^ , the vector field (f>Ar satisfies the approximate limit filtered system 

dt<^N + JnQl{<^, *) - '^^'l'^n = 0, 
(4.2.10) „ 

$jv|*=o = Jn<^°- 

Conjugating this equation by the semi-group C leads then to 

"• {'' (;) *") + j^ {'' (;) *~) + •'-«' {'' (0 "■'' (0 *) - "^i'^ (?) *" = "■ 

using the definitions (2.3.6) of Q^ and A^. Let us now rewrite this last equation in a convenient way 

= (g - Ql) U (^ '^N,C (^) <!>n] - ^(A' - A'l)C (^ <Pn 

+ (Id - J^)Ql U ( J) *' -C (^ $") + Q^ fc (^ i^N - $), £ (" J") ($jv + *; 

Because of (4.2.8) and (4.2.6), the last term in the right-hand side is expected to be small when TV is 
large, uniformly in e, and similarly for the third term, using the stability of the limit system proved in 
the previous chapter. So we are left with the first two terms, which as usual cannot be dealt with so 
easily since they do not converge strongly towards zero. However they are fast oscillating terms, and 
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will be treated by introducing a small quantity £(j)iq (which will be small when e goes to zero, for each 
fixed iV), so that 

{dt + -^L) (c {^-^ e^Ar") ^ -{Q - Ql) U {^-^ $w,£ (^-^ <^n\ + ^(A' ~~ A^)£ (^^ $^. 

Let us now define 

(4.2.11) ^N = - Y. (x^ _~s nAg(n^^Ar,n^$Ar) + ;/ ^ — — ^n^A'n^^Ar, 

and consider 

Let us prove the following result. 

Proposition 4-2.2. — For all but a countable number of (3, the following result holds. Consider a 
vector field $° = {t]o,uo) G L^(R x T), with {Up + 0/^)$° in _ff£ for some a > 1/2. Denote by 
$ the associate solution of (SWq). Then there exists a family (r/g^Ar, Ug ^r) = C (-) ^e.N, bounded in 
the space L^^(R+,L^) n L^q^(R+, iJ^), such that {Up + nx)(?7e,v, Ue,v) is uniformly bounded in the 
space L^^(R+, -ff£) H Lf^^(Il'^ , H2^ ), and satisfying the following properties: 



• 



$e.7V behaves asymptotically as <^ as e ^ and N ^ oo 



(4.2.12) VT > 0, ^hnjwa 11$,,^. - <I'IL.([o,t];L^(RxT)) = 0; 

• for all N G N, (rje^N, Ue,v) *s smooth: for all T > and all Q G R[-'f], 

(4.2.13) Q{xi){ris,N,u^,N) is bounded mi°°([0, T]; C°°(R x T)), uniformly ine; 

• ("UcNtUc^n) satisfies the uniform regularity estimate 

(4.2.14) VT > 0, sup lim ||V • u^,n\\l^([o,t];L'^(RxT)) < Ct; 

• (??e,Afj Wg.jv) satisfies approximatively the viscous Saint-Venant system (SWe) : 

(4.2.15) dt{r]e,N, Us.n) + -L{r]s^N, Me.At) + Q {{rje.N.U^.N), {T]e,N,Ue,N)) - iyA'{rie.N,Ue.N) = -Re, AT 

e ' ' ' ' ' 

where Re.N goes to as e ^ then N ^ oo: 

(4.2.16) lim lini (||i?e,Ar||ii([o,T]:L2(RxT)) +e|l^e.A|lL~([O.T]xRxT)) =0. 

Proof. — Let us define <i>Ar as in (4.2.3) and (J)n as in (4.2.11). We can write 

<^A=- ^ ^3-^-7^ nAQ(n^$Ar,n^$Ar) 






^-=- E 771— rn^^'n.'^ 



i(A — u) 
We will check that the approximate solution ^^^n defined by 

(4.2.17) <^e.N = ^N + S(pN 

satisfies the required properties. 
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It will be useful to notice that there are two positive functions A.±(N) such that if fj, belongs to 6n, 
then either /i = or 

< A_(iV) < ^ < A+{N) < +00. 
It will also be useful to recall that, considering the asymptotics of A = T(n, k, j) as fc or n go to infinity, 

the operator > — — is continuous from _ff£ to H'[~ for any given cr G R. 
— ^ A 

Finally we recall that the spectrum of L admits only and 00 as accumulation points. 

• The correction ^^v is defined as the sum of two terms. 

Let us consider the first one, 0^ . It can in turn be written 

J^^ -n ^(1) ^n w!,(i) 

The first part, Ilo4']^ , is easy to handle since Ho is of course continuous from _ff£ to _ff£ for any a. 

Let us now study Ilj_<l)j^'. Clearly (5(n^$Ar,n^$Ar) is in iJ£ for any a > 0. We infer that Il±(l>j^' 
belongs to 7f£ for any cr > 0. Indeed \i 11 + jj, y^ then \\ — fi — fi\ is bounded from below since /x and jj 
are in (Sat and the accumulation points of A are and 00. On the other hand if /i + /i = then we use 

the continuity property of N^ — :- recalled above. 

So we find that for all polynomials Q G R[X], 

g(xi)0«ei-([O,T];C°°(RxT)), 
as well as 



This obviously implies that for all fc e N 



Q{xi)C ( -) ^^^ e L°°([0,r];C°°(R x T)). 



t \ ,(1), 



VA^ e N, lim\\eQ{xi)C ( - ) 0]v ||l-([o.t]:C'=(R.xT)) =0. 

(2) 

The second term (f))^ is dealt with similarly, splitting it into two terms: 



a(2) 



iAes,iMeejv\{o} iAee 



Because of the relations (2.2.1) satisfied by the Hermite functions, it is easy to see that the first 
contribution can be rewritten as a finite combination of some eigenvectors of L (which are smooth 
functions rapidly decaying in xi), and the second contribution is dealt with again by using the fact 

that the operator N^ — ^ is continuous from HI to Hl"^. 

A/O 

We conclude that for all Q E 'R-[X] 

Q(xi)£ (^\ ^i^) e L-([0,r];C°°(R x T)). 

and thus 

V7V e N, Inn ||eg(xi)£ (^-) 0^'||loo([o,t];CMRxT)) = 0, Vfc G N. 
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Combining these results with (4.2.5), (4.2.6), (4.2.7) and (4.2.8) leads to (4.2.12) and (4.2.13). 

• The uniform regularity estimate (4.2.14) is obtained in a very similar way. Of course the regularity 
of the correction established previously shows that its contribution to V • Us^at converges to zero as e 
goes to in the sense of smooth functions rapidly decaying with respect to xi. Therefore the only 
point to be checked is that 

V • ic ( - j $Ar j is uniformly bounded in L^{[Q,T];L°^{Ii x T)), 

which is obtained as the regularity property stated in Proposition 3.6.1, remarking that the L°° bound 
comes from estimates which are stable by the truncation J^ and by conjugation by the semi-group C 
Indeed for almost all /3, provided that 

(n^ + np)$° e hi for some a > -, 

any weak solution $ to {SWq) satisfies the following estimates 

11^ • *'IIli([0,T];L~(RxT)) ^ '^T 

as well as 

V • (c (-^ $jv] < Ct, 

where Ct depends only on T G R^, ||$"||l2(rxT) and JKlIx + np)<i)°||/fcv (neither on N nor on e.) 

• It remains then to establish the equation satisfied by (%.Ar, Ug^^r). A direct computation provides 

A^i^t + A 

X^fi, 

ixe& ,iiJ.e& j\[ 

Egif(A-M-/I) 

iX€.& ,^^^,^fJ.G& j<^ 

By (4.2.10) we infer that dt^N is smooth and rapidly decaying (recalling in particular that IIoQl = 
due to Proposition 2.4.5 page 26), and thus the last two terms go to zero as £ ^- (for all fixed N). 
The previous identity can therefore be rewritten 
(4.2.18) 

edtcj^N = -C (-^A (Q - Ql) U (^ ^N,C (^ ^n) + i^C (-^-) (A' - A^)/: (^^ <Pn + r,,N 

where 

(4.2.19) VfceN,ViVeN,VQeR[X], Ihn ||Q(a;i)r,,Ar|Uoo([o,T];CHRxT)) = 0. 
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Now let US recall that 



dt{C{-]<pN 



11 £(ll*.. 



(4.2.20) 



iQ 



U)\C{-]<^N,c(-]'^N 



+ (/d- jAr)Qi(£( JU,/:Q)$ 



C[^UN,c(^UNJ-l^A'Cr-\<PM 






<p' 



We thus have, recalling that 



{Vs,N,Us,n) = -C - ($Ar + e(j)N), 



(4.2.21) 

dt{V£,N,Us^N) + -L{V6.N,Ue,N) + Q {{Ve,N, U^.n), (»7e.Af , "e,Af)) - '^A'{r]^N, Me.jv) 

£ ' ' ' 

= (/rf - Jn)Ql U (^ $, c (^) $") + g^ (c (^) (a>Ar - $),£ (i\ (a>^ + $) 

+ eQ (c r-j ^N, C r-\ (2$Ar + £0w) j - euA' (c (-\ c^n\ + VeM 

Note that the regularity estimates on (f>jv a-nd (/iat allow to prove that the two last explicit terms in 
the right-hand side go to zero as e ^ (for all fixed N), and therefore to incorporate them into the 
remainder r^jy. 

The stability of the limiting filtered system (S'VFo) allows to prove that the second term in the right- 
hand side of (4.2.21) goes to zero as A^ ^ oo uniformly in e. We have indeed 

Ql[C[-\ ($Ar - $), £ {-\ ($N + $; 

^/ V^/ / L2(RxT) 

and recalling that only the Kelvin waves can have resonances, 



\Ql ($7V - $, $Ar + $) ||l2(r,xT) 



WQl ($Ar - $, $Ar + $) ||l2(r.xT) < \\Ql {^k{'^N - '^),IIk{'^N + $)) |1l2(RxT) 

+ \\Ql (no($Ar - $),n_L($Ar + $)) ||l2(RxT) + WQl (Hx ($Ar - $), no($Ar + *)) ||l2(RxT), 

so by Proposition 3.3.3 page 34 and two-dimensional product rules on the Kelvin part (recall as in the 
previous chapter that i/^ and iJ£ spaces coincide in the case of Kelvin modes) we infer that 

(4.2.22) ||gL($Ar-$,$Ar + $)|U2(RxT) < C„ ^^(^W - $)|| ^= + 1 ||nx($Ar + $)|| ffg 

+ C||no($Ar - $)||L2(RxT)l|ni($Ar + $)||ffi 

+ q|no($jv + *)llL2(KxT)l|ni($jv - mni- 

So by (4.2.8) and (4.2.9) we conclude that 



N 



lim WQl (^jv - $, $Ar + $) ||l2([o,t];L2(rxT)) = 0. 



Let us estimate the first term in the right side of (4.2.21). We can write as above (recalling 

that g(no-,no-) = o) 



Q.[L[i\^,L(i^^ 



< WQl (nK$, Hk^l- + 2 WQl (no$, n^*) 



L2 



L2 



so we find that 
(4.2.23) 



QL[C[^U^,Cr-]<P^ 



L2 



<C lln^/: - $Ar||^i +||no$Ar||i, 
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and thus 



^lim^SUp||(/d- JAr)OL Uf-j $,/:(-)$) ||l2([0,T];L2(RxT)) =0. 



Note that in the case when /? belongs to Af (Theorem 5), equations (4.2.22) and (4.2.23) must be 
adapted using the third estimate of Proposition 3.3.3. 

FinaUy to prove that for aU iV e N, for all T > 0, the quantity 

eQl U (^) ($Ar - $), £ (^) ($Ar + $)") + e{Id - Jn)Ql U (^-) $, C (^) $ 

goes to zero as e goes to zero, in the space L°°{[Q,T] x R x T), wc simply notice that 

Ql U (^ ($Ar - $),£ (^ (a>Ar + $)") + (Id - Jn)Ql U Q <^,C (^) $ 

= JnQl U (^] ($V -$),/: (^] ($V + $)") + (Id - Jn)Ql U (-^ ^N, C (^-^ $Ar 

and the convergence result is obvious if one considers the right-hand side, simply because those terms 
are smooth for each fixed N . 

Combining all the previous estimates shows that (tj^^n, Ue,N) satisfies the expected approximate equa- 
tion (4.2.15), where Re,N satisfies the expected estimate (4.2.16) as well as 

(4.2.24) lim \ime\\Re,N\\L^{[a.T]xRxT) = 0- 

Proposition 4.2.2 is proved. D 



Equipped with that result, we are now ready to prove the strong convergence theorem. The method 
relies on a weak-strong stability method which we shall now detail. We are going to prove that 

(4.2.25) lim lirn\\{Ve,Ue) - {V6.N,Ue,N)\\L^{[0.T]xRxT) ^ 0, 

where {r^cN, Uc^n) is the approximate solution to {SW^) defined in Proposition 4.2.2. Note that combin- 
ing this estimate with the fact that (??e,Ar,Ue.Ar) is close to C (-) <i> provides the expected convergence, 
namely the fact that 



VT > 0, lim 



{Ve,Ue)-C[*-]<i> 



L2([0,T]xRxT) 



The key to the proof of (4.2.25) lies in the following proposition. 

Proposition 4-2. 3. — There is a constant C such that the following property holds. Let {rf^vP) 
and (r]'^,u'^) satisfy assumption (4-2.2), and let T > be given. For all e > 0, denote by {r]^,u^) 
a solution of (SWe) with initial data (77°,^"). For any couple of vector fields (rj,u) belonging 
to i°°([0, r]; C°°(R X T)) and rapidly decaying with respect to xi, define 

£e{t)^ - f {{T],-rif + {l + eij,)\ue--uf){t,x)dx + i^ f f \\/{u, ~ u)Wt' ,x)dxdt' . 
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Then the following stability inequality holds for all t € [0, T] : 

fsW<C£e(0)exp(x(i))+^eW 

+ C I e'^(*)-'^(*') f (dtri+-V-2L+y-{m)jiv-Ve){t\x)dxdt' 

+ C I eX(*)-'^(*') /"(1+er/e) (dtu+-{(3xiu^ + Vri) + {u-V)u-v/\u\ ■ {u- Ue){t' ,x)dxdt' , 
where LOe{t) depends on u and goes to zero with e, uniformly in time, and where 
x{t)^C j (||V-«|U^(RxT) + l|V«||i.(R^T))(i'Mi'- 

Let us postpone the proof of that result, and end the proof of the strong convergence. We apply that 
proposition to (?7, w) = {ri^,jq,u^,j^), where {rj^^N, u^^n) is the approximate solution given by Proposition 
4.2.2. We will denote by Xe.N and Sc^n the quantities defined in Proposition 4.2.3, where (?7,u) has 
been replaced by (rj^.NjUe.N)- 

Because of the uniform regularity estimates on (ry^^r, Ug^Ar), we have 

VT > 0, sup lim [||VUe,Ar||^2([o^T];L2(RxT)) + 11^ ' "£,w|Ui([0:T];L-(RxT)) j < Ct, 



N 

SO we get a uniform bound on Xg^^r : 



SUplim ||Xe,7v||L~([0,T]) < Ct- 



N 



e^O 



Then, from the initial convergence (4.2.2) we obtain that 

ViVeN, £e,7v(0)exp(x£.Ar(i))-^0 as £ ^ in L°°([0, T]). 
Moreover by Proposition 4.2.3 we have 

(4.2.26) dt{r]e.N,u^.N) + -L{ri^,N,Ue,N) + QiiVe.N, U£,iv), {■i]s,n, u^.n)) - iyA'{r]e.N,Ue,N) = Rs,n- 

e ' ' ' 

Let us estimate the contribution of the remainder term. We can write 

' e^-"^-'^-"^*') / R,,M ■ {{Ve.N - Ve), (1 + ef^,)(u,,M - u,)) {f, x)dxdt' = iflit) + I%{t), 
Jo J 

with 

Ie%it) =^ I eJ^^-^*^-^^^-^''^ J R,,NAVe,N~Ve){t\x)dxdt', and 

lj^%{t) <if I* e>^.,«(*)-x.,«(0 J ^^^^(1 + er,,)iu,,^ - u,)(t', x)dxdt' . 
The first term can be estimated in the following way: 

|-^e,Ar(i)| < CT\\Re.N\\L^([0,T];L2(RxT))\\lle,N - %|1l°°([0,T];L2(RxT)) • 

For the second term we can write 

\Ie,NW\ ^ C'tII ^1 + £'ns{Ue,N - "e) |1l°°([0,T];L2(RxT)) II V^l + e%-Re,Af |Ui([0,T];L2(RxT)) • 

Now we can write 



II Vl + S11eRs,N\\h(RxT) ^ ^"(11 -Re^A^ II L2(RxT) + ^lhellL2(RxT) II -R^.W |li4(R.xT))- 

Since 

£||-Re,w|lLi(RxT) ^ £||^e,Af||L=°(RxT)||^e,Af|lL2(RxT), 
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we infer that the quantity e^Ri,,N goes to zero as e goes to zero and TV goes to infinity, in the 
space L^([0,T];L^(R x T)), so in particular 

lim hnie2||i?5^^||^wr(,rp,.^4(RxT)) =0. 

FinaUy by the uniform bound on % in L°°([0, r];i^(RxT)) and by the smallness assumptions on Re,N, 
we deduce that 

' eX-"(*)-x-"(*') / R,^N ■ ((r^e.A - Ve), (1 + sv,){u,,N - u,)) (f, x)dxdt' 
J 

1 



< 2^\\Ve,N - '7elli~([0,T];L2) + \\^/l + erje{Ue^N - "e) lli-([o,T];L2)) + t^£,Ar (0, 

where 

lim lini||a;e,Ar(t)||ioo([o,T]) =0. 

We now recall that by Proposition 4.2.3, using (4.2.26), we have 

£e,Nit) < C£e,N{0)eJipiXe.N{t))+LU,{t) 

+ cj e>^-"W">^-"(*') / i?e,A • {{Ve,N - Ve), (1 + ef],){u,,N - u,)) (f, x)dxdt' 
where 

^e,N{t) = i (||(r/e - Tle.N){t)\\l2 + \\^JT+7^,{u, - U^.j^mWl^) + U j \\V {u, - U e, n) {f )\\l2 {f ) dt' . 

Putting together the previous results we get that lim lim £^^]s[{t) = uniformly on [0, T], hence that 

lim lim |J77e,Ar - ?7e||L~([O.T];L2(RxT)) == 0, 



lim lim\\^l + er]e{ue,N - '"£)||l~([o.t];L2(rxT)) =0, 

A— >CXD £— >0 VI ' J 

By interpolation we therefore find that 

lim lini (11%. at - rie\\L^([0,ThL^(TLxT)) + \\Ue,N - '"s|lL2([0,T];ffi(RxT))) = 0, 

hence (4.2.25) is proved. 

To conclude the proof of Theorem 6 it remains to prove Proposition 4.2.3. As the energy is a Lyapunov 
functional for (SWe), we have 

£e{t) - £e{0) < Ij^J (illt - W^) + (1 + ^Ve){l\u\^ - U ■ U,)\ {f , x)dxdt' 

+ v{Vu - 2Vue) • Vu(t', x)dxdt' 

^11 {9tViv - Vs) + (1 + £Ve)dtU ■ {u - Me)) {t' , x)dxdt' 

dt-qe'n_+ dt{{l + e-qe)ue) 'U- ^^tVelu.l'^)) {f , x)dxdt' 
/ 1^ (Am • {u — Ue) — Aue ■ u) (/:', x)dxdt' . 
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Using the conservation of mass and of momentum we get 

£eit) -S-M< I I {dtviv - Ve) + (1 + e7ie){dtu - lyAu) ■ (u - u^)) (t', x)dxdt' 
t 



J Jl"^- {{l+eve)ue) (iZ - |l^l') {t\x)dxdt' 





/ (- "^^{Pxiuj- + V77£) + V- {{I + ef]e)ue (E) Ue) ] ■u{t',x)dxdt' 

+ / £vrj^/S.u- [u — u^){t' ,x)dxdt' . 



Integrating by parts leads then to 

(4.2.27) 



Se{t)-Ee{Q)< f f (dtfl+^y-U + ^-imn{v-Ve){t',x)dxdt' 

+ hl + STje) idtu+ -(/3a:iu-^ + Vry) + {u-V)u-iyAu) ■ {u-u^){t\x)dxdt' 

- / I {l + er]e)Du:{u-Ue)'^^{t',x)dxdt' 

- [ 2^??^ -li+iR- %)V • (7?m) + Veu ■ VrA (f, x)dxdt' + R„ 



where 



Re(t) = / I evrj^Au- {u — Us){t' ,x)dxdt' . 



(4.2.28) 



/o 
The last term is rewritten in a convenient form by integrating by parts 

— ?7g V ■ u+ {rj — rj^y^ ■ (rju) + rjgU ■ Wrj J (t' , x)dxdt' 

-rjl^ ■ u+ [rj — 'rj^){u ■ Vrj + ryV • u) + rj^u ■ V77 j (i', x)dxdt' 

/ ( -rJlV -u+iri- rie)r]y ■ u+ -u- Vrf ) (t', x)dxdt' 

* /• 1 

/ oiVs ^ v) ^ ■ u(t' ,x)dxdt' . 

Plugging (4.2.28) into (4.2.27) leads to 

£,{t)~£M< I f (dtri+^V-u + V-{riu)j{ri-Ve)it\x)dxdt' 
(4.2.29) + f f{l + er],)(dtu+-{Pxiu-^ + \7Ti) + {u-V)u-i^Au\ ■{u-u^){t',x)dxdt' 

- / {l + er)e)Du:{u-u^)^^{t',x)dxdt'- / -{t]^ - ^f^ ■ uit' ,x)dxdt' + Reit). 



4.2. STRONG CONVERGENCE OF FILTERED WEAK SOLUTIONS TOWARDS A STRONG SOLUTION 77 



In order to get an inequality of Gronwall type, one has to control the right hand side in terms of £e 
We start by estimating the flux term. We have 

(1 + er]e)Du : {u - Ue)^^{t', x)dxdt' 



< / (||Vu||l2(r.xT) +e||?7£|U2(RxT)||Vu||L°o(RxT)) \\u- Us\\\ir^^T,Jt')dt' 
JO 

<C I (||Vu||l2(RxT) +e||f?e||L2(RxT)||VM||L^(RxT)) ||w-'"e|lL2(RxT) 

Jo 

X \\u- Ue\\H^(^-R.y,T)i^')d't' 

and 

lll^-"ellL2(RxT) < |lv/l + e%(Ue-u)|li2(RxT) + ehelU^lRxT) II " - "ell L2(rxT) 11" - "e llffi (RxT) 

which implies 



lll*-'«ellL2(RxT) < 211^/1 + e?7e(Ue-u)||^2(RxT) + 16e^ he II L2(RxT) 111* " '«e|||i(R,xT) ' 



Therefore, using the uniform bounds on 77^, ^1 + erj^u^ and on Ug given by the energy estimate, we 

gather that 

ft 

{l + eT]e)Du: {u-Ue)'^^{t',x)dxdt' 



<C / {\\yu\\L^ + e\\Wu\\L^)\yi + er^eiue ~ u)\\lA\u- ue\\m{t')dt' 

(4.2.30) -^o 

+ Ce / {\\Vu\\l2 + e\\Vu\\L^^)\\u- ue\\%,{t')dt' 
Jo 

< ^ /"||w-ws|||i(t'yi' + - /"l|Vw||i2||Vl + £%(«e-ll)lli2(t')c^^' + ^sW- 

We also have 

(4.2.31) -y" j\{ile-rifV-u{t\x)dxdt' <]^j ||V • 2i|U~(RxT)h - %||i2(RxT)(i')rfi', 
so we are left with the study of the remainder i?^. We have 

Re{t) < £'^ll»7ellL=°(R+;L2(RxT)) / II Am||l4(rxT) || " - "e IIl'I(RxT) (i')c^i'- 

"'0 



'0 

The above estimate on ||u — Me||L2(-RxT) implies in particular that ||u — Ue||i2(RxT) is bounded 
in i^([0,T]), hence we get that ||m — ■Ue||i4(RxT) is also bounded in L^([0,T]). So we infer directly 
that Re{t) goes to zero in _L°°([0,r]) as e goes to zero. That result, joint with (4.2.30) and (4.2.31) 
allows to deduce from (4.2.29) the following estimate: 



^£e{t)~£AO)< I f (dtv+^y -u + y ■ imn {V-Ve){t',x)dxdt' 



2 

rt 



+ / / (1 + e?7e) [ atM+ -(/3xiu^ + Vry) + (u • V)m- J/Au j • (u-Ue)(t',a:;)dxdi' 

+ ^J l|Vli||i2|| Vl + £VeiUe - u)\\h(R^T)it')dt' + \j \\V ■U\\LAV-Ve\\l2{t')dt' + uS) 

thus applying Gronwall's lemma provides the expected stability inequality. Theorem 6 is proved. D 
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4.3. Strong convergence of filtered weak solutions towards a weak solution 

The aim of this section is to prove an intermediate convergence resuh, in the sense that we will seek 
a strong convergence result of the filtered weak solutions, towards a weak solution of the limit system; 
thus no additional smoothness will be required on the initial data other than L^(R x T). As explained 
in the introduction of this chapter, the lack of compactness in the spatial variables of 77^ will prevent 
us from obtaining at the limit the expected system {SWq): a defect measure remains at the limit, 
which we are unable to remove. In order to gain some space compactness and to get rid of that defect 
measure, we propose in the final paragraph of this section (Paragraph 4.3.4 below) an alternate model 
which takes into account capillarity effects, and for which one can prove the strong convergence of 
filtered solutions towards a weak solution of (SWq). 

The first result of this paragraph is the following. 

Theorem 7 (strong convergence towards weak solutions). — Let (77", u°) e L^(R x T) and 
(77°, u^) satisfy (4.1.1). For all e > 0, denote by {r/^,Ue) a solution of (SWe) with initial data (77", w^), 
and by 

Up to the extraction of a subsequence, $e converges strongly in Lf^^{'R'^; Hf^^(R x T)) (for all s < 0) 
towards some solution $ of the following limiting filtered system: for all lA G 6, there is a bounded 
measure v\ G A^(R+ x R x T) (which vanishes if X = 0), such that for all smooth ^^ G ¥^ci{L — i\Id), 



I ^■^l{x)dx-v I I A'^^-^l{t',x)dxdt' 



+ f fQL{^,<5)-nit',^)dxdt'+f fw-{^iyv^{dt'dx)^ J<5°-^l{x)dx, 
where Ql and A'^ are defined by (2.3.6) page 19, and where <^^ — [if^u^). 

Remark 4-3.1. • Note that, by interpolation with the uniform L^q^(R+; iJ^(R x T)) bound on u^, 
we get the strong convergence of u^ in i^|^^(R+; i^(R x T)) ; up to extraction of a subsequence. 



£|iU~' 



^ m LL(R+). 

L2(RxT) 



• As explained above, the presence of the defect measure v\ at the limit is due to a possible defect of 
compactness in space of the sequence {rig)gyQ. As the proof of the theorem will show, that measure is 
zero if one is able to prove some equicontinuity in space on rj;,, or even on erj;,. Since we have been 
unable to prove such a result, we study in the final paragraph of this section a slightly different model, 
where capillarity effects are added in order to gain that compactness. Note that the model introduced in 
Paragraph 4-3.4 is unfortunately not very physical due to the particular form of the capillarity operator 
(see its definition in (4.3.9) below). 

Theorem 7 is proved in Sections 4.3.1 to 4.3.3, and the result in the presence of capillarity is stated 
and proved in Section 4.3.4 . 
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4.3.1. Strong compactness of IlA^e. — Let us prove the following lemma. 
Lemma 4-3.2. — With the notation of Theorem 7, the following results hold. 

• For all iX e &\ {0}, Hx^e is strongly compact in L^([0, T]; iJ''(R x T)) for all T > and all s G R; 

• no$e is strongly compact in L^{[0,T]; Hi^^(R x T)) for all T > and all s < 0. 

Proof. — • For all A 7^ 0, we recall that by Proposition 2.2.3 page 12, the eigenspace of L associated 
with the eigenvalue i\ is a finite dimensional subspace of H°° (R x T) . Therefore the only point to be 
checked is the compactness with respect to time, which is obtained as follows. 

Let {n,k,j) e N x Z x {—1,0,1} be given, such that A = T{n,k,j) ^ 0, and let '^n,k,j be the 
corresponding eigenvector. Multiplying the system (SW^) by ^n.fcj (which is smooth and rapidly 
decaying as |xi| goes to infinity) and integrating with respect to x leads to 

dt / {vei^n.k.j)o + m, ■ 'I-; fe^^) (t, x) dx + '^^"^ '-^^ / (?7e(*„,fcj)o + m, ■ *^,fc,j-)(t, x) dx 

+i^Jwu, : V*;_fc_^-(t,a;) dx - J m, ■ {u, ■ W)^'^^,,^j{t,x) dx- ij ry^V • ^'^^^^ dx = 

where ^n.fc.j denotes the complex conjugate of "^n.k.j, or equivalcntly 

dt Uxp nJli^lihlJlj /"(77,(§„,fc,,)o + me • K^k^j)it,x) dxj 

(4.3.1) +^|v(^exp(^^^^i^^)«,) :V<,,^,.(i,x)dx 

_ J ^^p ^ »tT(n,fc,.7) ^ ^^^ ^^^ V)§;,fc,, + ^T^^V • §;,,,,.) (t, x) dx = 0. 

By the uniform estimates coming from the energy inequality we then deduce that 

/ f itrln k j)\ f - - \ 

dt exp ^ — ^ — j / (?7e(^n,fcj)o + ^T^s ■ ^n k j){t,x)dx J is uniformly bounded in e. 

Therefore the family 

exp ( — I Il\{rie,ms) ] is compact in i^([0,T]; iJ''(R x T)) for anys G R, 
V '^ / /e>0 

and since erj^u^ converges to in L^(R+; H^{Yi x T)) for all s < 0, we deduce that 

exp ( — j Iix{Tle,u^) = nA<i>e is compact in L^([0, r];i7''(R x T)) for anys e R. 

• For Ilo'I'e = no(77£,Me) the study is a little more difficult since the compactness with respect to 
spatial variables has to be taken into account. By the energy estimate we have the uniform bound 

$e is uniformly bounded in Lfo^(R+, L^(R x T)). 
We recall that we have defined in Section 3.1 (Definition 3.1.1 page 30) the space 

Hl = U^L\B.xT)/ Y. (l + n + A:')^(V'l*n,fc,,)i2(R,T)<+^ 
I n,k,ji£S 
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where 5 = N x Z x { — 1, 0, 1} or equivalently (see Proposition 3.1.2 page 30) 

Hi = U e l2(R X T) / (Id - A + (3^xly/^iP e L^(Ti X T)| 

As (^n,o,o)neN IS a Hilbertian basis of KerL, we have for all T > and all s < 
J2 (^sl*".0.0)i2(RxT) *«.o,o - no$e 



n<N 



^ as iV ^ oo uniformly in e. 

L2([0,T]:HJ) 

Let il be any relatively compact open subset of R x T. Proposition 3.1.2 page 30 implies that, for 
all s > 

and conversely for s < 0, 

(4.3.2) i?"(Rx T) Ci?i C i/"(0). 

Thus for all s < and all T > 0, we have 



as A^ ^ oo uniformly in e. 



J2 (*e|*n,0,o)i2(RxT)*«.0,0-no$e 

L2([0,T];fl-=(n)) 

Moreover the same computation as previously shows that for any n e N, 

dt ( / {Vefinfifi + me ■ M„,o,o)(t, x)dx ) + >^ Vm^ : Vu„,o,o(t, x)dx 
(4.3.3) ^"^ /J 

- nie ■ {ue ■ W)unfifl{t,x)dx = 0, 

and, since erj^u^ converges to in L^(R+; H^iJEL x T)) for any s < we get 

2_j nn,o.o('7e, We) IS Compact in i^([0,r] x R x T). 

n<N 

Combining both results shows finally that 

no$£ is compact in L2([0,r]; 7?f„^(R x T)) 
for all T > and all s < 0. Lemma 4.3.2 is proved. 



D 



As the spectrum of L, 6 is countable (see Chapter 2), we are therefore able to construct (by diagonal 
extraction) a subsqucncc of $£, and some $a G Kcr(L — iXId) such that for all s < and all T > 

(4.3.4) Va e 6, nA$e ^ $A in ^'([O, T]; Hf„,(R x T)). 

Note that the $a defined as the strong limit of IIa^s can also be obtained as the weak limit 
of exp (— ) (r/g, Ug). We have indeed the following lemma. 

Lemma 4-3.3. — With the notation of Theorem 7, consider a subsequence o/($e)£>o, and some $a 
in Ker(L — iXId) such that for all s < and all T > 

VzA e 6, nA*e ^ 3>A in L2([0,T];i/f<,,(R x T)). 

Then, for all iX G &, e~^ (jig^u^) converges to $a weakly in L^([0,T] x R x T). In particular, for 
alliXe 6, the vector field <^'^ belongs to L'^{[0,T]; H^(R x T)). 



4.3. STRONG CONVERGENCE OF FILTERED WEAK SOLUTIONS TOWARDS A WEAK SOLUTION 81 



Proof. Denote by (,.,..) any weak limit point of the sequence exp {^-^) (,.,..) (recall that 

the sequence is bounded in L^^^(R+ x R x T)). Let x and ij^ be any test function and vector field 
in I?(R+ X R X T). Multiplying the conservation of mass in {SW^) by ex exp (^i^) and integrating 
with respect to all variables leads to 

erie exp f — j f dtX H X j + (1 + sVehe exp f — j • Vx j dxdt = 0. 

Because of the bounds coming from the energy estimate (1.4.5) page 6, we can take limits in the 
previous identity as e goes to to get 

{u\ ■ Vx + i\rj\x)dxdt = 0. 

Similarly, multiplying the conservation of momentum by eV'exp ( — \ and integrating with respect to 
all variables leads to 

fe(l + e?7e)ueexp f — j {dti}}^ 1/" j +£cxp f — j (1 + ei^e)ue ■ [u^ ■ V)ip 

+I3xi exp f — ) (^ + £%)«£ • -0^ + (1 + -%)% exp f — j V- ■0 

fit\\ ^ \ , , 
+evu^ exp I — I • i\ip \ dxdt = 0. 

Once again the bounds coming from the energy estimate (1.4.5) will enable us to take the limit as e 
goes to 0, to get 

(?7aV • tj) + [3xxU\ ■ ip + iXu\ip)dxdt ~ 0. 

It follows that {ri\{t),u\{t)) belongs to Ker(L — iXId) for almost all t E R+, and we conclude by 
uniqueness of the limit and L^ continuity of 11^ that $a = {v\t'^\)- The lemma is proved. D 



4.3.2. Strong convergence of <&£. — As a corollary of the previous mode by mode convergence 
results, we get the following convergence for <!>£. 

Lemma 4-3. 4- — With the notation of Theorem 7, the following results hold. Consider a subsequence 
of ($e); and some ^\ e Ker(_L — iXId) such that as constructed in (4.3.4), for all s < and all T > 

VzA e e, nA$e ^ $A in L2([0,T];i/f<,,(R x T)). 

Then, 

$e ^ $ = ^ $A weakly m Lf^^{Il+;L^{Il x T)), 
iAee 
and $£ ^ $ strongly in Lf^^{K'^-,Hi^^(R x T)) for all s < 0. 

Moreover, defining Kn as in (3.4-1) page 40, we have for any relatively compact subset Q o/R x T, 
for all T > and for all s < 0, 

(4.3.5) ||(W- is:Ar)$£||i2([o,T];H=(!:2)) + Wild - KN)C{-)'i>,\\L2([a,T]:H^{n)) -^0 as N^oo, 

uniformly in e. 
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Proof. — The first convergence statement conies directly from tlie uniform bound on <i>g in the 
space Lf^^(R'^;L^(R x T)) and the L^ continuity of Ux- 

In order to establish the strong convergence result, the crucial argument is to approximate (uni- 
formly) <i>e by a finite number of modes, i.e. to prove (4.3.5). The main idea is the same as for the 
approximation of Ha^e in Lemma 4.3.2. We have for all T > and all s < 



Y. (*",fc J l^e)i2(RxT) *",fc,J - ^e 

n<N,\k\<N 



as A^ ^ oo uniformly in e, 



Lmo,T];Hl 



and similarly 



E 



-iT{n,k,j)l 



(*n,fcjl* 



3\^^'L^{-RxT] 



*n,fcj - 'C(-)$e 



n<N,\k\<N 



as A^ ^- OO, 



L2([0,T];H£) 



uniformly in e. Therefore for all relatively compact subsets O of RxT, the embedding oi HI into H''{il) 
recalled in (4.3.2) implies that both quantities 



J2 (*n,fc,jl*e)i2(KxT)*»,fc,J 



$. 



■i<N,\k\<N 



and 



2(RxT) ^n,k,j 



/:(-)$s 



n<N,\k\<N 

converge strongly towards zero in L^([0, T]: H^(fl)) as N goes to infinity, uniformly in e. Finally (4.3.5) 
is proved. 

The strong convergence is therefore obtained from the following decomposition: 

$£ - $ = (/rf - XAr)$e + KN{<^e - *) - (W - Kn)<^ 

The first term converges to as A^ ^ oo uniformly in e in L^^^(R+, i/f^^(R x T)) for all s < 
by (4.3.5). By Lemma 4.3.2, the second term (which is a finite sum of modes) converges to as e ^ 
for all fixed N in Lf^^(R+;H%K x T)) for all s < 0. The last term converges to as A^ ^ oo 
in i^^^(R+; i^(R x T)). Thus taking limits as e ^ 0, then as A^ ^ oo leads to the expected strong 
convergence. D 



4.3.3. Taking limits in the equation on n^^e. — The next step is then to obtain the evolution 
equation for each mode $>, taking limits in (4.3.1) and (4.3.3). In the following proposition, we 
recall that the first result (concerning the geostrophic motion) relies on a compensated compactness 
argument, i.e. on both the algebraic structure of the coupling term and the particular form of the 
oscillating modes, which implies that there is no contribution of the equatorial waves to the geostrophic 
flow. That result was proved in Section 4.1 (see also Proposition 2.4.5 page 26). Here we will prove 
the second part of the statement, concerning the limit ageostrophic motion. 

Proposition 4-3.5. — With the notation of Theorem 7, there is a subsequence o/ ($e) such that the 
following result holds. Consider a family {^\)i\ee such that $a G Ker(L — iXId) and such that for 
all s < and all T > 

VzA e 6, nA$e ^ *A m L2([0,T];i/f<,,(R x T)), 



as constructed in (4-.S.4)- 
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Then, $o = (^Oi^o) satisfies the geostrophic equation : for all (77*, u*) belonging to KerL and satisfy- 
ing u* e i/^(Rx T), 

('7o'7* + uo,2'U'2){t, x) dx + v I / Vuo,2 • ^^u^it' , x) dx dt' = / [rfrj* + u^U2){x) dx. 

Moreover for A 7^ 0, $a = (^A'^a) satisfies the following envelope equation : there is a measure v\ 
in 7V{(R+ X R X T), such that for all smooth $^ = (<I>^ „, i^xY) € Kcr(i - iXId), 



(^x-<^*x{t,x)dx + iy jV'^\:V{'^X)'{t\x)dxdt'+j I V ■ {-^D' v x{dt' , dx) 

+ Y. f Q('^f.,'^^)-n{t',x)dxdt'=. f<^'>-^l{x)dx, 

X = H + ^i 

where Q is defined by (2.3.4) P<^9^ ^8. 

Proof. — Let us first recall that for A 7^ 0, Kcr(L — iXId) is constituted of smooth, rapidly decaying 
vector fields, so that it makes sense to apply Ilx to any distribution. 

Starting from (4.3.1) we get that for all smooth $^ = ($^ q, ($^)') G Ker(i - iXId) 

exp (^^) (,7s$io + ^e ■ inrnt, ^)dx - Jiv^eKo + ^° • {ny){^)dx 

(4.3.6) +1^ f f\/(exp (—^ u,\ : ^{^D' it' ,x)dxdt' 

- j'j exp (^^^ (m, . [u, ■ V)(|.l)' + iry^V • ($^)') (f, x)dxdt' = 0. 

Taking limits as e ^ in the three first terms is immediate using Lemma 4.3.3 and the assumption on 
the initial data. The limit as e ^ in the two nonlinear terms is given in the following proposition. 

Proposition 4-3.6. — With the notation of Proposition 4.3.5, we have 

1^ I exp (^^) m, . (u, ■ V)(#l)'(i', x)dxdt' ^ fj ^ $;,•($;,• V)($l)'(t', x)dxdt' , 



^^^-^i=A 



and 



i r /exp (—\ ^IW ■ {^iy{t',x)dxdt' -^\f I E *M.o$A.oV • {^iy{t',x)dxdt' 



/j, + /j, = A 



t 

V • i^lYvxidt'dx). 




Before proving that result, let us conclude the proof of Proposition 4.3.5. It remains to check that 
-/ E (*;,-(*;i-V)(l>I)' + $^,o*A,oV-(il)')dx 



^A,^^^,^^^G © 



/ E ((i>;-v)$;,-(i>i)' + v.(<i>,,o<i';,)i>io)dx, 

J ,.-Lr. — \ 



IJ, + IJ, = X 



84 CHAPTER 4. CONVERGENCE RESULTS 



Since 

/j.-\- jl^ X iJ,-\- p.^ X 

Clearly one has 

so since /z and p, play symmetric roles, we just need to check that 

fi + il = X ^ -^ /^ + ,i = A "^ 

Recalling that 
we have 



V • {<^^,o<^';,)n.,o dx = i(p +p) j $^,0$;^ • my dx+j $^.0/32:1$;^ • iny dx. 

Then we write 

so that 

J V • {^^,o<^'^)n.o dx^ JV- $;$;^ • (l.^)' dx + ifij $^,o<i>;i • (l>D' dx 

- J <Pf,.ol3x,<P'^^ ■ i^iy dx. 
Exchanging the roles of /i and /t in the first integral we get 

E [^■{<^f.fi^'f.)Kodx= J2 (/v •$;,$;.($*)' da; 

+zji J <i>^,o^'^ ■ iny dx- J $^,o(.A$A - ^*M-o) • iny dx) , 






so 



E [^■i'^^,o<^'f.)Kodx^ E {[^■'^'f.^'^-inydx+ f<^^,oy<^f,^o-{nydx 



/^ + /^^A /^-j-/^^A 

The result finally follows from the fact that, by symmetry, 



E /*M.oV<i>A,o-(*I)'rfx = i J2 /v(a>^,o<i'p.o)-($A)''^^ 



(1 + fi^ X /j + ^i^A 



which finally implies that 



E ( /<i>;,-(«'l)'V-$;, da; -i/'^^^o'I'A^oV •($!)' dx) - ^ f \/ ■ {^^.,^'^}^lo dx. 

u. + fi = X "^ -^ M + A = ^ 






Now let us prove Proposition 4.3.6. The idea is to decompose $e on the eigenmodes of L, by writing 



ixee 
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Note in particular that by (4.3.5), for any s < 0, 

{T^,,u,){t)-c(^^ Kn^S) ^ in LL(R+;i?L(R X T)) 
as N goes to infinity, uniformly in e. Let us also introduce the notation 

^e,N = Ci--j{T]e,N,Ue,N)^KN^s, and 



We will start by considering the first nonlinear term in Proposition 4.3.6, namely 

f /"cxp ('-^^ me ■ (u, ■ V)i^iy{t',x)dxdt'. 
We can notice that 

I exp(- — \me-{ue-V){^iy{t',x)dxdt'^ / cxp(- — \ e-q^Ue ■ {ue ■ \7)i^*x)' {f ,x)dxdt' 

'"exp (—\ u, ■ {ue ■ \/)my{t',x)dxdt'. 



IQ J \ £ / 

The uniform bounds coining from the energy estimate imply clearly that the first term converges to 
as £ ^ 0. Then we can decompose the second contribution in the following way: 

/ /exp r-^j u, ■ {ue ■ \/){^iy{t',x)dxdt' 
^ exp ( — ) u, ■ {u, ■ \/){^iy{t',x)dxdt' 

Jo JiR\l-R.,R])xT \ £ / 

(4.3.7) + / / exp ( — ) (we - Ue,N) ■ {ue ■ V)(l>^)'(t', x)dxdt' 

Jo Jl-R.,R]xT V £ / 

+ / / exp ( I Ue.Tv • {{ue - u^.m) ' V) {^\y{t' , x)dxdt' 

Jo J[-R,R]xT \ £ / 

+ exp(- — ) u,,N ■ (ue.M ■ V)(l>^)'(t', x)dxdt'. 

Jo J[-R,R]xT \ £ / 

Let us consider now all the terms in the right-hand side of (4.3.7). The uniform bound on Ug and the 
decay of (f>^ imply that the first term on the right-hand side converges to as i? ^ oo uniformly in e. 

By the inequality 

/ exp ( ^^ — ) {u, - u,,n) ■ {u, ■ V)(l'^)'(t', x)dxdt' 

J[-R.,R]xT \ £ / 

< C\\Ue - Ue,Jv||L2([0,T];H=([-_R,i?:]xT))||We||L2([0,T];Hi(RxT))||$AllM'2.°°(RxT), with - 1 < S < 0, 

we deduce that the third term converges to as A^ ^ oo uniformly in e. 

Now let us consider the third term on the right-hand side. Since u^^n corresponds to the projection 
of $£ onto a finite number of eigenvectors of L, wc deduce that 

VN e N, 3CAr,V£ > 0, \\Ue.N\\L^(R+;H^{IlxT)) < Cn- 
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Thus 

/ exp I Ue.AT • ((«£ - Ue.M) ' y)(^*x)'{t',x)dxdt' 

J[-fl,fl]xT V ^ / 

< CnW^s - Ws.Af ||L2([o^T];H=([-fl,_R]xT))ll*Allw2.°=(RxT) 

and, for all fixed N and i?, the fourth term eonverges to as M ^ oo uniformly in e. 
It remains then to take limits in the last term of (4.3.7). It can be rewritten 

/ exp ( ) Ue.N ■ {us,M ■ V)($^)'(t', x)dxdt' 

J[-R,,B]y.T \ ^ / 

^ exp f 't'(>^-^-f^) \ ($^^^)' . ($^^^^ . V)($l)'(t', x)dxdt'. 
A-R-MxT.^.^r^ee V ^ / 

This in turn can be written in the following way: 

"''^^ " '' " ^^ ^ K.,.N ■ (K.^.M ■ v)($l)'(t', x)dxdt' 



J2 exp 

^[-fl,fl]xT^^^^^gg 



^ exp ( '''^^ j" ^^ ) ($i_^_^ - cl>;^,^) . (ci,;^^ . V)($l)'(t', a:)dxdt' 



-'[-fl,i?]xT^^_^^gg 

^ exp f '^'^^'^'^^ ) $;^ . ((d,;^ ,, - $;,^) . V)($l)'(t', :.)d:rdt' 
Jl-R-MxT^f,,r^ee \ ^ / 

^ exp f 't'{\-f.-f^) \ ^^^ , ^^^^^ , V)i^lYit',x)dxdt'. 

-'hfl^^JxT^^^^^gg V e / ^ 

We have denoted 

^M.Af = n^'^w, where $Ar = i^Ar$. 
The first two terms on the right-hand side go to zero as e goes to zero, for all given N, M and R, due 
to the following estimates: for —1 < s < 0, 

E K'^Un - k.n) ■ iK,^,M ■ v){nnt',x)\dxdt' 

< CN,M\\^'e^N ^ ^'N\\L^i[0,T]:H'{[-R.,R]xT))\\^e,M\\L'=^{[O.T]-H^{RxT))\\^*\\\w^-°°i'R.xT)i 

and similarly 

/ / E K-N-i{<^UM-'^kM)-^)innt',x)\dxdt' 

Jo -'[-R.-RlxT^^^^-gg 

< CN,M\\^'e,M ~ *MllL2([0,T];ff«([-fi,fi]xT))ll*i>e,w|lL°°([0,T];Hi(R.xT))ll*Allw2.~(RxT)- 

Finally let us consider the last term, which can be decomposed in the following way: 

^ exp f 't'i^-^^-f^) \ $^^ . ($^^ . S7){nnt',x)dxdt' 
^[~fl,fl]xT^^^,-ge \ s J 

J[-R,B]xT ,^,,,^£6 

^ exp ( ^t'i^-^^-~^) \ ^^^^^ . ($1^ ^^ . v)($*)'(t',x)dxdt'. 

J[-fl,K]xT ,^,,^£6 V £ / 
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For fixed N and M, the nonstationary phase theorem (which corresponds here to a simple integration 
by parts in the t' variable) shows that the second term is a finite sum of terms converging to as 
e -^ 0. And the first term (which does not depend on e) converges to 

ft 



I I Y. <^',-{^'f.-^)m'{t',x)dxdt' 



l_i + l_i = X 



as TV, M, R -^ oo, because $^ converges towards $' strongly in i^([0, T]; L^(R x T)) when N goes to 
infinity, and then by Lebesgue's theorem when R goes to infinity. 

Therefore, taking limits as e ^ 0, then M -^ oo, then A^ -^ oo, then i? ^ oo in (4.3.7) leads to 
^* j exp (^^^ m, . {u, ■ W){nnt', x)dxdt' -^ j'j ^ *;,•($-• V)($l)'(t', x)dxdt'. 

iX,ifj.,iil^& 

Finally let us consider the second term of the proposition, namely 

The first step of the above study remains valid, in the sense that one can write 

1 f /exp f ^) ^,2V . {nnt',x) dxdt' ^\fj exp f ^) r^^V • ($!)'(*', x) dxdt' 

^ Jo J \ ^ / ^ Jo Jll\[-R,R]xT V £ / 

+ 1 f I exp f ^) ryfV • {^l)'{f.x) dxdt', 

^ Jo J[~R,R]xT \ £ / 

and the first term converges to zero uniformly in e as i? goes to infinity, due to the spatial decay of 
the eigenvectors of L. For such a result, a uniform bound of % in iy°°(R+;L^(R x T)) is sufficient. 
However the next steps of the above study do not work here, as we have no smoothness on ry^ other 
than that energy bound. In order to conclude, let us nevertheless decompose the remaining term as 
above, for any integers N and M to be chosen large enough below: 



(4.3.8) 



U f exp('^)r^^,V-m'it',x)dxdt' 

^ Jo J\-R,R]xT V £ / 

exp f ^ j (rye - rie,N)rieV ■ ($^)'(t',a;) dxdt' 
exp ( ^^ — ) rie,N{rie - ??£,m)V • {^D' [f , x) dxdt' 



if 

2 Jo J[-R,R]y.T 

1 /■* /■ fit'y 

2 Jo J[-R,R]xT \ £ , 

1 /■* /■ fit'\\ 

exp f — j Tie^NVe.M^ ■ i^\Y(t\ x) dxdt'. 



2 Jo J[-R,R] 



IxT 



The sequence — exp 1 (77^ — % jv)^ is uniformly bounded in TV G N and e > in the 

2 \ e J 

space Lj^^^(R+ x R x T), so up to the extraction of a subsequence it converges weakly, as e goes to 
zero, towards a measure wa.w, which in turn is uniformly bounded in A^(R+ x R x T). Denoting 
by vx its limit in A^(R+ x R x T) as iV goes to infinity, we find that 

Iff exp (—\ (r/e - %,jv)%V • (l'^)'(t',a;) dxdt' ^- f f V • {^X)'vx{dt'dx) 

^ Jo J[-R.,R]xT \ £ / JQ J[-fl,fl]xT 
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as e goes to zero and A'^ goes to infinity, which in turn converges to 



V-{^l)'v^{dt'dx) 



as R goes to infinity, due to the smoothness of V • ($^)'- 

Note that as 6 is countable, one can choose a subsequence such that for all iA e 6, the se- 
quence — — exp I {rj^ — rie^N)'r]e converges towards v\ as e goes to zero and N goes to infinity. 



Now let us consider the two last terms in (4.3.8). We recall that yy^jv corresponds to the projection 
of $e onto a finite number of eigenvectors of L, so it is smooth for each fixed N. In particular we can 
write, for any s < 0, 



/ exp ) r]e,N{r]e - ?7£.m)V • ($^)'(i',a;) dxdt' 

J[-R,R]y.T \ £ 



< CNWrie — ??e,M||L=°([0,T]:H»([-_R,fl]xT))||'I'Allvi/2,oo(R_xT)- 



So letting M go to infinity we find that this term converges to zero uniformly in e for each fixed N 
and R. 





e 




lt'{\ 


-M- 


-P) 




e 




it'{\ - 


-M- 


-n) 




£ 




it'{\ - 


-M- 


-n) 



Finally for the last term of (4.3.8) we write similar computations as for the first nonlinear term in 
Proposition 4.3.6. We have indeed 

/ exp ( I r?e,Ar%.MV • (^*x)'{t' , x) dxdt' 

J[-R,R]y.T V £ / 

V E e'^P f '^'^^'^'^^ l ('i>M.e,^)o(<i>M.e.M)oV • (l>l)'(t', x) dxdt' 

Y^ exp ( '^ ''' '" ""' ) ($^,e,N - *p,w)o($a,£.m)oV • ($a)'(^'> 2;) dxdt' 
V exp ( "" ^^' "^ "^^ ) ($^.N)o($A.e.M - $a.m)oV • ($D'(t', x) dxdt' 
/ V exp ( "" ^^' "^ '^' ) ($^.w)o($a.m)oV • (^a)'(^', 2;) dxdt'. 

■>'[-«,«] XT, ^^^^gg V e / 

The two first terms in the right-hand side are easily shown to converge to zero as e goes to zero, for 
each fixed N and M. We have indeed 

V exp ^^'(^-A^-A^) j (,1,^ ^^ _ <I>^,Ar)o($A.e.M)oV • ($l)'(t', x) dxdi' 
.^[-fl.fllxT, ,-gg V £ 



< CM||<&M,e,^f ~ ^M,wlU°=([0,T]:ff»([-fl,fl]xT))II^Allw2.°°, 



and similarly 



V exp f '^'^^ ^ ^^ ) (<i>^.w)o($A,e.M - <i>M.M)oV • (l>l)'(t', x) dxdt' 

-^[-fl.flJxT^^^-gg V £ / 

< C'7v||<&^,e,M - *I>/i,M|lL=°([0,T];ff«([-fi,fi]xT))II^Allw2.' 
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Finally the last term on the right-hand side is dealt with by a nonstationary phase argument, and we 
have as above, as e goes to zero and then M, N and R go successively to infinity, 

1 /■* /■ .(-^ (it'{\-n-iiy 



^ , , E ^^P ^ ^-^ (*m,a^)o($a.m)oV • (I.^)'(i', x) dxdt' 



^ Jo J ^ + ^ = ;, 



Proposition 4.3.6 is proved, and therefore also Theorem 7. D 

4.3.4. The case when capillarity is added. — In this final paragraph we propose an adaptation 
to the Saint- Venant model which provides some additional smoothness on 77^, and which enables one 
to get rid of the defect measure present in the above study. The model is presented in the next part, 
and the convergence result stated and proved below. 

4.S.4-1- The model. — Let us present an alternative to the Saint-Venant model studied up to now, 
which presents the advantage of providing the additional smoothness of erj^ which is lacking in the 
original system. Its disadvantage however is that there is no real physical meaning to the capillarity 
operator we use in that model. With the notation of Chapter 1, we choose indeed the capillarity 
operator 

(4.3.9) K{h) = ni-Af^h, 

where k > and a > 1/2 are given constants. After rescaling as in Chapter 1, we find the following 
system: 

dtV + -V- f (1 + er))u\ = 0, 

(4-3-10) dtu + u-Wu+ — u-L + -Vn ^Au + eKW(-Af"T] = 0, 

e e 1 + erj 

V\t=o = ff. "|t=o = '"°- 

In the next part we discuss the existence of bounded energy solutions to that system of equations (under 
a smallness assumption), and the following part consists in the proof of the analogue of Theorem 7 
in that setting. One should emphasize here that the additional capillarity term that is added in the 
system will not appear in the limit, since it comes as a 0(e) term. Moreover it is a linear term, so it 
should not change the other asymptotics proved in this chapter. However its unphysical character (as 
well as the smallness condition on the initial data) made us prefer to study the original Saint-Venant 
system for all the convergence results of this chapter. 

4-3.4-2. Existence of solutions. — The following theorem is an easy adaptation of the result by D. 
Bresch and B. Dcsjardins in [2] (see also [23] for the compressible Navier-Stokes system), we give a 
sketch of the proof below. 

Theorem 8 (existence of solutions in the case with capillarity) 

There is a constant C > such that the following result holds. Let (77^, u^) be a family of H^"^ x i^(R x 
T) such that for all e > 0, 



\j{{rr.f + «£^I(-A)".7.T + (1 + ery,")|K,T) W 



dx < £". 
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If £^ < C, then there is a family (rj^^u^) of weak solutions to (4-3.10), satisfying the energy estimate 
i f{vl + Ke2|(-A)"ry,|2 + (1 + e77,)|w,,|2) {t,x) dx + y j j \Vu,\^{t', x) dxdt' < f «. 

Proof. — Weak solutions can be constructed by a standard approximation scheme obtained by regu- 
larization : compactness on the approximate solutions comes from the a priori bounds derived from 
the energy inequality, which is obtained formally in a classical way by multiplying the momentum 
equation by u, using the mass conservation and integrating by parts. It allows to derive immediately 
the following a priori bounds (denoting by rj and u approximate solutions) : 

rje L°°(R+;L2(Rx T)) 

er/eL°°(R+;ij2«(Rx T)) 

(1 + e??)|u|2 e L°°(R+; Li(R x T)) 

|Vw|2eLi(R+;Li(RxT)). 

Since a > 1/2, the first bound implies in particular that 

erje L°°(R+ x R x T), 

and in particular if £q is small enough (compared to the reference height which is 1 here) , then 1 + srj 
is bounded from below. We infer that u is bounded in L°°(R+; ^^(r x T)) x L2(R+, iJi(R x T)). 

It is standard (see [2], [23]) to deduce that 1 + erj is compact in _L^^^(R+; H^(Il x T)), and that u is 
compact in Lf^^(Ii+;L'^{'R x T)). 

Taking the limit in the non linear terms is now possible : we need indeed to deal with the following 

nonlinear terms : 

u-Wu, W- ((1 + erjju) and Au. 

1 + srj 

The compactness of 1 + e?7 and u derived above allows to deal with the two first terms in a standard 

fashion. For the last one we just have to recall that 1/(1 + er/) is bounded in i°°(R+; _ff^"(R x T)). 

This completes the proof of Theorem 8. D 

4.3.4.3. Convergence. — In this section our aim is to show that the capillarity term enables us to 
get rid of the defect measure present in the conclusion of Theorem 7 page 78. As the proof is very 
similar to that theorem, up to the compactness of 77^, we will not give the full details. The result is 
the following. 

Theorem 9 (strong convergence in the case with capillarity). — Under the assumptions of 
Theorem 8, denote by [rj^^u^) a solution of (4.3.10) with initial data (rfl,u^^), and define 

Up to the extraction of a subsequence, $e converges weakly in Lj^^ill^ ; Hi^^(Ti x T)) (for all s < Q) 
toward some solution $ of the following limiting filtered system: for all iX in & and for all smooth $^ 
in Ker(L — iXId), 

I ^■^l{x)dx-v I I A'i^<P-^l{t',x)dxdt' + I I QL{^,<P)-^l{t',x)dxdt' ^ j <P° ■ ^l{x) dx, 

where $" = (ry*',^*'). 
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Proof. — We will follow the lines of the proof of Theorem 7. In particular all the results of Sec- 
tions 4.3.1 and 4.3.2 arc true in this situation and wc will not detail the proofs. So the point, as in 
Section 4.3.3, consists in taking the limit as s goes to zero, of the equation on IlA^e. 

Equation (4.3.6) page 83 can be written here as follows: for all smooth $^ = ($^ g, ($^)') belonging 
to Ker(L-iA/rf), 



(4.3.11) 



exp ( ^ ) (r;e$Io + "e • ($D')(i, x)dx ~ / {v"M,o + 4 ' {'^'>d'){^)dx 



-en I /"(-A)" exp (—\ r^^V • (-A)"(l>^)'(t', x)dxdt' 

-/7tt^^('=^''(^)"')-*>'"'^"'^-'"' 

Remark 4-3.7. — We have chosen to keep the unknowns {rj^^U;,) and not write the analysis in terms 
of {rii;,m^) as previously (recall that nUg = (1 + erj^)u^): the study of m^ rather than Ug is indeed 

unnecessary here as the factor which appears in the diffusion term in the equation on u^ can 

1 + £% 
he controled in this situation, contrary to the previous case. The advantage of writing the equations 

on (rjg^Ui;) is that there is no nonlinear term in rj^, contrary to the previous study, hut of course the 

difficulty is transfered to the study of the diffusion operator; the gain of regularity in rj^ will appear 

here. 



Taking limits as £ ^ in the two first terms is immediate. For the third term, we simply recall that rj^ 
is bounded in L°°(R+; L^(R x T)) and £T]e is bounded in L°°(R+; 7J^"(R x T)), so e% goes strongly 
to zero in L°°(R+; iJ''(R x T)) for every s < 2a. Since <i>^ is smooth, we infer that 

EK f /"(-A)"exp( — UeV-(-A)"($^)'(t',a;)da;di'^0, ase^O. 
Let us now consider the fourth term, 

It is here that the presence of capillarity enables us to get a better control. Let us write 
= z/ /" /"v (exp (—\ u,\ : V {^D' {t' , x)dxdt' 

'•'[h ("- (?) "') ^ ' (ifk'*^'') "' "="'^'"' 

Clearly the first term on the right-hand side converges towards the expected limit: we have 



V 



f fv (cxp(—juA ■.V{^iy{t',x)dxdt' ^ly I /" V$';, : V($^)'(t', a;)dxdi', as£^0. 
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To study the second one, we can notice that 



^^^ (',T>*^' 1 _ V7 ( ^^^ \ l'^^*\' I ^^^ V7/',T>*^' 



(*a)' = V — ^ ($1)' + — ^V(<i>l)', 



1 + £% / \ 1 + e?7e / 1 + £?7, 

and since the second term on the right-hand side is obviously easier to study than the first one, let us 
concentrate on the first term. We have 

l + er]e 1 + e?7e (1 + e%)2 
Since er]^ is bounded in L°°(R+; iJ^"(R x T)), we infer easily, by product laws in Sobolev spaces, that 

e'^rjeVrie is bounded in L°° (R+ ; iJ'" (R x T)), for some cr > 0. 
But on the other hand rj^ is bounded in L°°(R+; iy^(R x T)), so we have also 

e^ry^Vr^e^O in L°°(R+; ij2"-2(R x T)). 
By interpolation we gather that 

e^-qe^-ne ^0 in L°°(R+; L2(R x T)), 
and the lower bound on 1 + erj;, ensures that 

^'''"^''" in L~(R+;L2(RxT)). 



(1 + £77^)2 



The argument is similar (and easier) for the term , so we can conclude that 

1 + e?7e 



t r V / /-•+' 



A (exp ( — j uA ■ {^iy{t',x)dxdt' ^v j /"v$A : V {^*^)' {f , x)dxdt' . 



10 J ^ + £Ve 
Finally we are left with the nonlinear terms: let us study the limit of 

/ /exp ^^ j {u, ■ V)ue • {^iy{t',x)dxdt' - I /exp (^ j r],u,\/ ■ ^l^{t' ,x)dxdt' . 

The study is very similar to the case studied in Section 4.3.3 (see Proposition 4.3.6), so we will not 
give all the details but merely point out the differences. First, one can truncate the integral in xi G R 
to xi G [— -R, R], where i? is a parameter to be chosen large enough in the end. As previously that is 
simply due to the decay of the eigenvectors of L at infinity. So we are reduced to the study of 



J[-fl,fl]xT 



exp ( j {ue ■ V)ue ■ {<^*x)'{t' ,x)dxdt' and 

/ / cxpf^^ \'q^Ue-V^l^{t',x)dxdt'. 

Jo J\-R,R]>cT V g / 



/O J[-B„B]xT 

The limit of the first term is obtained in an identical way to Section 4.3.3, since Ug satisfies the same 
bounds, so we have 

/* / exp /— ) {u, ■ V)u, • ($l)'(i', x)dxdt' -^ f [ Y. (K- V)*A • (^A)'(i', ^)dxdt' , 

Jo Jl-R.,R]xT V £ / Jo J ^+^=;, 

as e goes to and R goes to infinity. 
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Now let US concentrate on the last nonlinear term. With the notation of Section 4.3.3, we can write 

/ exp ( ^^ — ) T]eUe ■ V$^ o{t', x)dxdt' 

J[-R,R]kT \ ^ / 

/ exp I 1 (?7e - rie.N)ue ■ V$^ Q{t',x)dxdt' 

J[~R.,R]xT \ ^ / 

/ exp I I rje.Niue -Us.m) ■V^*xoit',x)dxdt' 

J[-R,R]xT \ ^ / 

/ exp I ri^.NUe,M ■ Vl>;^ o(*'> x)dxdt' . 

J[-R,R]y.T \ ^ / 

The first two terms on the right-hand side converge to zero, due to the following estimates: for 
some -1 < s < and for aU t e [0, T], 

/ exp ( 1 (?7e - ■qs,N)ue ■ V$^ o(*'' x)dxdt' 

J[-fl,fl]xT V ^ / 

< CrWrie - '7e,Af||L°=([0,T];H=([-i?,fl]xT))||Me||L2([o,T];Hi([-it:,fl]xT))|l*i>Allw2,oo(R_xT), 

and similarly 

/ exp I I ri^.N{ue -Us,m) ■ ^7^1^(1' ,x)dxdt' 

J[-R,R]xT \ ^ / 

< CT,N\\Ue - Ue.M\\L-^{[Q,T]:H'{[-R.,R]xT))\\^*\\\w^.^i'R.xT)- 

Finally the limit of the third term is obtained by the (by now) classical nonstationary phase theorem, 
namely we find, exactly as in the proof of Proposition 4.3.6, that 



* f /it'X 

/ exp ) ris,NUe,M ■ V$;^ o(i', x)dxdt 

J[-R,R]xT \ £ 



as e goes to and M, N and R go to infinity. 
That concludes the proof of the theorem. 



Jo J ...r.-x 



^i-|-^^=A 



D 
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A', diffusion operator in the Saint- Venant system, p. 18 
A'^, diffusion operator in tlie limit filtered system, p. 19 

£g(t), modulated energy, p. 73 

£e,N{t), modulated energy applied to a sequence of approximate solutions, p. 75 

T , Fourier transform, p. 5 

JF27 Fourier transform with respect to X2, p. 9 

/(n, fc), coefficients of / in the Hermite-Fourier basis, p. 9 

H'^ , inhomogeneous Sobolev space, p. 5 

W , homogeneous Sobolev space, p. 5 

H^{^), for s > 0, Sobolev space on a bounded set Q, with Dirichlet boundary conditions, p. 5 

H^''{n), for s > 0, dual space of -ffo(f^), P- 5 

HI, weighted Sobolev space, p. 30 

Jn, spectral truncation, p. 67 

Kn, truncation operator in the {'^n,k,j) basis, p. 40 
K, a regularizing kernel, p. 60 

£P-°°, weakly convergent series, p. 35 

L, singular perturbation, p. 7 

£, semi-group generated by L, p. 18 

Ns, pseudo-differential operator on KerL, p. 44 
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n„,fcj, projection on '^n,k,j, P- 17 

IIa, projection on Ker(L — iXId), p. 17 

Ho, projection on KerL, p. 8 

IIj^, projection on (KerL)-'-, p. 31 

Hk, projection on Kelvin modes, p. 16 

IIm, projection on mixed Rossby-Poincarc modes, p. 16 

lip, projection on Poincarc modes, p. 16 

Hr, projection on Rossby modes, p. 16 

Q, quadratic operator in the Saint- Venant system, p. 18 

Ql, quadratic operator in the limit filtered system, p. 19 

Ql, purely ageostrophic quadratic operator in the limit filtered system, p. 36 

5, set of indices {n,k,j), p. 30 

S* , set of indices {n,k,j) corresponding to (Keri)^, p. 35 
(5We), the shallow water system, p. 57 
(SWq), the limit system after filtering, p. 57 

6, spectrum of L, p. 16 

&K, subset of & corresponding to Kelvin modes, p. 16 
&p, subset of & corresponding to Poincarc modes, p. 16 
6p, subset of 6 corresponding to Rossby modes, p. 16 
&N, subset of 6 defined by a frequency truncation, p. 68 

i'''n,k.j, eigenvalues of i, p. 9 

W^'°°, Sobolev space, p. 12 

$0j first coordinate of the three component vector field $, p. 18 

$' ~ ($1, $2), two last coordinates of the three component vector field <i>, p. 18 

$' ~ ($2, —^1), image of <&' by a rotation of angle 7r/2, p. 3 

$, complex conjugate of $, p. 26 

<i>g = £ (— -) (%,Ue), where {rj^^u^) solves the Saint- Venant system, p. 78 

<I>, a solution to the limit system {SWo), p. 78 

$A, an element of Ker(L — iXId), p. 80 

<i>Ar, an approximation of $, p. 67 
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<I>(7V)7 an approximate Leray solution of (SWq), p. 40 

^e.N, an approximation of $g, p. 85 

^e.\,N, the projection of ^^.n onto Ker(L — iXId), p. 85 

$, a solution to the gcostrophic equation, p. 43 

ip , coefficients of $ G Keri in the (^n,o,o) basis, p. 43 

4>Ni a corrector to $7v, P- 69 

'fn,k,j, coefficients of $ in the (^n.fe.j) basis, p. 15 

■0„, Hermite functions, p. 8 
^Pn.fcj, eigenvectors of L, p. 10-11 



